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New measurements of the equilibrium in the exchange reaction between methanol and 
heavy water, in conjunction with data in the literature, are considered in relation to the 
vibrations of the hydroxyl group, including the torsional degree of freedom. An interpretation 
is presented which allows for a torsional frequency of 250 cm~, a value which is shown to be in 
reasonable agreement with the entropy from thermal data. 


HROUGH the courtesy of Professor E. F. 
Barker of the Physics Department, new 
data on the infra-red absorption of methanol 
and methanol-d have been made available in 
advance of publication. In addition to frequencies 
which can be ascribed to motions of the methyl 


_ and methoxy groups, they include the OH and 


OD valence vibrations at 3683 and 2720 cm— 
and the distortion of the C—O—D angle at 863, 
but the corresponding distortion of the C-O—H 


angle unaccountably fails to- appear. The tor- 


sional motion produces a broad complicated 
band system which is spread over a large part 
of the low frequency region. 

The data are not adequate to permit a calcu- 
lation of the equilibrium constant of the reaction, 


CH;,;OH+DOH =CH;,0D+HOH, 


by the statistical method. However, the constant 
can be obtained experimentally to a good degree 
of accuracy, and it is the purpose of this paper to 
determine its value and apply it toward the 
search for the missing deformation and torsional 
frequencies. The exchange constant turns out to 

Insensitive to variations of the torsional 


frequency. For this reason, little specific informa- 
tion about the torsional motion is obtained, but 
the deformation frequency is consequently quite 
closely determined. 


EXPERIMENTAL 


A recent measurement by Okamoto! gives the 
distribution ratio of deuterium between methanol 
and water as 1.00+0.01 at 0°C, corresponding 
to an equilibrium constant for the exchange 
reaction of 0.500.005. This result was obtained 
from two sets of experiments, in which Hz or D2 
was brought to equilibrium with water or 
methanol with the aid of platinum black, so 
that no use was made of water-alcohol solutions, 
and consequently any peculiar activity effects 
which such solutions might introduce were 
eliminated. 

In our measurements, equilibrium solutions 
were made by mixing equimolar quantities of 
methanol and dilute deuterium oxide, and the 
alcohol and water were separated by fractional 


1G, Okamoto, Sci. Papers Inst. Phys. Chem. Research 
(Tokyo), 31, 211 (1936). 
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TABLE I. Methanol exchange equilibrium at 80°C. 


CH;0D| DOH 
DOH | DOH | Dirrer- | MOLES WATER |CH;OH| HOH 
INITIAL} FINAL ENCE MOLEs ALC. | K100 | 
2.644) 1.798) 0.846 0.982 0.848 | 1.811) 0.468 
7.610} 5.086} 2.524 0.990 2.575 | 5.194} .496 
7.610) 5.116} 2.494 0.985 2.529 | 5.225) .484 
9.844! 6.762} 3.122 1.027 3.322} 6.955} .478 
17.44 |11.90 | 5.54 1.005 S91 [12:52 | 
17.44 |11.86 | 5.58 0.990 §.85 |12.47 | .469 
24.90 |16.92 7.98 1.012 8.79 {18.20 | .483 
24.90 |16.74 8.16 1.005 8.94 |17.99 | .496 


Mean value, =0.481+0.015. 


distillation, using a 90 cm vacuum-jacketed 
fractionating column. The equilibrium constant 
was calculated from the decrease in density of 
the water, as shown by the balance temperature 
of a glass diver. Initial DO concentrations from 
1.3 percent to 12.5 percent were employed. 
The conversion of balance temperatures to D.O 
percentages assumes that there is no significant 
volume change upon mixing D.O and H.O, and 
that the density difference between the two kinds 
of water is constant in the working temperature 
range. The second assumption might introduce 
an error which would appear as a trend in the 
resulting constant, such that the equilibrium 
quotient would be a function of the total 
deuterium content. No trend was observed, and 
the error introduced is therefore within the 
limits imposed by the total experimental error. 
The temperature measurements upon individual 
samples of water were reproducible to within 
+0.02°, corresponding to an uncertainty of 
+0.006 percent D,O at the highest concentra- 
tions. The spread of the results is somewhat 
larger than this figure would allow, probably 
because of differences in the completeness of 
separation of the water from the alcohol, and 
because of fractionation of the water during the 
distillation so that the sample analyzed was not 
quite representative. 

The tact that the separation was not iso- 
thermal was not serious because the heat of 
reaction is very small. The result is assumed to 
be correct at 80°C. In the calculation. of the 
equilibrium quotient, the ratio of HzO to DOH 
was obtained for each concentration from the 
constant K = 3.37 for the reaction, 


D,O+H,0=2D0OH. 


O. HALFORD AND B. 
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The ratio is not sensitive to the selected value of 
this constant, but since the constant is a property 
of the vapor equilibrium, it should be remem. 
bered that its application to the liquid phase 
reaction assumes that the partial pressure of 
DOH is the geometrical mean of the partial 
pressures which D,O and H:O would have at 
the same concentration. 

The results are summarized in Table I, which 
shows, in successive columns, the initial and final 
total deuterium content (expressed as DOH 
percent of the water content), the decrease 
produced by reaction, the ratio of total moles of 
water to total moles of methanol, the ratio of 
DOH to H.O, and the equilibrium quotient 
K(J)353 for the liquid phase reaction at 80°C. 

The relation between the measured quantity 
K(l) and the vapor phase constant K(g) must 
now be considered, since only the latter can be 
calculated from available data. The transforma- 
tion would require that K(/) be multiplied by a 


factor 


p°(HOH) X °(CH;OD) 
p*(DOH) X p°(CH;,OH) 


in which the symbol p° refers to the vapor 
pressure of the pure substance, in order to give 
the first approximation to K(g). From the vapor 
pressure of the three kinds of water published by 
Wahl and Urey,’ the correction factor to be 
applied for the HOH/DOH ratio would increase 
the constant by about 3 percent. It is probable 
that this correction would be at least partially 
compensated by the corresponding effect with 
the methanols. Thus the vapor phase quotient! 
should be slightly greater than the measured 
quotient K(l). However, the principal error, 
incomplete separation of the methanol and water, 
would operate to give a high result. It is known 
from the spread of the results that the error 0! 
measurement is of roughly the same magnitude 
as the estimated difference between the liquil 
and vapor quotients. It appears, therefore, that 
the mean value of K(J) from Table I is probably 
very close to the correct vapor phase constant. — 
At 0°C, the temperature of Okamoto’s meas 
urements, the HOH/DOH ratio is about 11 


2M. H. Wahl and H. C. Urey, J. Chem. Phys. 3, 4!! F@ 
(1935). 
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If it is assumed that the corresponding ratio for 
the methanols is nearer to 1.0, say 1.05, the 
quotient for the vapor equilibrium at the ice 
point should be 0.525, which differs from the 
result for the liquid state by about five times the 
estimated experimental error. 


The statistical constant 


The calculation of the equilibrium constant of 
the exchange reaction from the masses, moments 
of inertia, and fundamental vibration frequencies 
by the well-established statistical method has 
been considered in detail by a number of 
authors.2 The paper by Urey and Greiff gives 
the fundamental equation except for the question 
of the character of the rotational function for 
molecules with more than three rotational de- 
grees of freedom. This has been treated rigorously 
by Eidinoff and Aston,* who have worked with 
certain special cases, and has been generalized 
by Kassel.° In addition, the function has been 


derived for ethane by Mayer, Brunauer and 


Mayer,’ and approximate methods of calculation 


have been employed by Halford’ and by Wirtz.® 


vs. 3, 411 


The equilibrium quotient is given by the 
expression, 


_L(ale.) _f(CHsOD)xf(HOH) 
L(aq.) f(CH;OH) Xf(DOH)’ 


K(g) 


in which f is the quantity frequently called the 
partition function and L may conveniently be 
designated as a partition ratio, such that 


Mp}? 
| 
Mu Ru 


p(1 


—e-hevilkT) 


*H. C. Urey and D. Rittenberg, J. Chem. Phys. 1, 137 
(1933). R. H. Crist and G. A. Dalin, J. Chem. Phys. 1, 
677 (1933). D. Rittenberg, W. Bleakey and H. C. Urey, 
J. Chem. Phys. 2, 48 (1934). H. C. Urey and L. J. Greiff, 
J. Am. Chem. Soc. 57, 321 (1935). K. Wirtz, Z. physik 
Chem. B31, 309 (1936); B34, 121 (1936). P. A. Small, 
a Faraday Soc. 33, 820 (1937). T. Jones and A. 

herman, J. Chem. Phys. 5, 375 (1937). 

M. L. Ejdinoff and J. G. Aston, J. Chem. Phys. 3, 


379 (1935), 


iL. S. Kassel, J. Chem. Phys. 4, 276 (1936). 
Ch J. E. Mayer, S, Brunauer and M. G. Mayer, J. Am. 
rem. Soc. 55, 37 (1933), 
a 0. Halford, J. Chem. Phys. 2, 694 (1934). 
K. Wirtz, Zeits. f. physik. Chemie B34, 121 (1936). 


TABLE II. Methanol moments of inertia 
(approximate, X 10%). 


AXIS CH;0OH | CH;0D 
Transverse, in COH plane 33.16 | 34.59 
Transverse, perpendicular to COH plane | 34.55 | 37.27 
Methyl group 5.268 | 5.268 
Hydroxyl group, axis through its C.G., 
parallel to CO bond 1.344 | 2.537 
Product = R’ x 0.811 | 1.723 


M is the molecular weight, R is Kassel’s de- 
terminant which reduces to the product of the 
principal moments of inertia if there is no 
internal rotation, and h, c, k and T have their 
usual significance. The subscripts H and D refer 
to the protium and deuterium compounds, re- 
spectively. The corresponding ratio for water 
contains a symmetry factor of two and in this 
case R is perhaps more conveniently expressed in 
terms of the principal moments of inertia. 


The water partition ratio 


The fundamental frequencies of the H2O and 
DOH molecules have been taken from a paper by 
Barker and Sleator.® For the ratio of the moment 
products R it has seemed safer to use the result 
of a calculation which assumes that the valence 
angle and the interatomic distances are the same 
in the two kinds of water. The values of L(aq.) 
at 0°C and 80°C are, respectively, 70.176 and 
35.807. These numbers are somewhat sensitive to 
the values of the constants / and k, but the 
effect of the sensitivity upon the equilibrium 
quotient is negligible. The constants used here 
have been taken from a table in Pauling and 
Wilson.!° 


Rotational functions of the methanols 


In the evaluation of Kassel’s determinant for 
the methanols, the C—O—H angle has been 
taken as 105° and the interatomic distances as 
C—H=1.09, C—O=1.42 and O—H=0.96 ang- 
strom units. The values of R for methanol 
and methanol-d, respectively, are 0.7989 and 
1.6671(X10-'*), and the square root of their 
ratio, which enters directly into the exchange 
quotient, is 1.444,. 


9E. F. Barker and W. W. Sleator, J. Chem. Phys. 3, 
660 (1935). 

10 Pauling and Wilson, Introduction to Quantum Me- 
chanics, McGraw-Hill Book Company (1935), p. 439. 
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It becomes desirable to compare this result 
with one which may be reached by an approxi- 
mate method, as a partial justification of the 
ratio used further on in the calculations con- 
cerned with the possibility of torsional vibration. 
In addition to the three principal moments of 
inertia which the methanol molecule would 
possess as a rigid body, a reduced internal 
moment is assigned, such as could be ascribed 
to a molecule like CH;CDs, equal to the ratio of 
the product to the sum of the moments of two 
coaxial rotators. Since the sum is equal to one of 
the principal moments, the quantity calculated 
as substitute for R in the exact equation becomes 
the quadruple product of the remaining principal 
moments and the two internal moments. Table II 
shows the results of the calculations. 

The square root of the ratio of the products 
becomes 1.457, about one percent higher than 
the correct value. Consequently, no significant 
error would be introduced into the exchange 
quotient by the use of these approximate 
moments. If the entropy of methanol-d were 
calculated with the moments of Table II, it 
would be high by 0.07 unit at 25°C. 


The hydroxyl deformation frequency 


If free rotation is assumed, the expression for 
the equilibrium quotient may be put into the 
form 


K(g)-L(aq.) = V, 
where 2.303 log V= (he/2kT) ivi. 


The factors (1—e~"*”/*7) are omitted because for 
all the frequencies with which this calculation is 
concerned they are negligibly different from 
unity. The one unknown frequency is evaluated 
by these equations directly from the equilibrium 
quotient. Numerical solution yields V353= 11.382 
from the new data presented here and V273 
= 23.188 from Okamoto’s data. The same result, 


1199 within a wave number, is obtained for }>»; 
4 


from both values of V, but the close agreement 
is obviously fortuitous, in view of the probable 
magnitude of the vapor pressure corrections 
associated with the low temperature measure- 
ment. The algebraic sum of the known methanol 
frequencies, including a correction for the shift 
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of v; from 1034 to 1040, is then introduced to 
give for the missing deformation frequency 
1105 

This is in very close agreement with the value 
selected from the Raman spectrum by Bates, 
Anderson and Halford," and independently by 
Mizushima, Morino and Okamoto,” at 1110. 

The spread of the experimental results of 
Table I corresponds to an _ uncertainty of 
+15 cm-'. Comparison with the observed OD 
deformation at 863 in conjunction with known 
isotope shifts for other hydrogen frequencies 
suggests that the upper limit of 1120 is the most 
reasonable result. Ratios of protium to deuterium 
frequencies are: limiting value 1.414, oxygen- 
hydrogen valence vibration 1.354, carbon-hydro- 
gen valence vibration 1.34, carbon-hydrogen 
deformation from methane 1.308, oxygen-hydro- 
gen deformation as above 1.280+0.018. The 
higher value is further indicated by the frequent 
observation that Raman lines for liquids are a 
few wave numbers below the correct frequencies. 

The above calculation suggests practically 
unrestricted internal rotation in the methanol 
molecule. In order to allow for deviation from 
free rotation, it would be necessary to adopt 
one of the lower values of the deformation 
frequency or to increase the equilibrium quotient. 
However, as will be shown below, a precise inde- 
pendent deformation frequency and a more 
accurate study of the equilibrium than is feasible 
at present would be required before this con- 
clusion could be confidently drawn. In order to 
illustrate this point, an approximation to the 
equilibrium constant will be derived from the 
entropy, which is itself different from the free 
rotation value by about one entropy unit in the 
ordinary temperature range. 


The entropy of methanol 


The entropy of methanol, as a hypothetical 
perfect gas at atmospheric pressure and 25°C, 
based upon free rotation and the new set o 
fundamental frequencies, is 57.72, as against the 
value 58.38 recently published by Kassel.” 


J. R. Bates, L. C. Anderson and J. O. Halford, J. 
Chem. Phys. 4, 535 (1936). J. O. Halford, L. C. Anderson 
and G. H. Kissin, J. Chem. Phys. 5, 927 (1937). 

12 Mizushima, Morino and Okamoto, Bull. Chem. So: 
Japan 11, 699 (1936). 

13 1,, S. Kassel, J. Chem. Phys. 4, 493 (1936). 
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The difference arises primarily from the new 
frequencies. The third law value obtained by 
combining the result of Kelley,“ S(L)o93= 30.26 
+0.2, with the heat of vaporization measured by 
Fiock, Ginnings and Holton,'® and the vapor 
pressure from International Critical Tables'® is 
56.6;+0.2. The excess above 55.72 is the con- 
tribution of the torsional motion. The harmonic 
frequency required to bring the entropy to the 
upper experimental limit of 56.83 is 370 cm™. 
If any reasonable smooth curve of entropy 
against frequency is drawn in such a way as to 
approach the harmonic oscillator curve at high 
frequencies and the free rotation contribution of 
2.00 entropy units at zero frequency, it will 
indicate that a frequency not far from 250 cm 
in the lowest vibrational state corresponds to the 
harmonic value of 370. 


The case of torsional vibration 


The statistical expression for the equilibrium 
constant is rewritten with the contribution of 
the torsional motion expressed separately as W, 
so that 


K(g)-L(aq.) V- W. 


Only three rotational degrees of freedom are 
now included in each moment product R. The 
factor W is obtained from the torsional fre- 
quencies, y for CH;OH and y’ for CH;0OD, by 
means of two equations, 


(1 
and y=1.265y", 


in the second of which the quantity y/y’ is set 
equal to the ratio of the square roots of the 
reduced internal moments of inertia from Table 
II. Although W is here expressed by means of 
a function derived from the harmonic oscillator, 
it appears to depend primarily on the frequency 
of the ground state. This follows from the fact 
that it approaches y/y’, the ratio for free rota- 


Kelley, J. Am. Chem. Soc. 51, 181 (1929). 
. Fiock, Ginnings and Holton, Bur. Standards J. 
esearch 6, 881 (1931). 
International Critical Tables, vol. 3, p. 216. 
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tion, as the frequency approaches zero, and 
becomes equal to e’‘v-¥/2*? for high frequencies, 
in much the same way as a theoretically correct 
expression in terms of the actual frequency 
would be expected to behave. 
For this reason the computation of W is based 
on the value y= 250 estimated from the entropy. 
The introduction of the Raman frequency 1110 
for the OH deformation gives V3;3= 11.493 and 
Ve73= 23.480. For the exchange constant, then, 
based on the above interpretation of the entropy 
and the function W, K(g)353=0.495 and K(g)o73 
= 0.522. The first of these lies within the spread 
of the experimental values and the second comes 
within the range of the vapor pressure correction. 
It is evident that the exchange constant would 
be a good measure of the torsional frequency 
and the deviation from free rotation only: if it 
were known to three significant figures. Although 
the accuracy could be improved, it is difficult 
to see how the required degree of precision could 
be attained. Agreement between the entropy and 
the exchange reaction occurs only in the region 
near the upper limit of the experimental values 
of both quantities, and then only when the 
deformation frequency is not too close to the 
most probable value. This study suggests, but 
does not prove, that either the exchange con- 
stant or the entropy is low. 
A minimum potential barrier of 2500 cal. may 
be estimated for methyl alcohol from the entropy 
by means of the tables recently published by 
Pitzer.!? This is nearly as high as the estimate of 
3000 for ethane, and appears to be too high, 
relatively, for the number of hydrogen atoms 
involved in the effect. It would be unwise, how- 
ever, to conclude from this comparison that the 
entropy is low, because the underlying factors 
are not well enough understood. In this con- 
nection, the surprisingly high potential barriers 
for ethyl alcohol and isopropyl] alcohol, 9000 and 
6000 cal., which Schumann! and Aston have 
shown to be consistent with equilibrium data 
and experimental entropies, are particularly 
interesting. 
17K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 


18S. C. Schumann and J. G. Aston, J. Am. Chem. Soc. 
60, 985 (1938). 
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The Absorption Spectra of Cyclopropane, Ethylene Oxide, Ethylene Sulfide, and 
Ethylene Imine in the Photographic Infra-Red* 


EuGENE H. EystTEr 
University of Michigan, Ann Arbor, Michigan 
(Received July 18, 1938) 


The absorption spectra of cyclopropane, ethylene oxide, ethylene sulfide, and ethylene imine 
have been investigated in the region between 7000 and 12,000A and the stronger bands have 
been assigned as third harmonics of the fundamental C—H stretching frequencies of these 
molecules. Rotational fine structure was not observed under high dispersion, but the well-defined 
band envelopes of the three asymmetric rotators provide interesting examples of these little 


investigated band types. 


INTRODUCTION 


LTHOUGH the majority of compounds of 
chemical interest are asymmetric rotators, 
comparatively little is known about the types of 
band envelopes which may be expected in the 
absorption spectra of any molecules of large 
asymmetry. The only molecule of considerable 
asymmetry whose spectrum has been investi- 
gated in great detail is water. This is a spindle- 
shaped rotator, and unfortunately the simplicity 
of the molecule prevents the occurrence of one 
of the three elementary band types. Still less is 
known about the spectrum of the asymmetric 
disk-shaped molecule, and hence it has been 
thought of interest to investigate a number of 
disk-shaped molecules of considerable asym- 
metry, and of sufficient complexity to insure the 
appearance of all of the elementary band types. 
The ethylene oxide and ethylene sulfide molecules 
were chosen because their moments of inertia are 
not too large, and their C—H frequencies are 
high enough to bring their third harmonics well 
within the range of the present Eastman I-Z 
plates. Ethylene imine was included in order to 
observe the envelope of the hybrid third N—H 
harmonic. Cyclopropane, the related symmetric 
rotator, whose fundamentals have been carefully 
investigated, was also studied to assist in making 
the vibrational assignments. 


EXPERIMENTAL 


The cyclopropane employed in these experi- 
ments was a pure commercial product prepared 


* Contribution from the Gates and Crellin Laboratories 
of Chemistry, California Institute of Technology, No. 659. 


by the Ohio Chemical Company for use as an 
anaesthetic, and was used directly from the 
ampoule. Photographs were taken at pressures 
ranging from 120 mm to 1 atmosphere, the gas 
being confined in a 6-meter absorption tube. 

Commercial ethylene oxide, obtained from the 
Ohio Chemical Company, was also used directly 
from the tank, photographs being taken at 
pressures ranging from 120 mm to 1.5 atmos- 
pheres in the 6-meter tube, and at a pressure of 1 
atmosphere in the 72-foot absorption tube of this 
laboratory. 

The ethylene sulfide was prepared from ethyl- 
ene dithiocyanate kindly furnished by Dr. E. 
R. Buchman, according to the method of 
Delépine and Eschenbrenner.! The product was 
simply distilled, as it is decomposed by the 
ordinary drying agents, and photographed at a 
pressure of 200 mm in the 6-meter tube. 

The ethylene imine sample was prepared from 
redistilled monoethanolamine by the method of 
Wenker,® and the product, twice distilled and 
dried over metallic sodium, was photographed 
at pressures ranging from 240 to 40 mm in the 
6-meter tube. A second sample, kindly prepared 
by Mr. Herbert Sargent, Jr. according to the 
same method, was photographed at a pressure 
of 1 atmosphere in a 3-meter absorption tube. 

Low dispersion spectrograms were first taken 
on a Bausch and Lomb glass Littrow spectro- 
graph, and the strong bands were then photo- 
graphed under the high dispersion of the 21-foot 
grating spectrograph. 


1 Delépine and Eschenbrenner, Bull. Soc. Chim. France, 
33, 703 (1923). 
2 Wenker, J. Am. Chem. Soc. 57, 2328 (1935). 
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INFRA-RED ABSORPTION OF 


VIBRATIONAL STRUCTURE 


From the standpoint of their fundamental 
frequencies, the molecules cyclopropane, ethylene 
oxide and sulfide, and ethylene imine form an 
interesting sequence. The four C—H stretching 
frequencies 7s, 1a(2), os, ¢a(2) of cyclopropane,’ 
go over to the four nondegenerate frequencies 
Ts, Ta, Os, Ga Of ethylene oxide as the symmetry 
is reduced from that of D3, to C2, and the 
number of H atoms from six to four. In ethylene 
imine, which undoubtedly has the symmetry C,, 
with the symmetry plane normal to the ring, 
the four C—H frequencies will remain rather 
similar to those of ethylene oxide, and may be 
characterized as 7’, 7a’, os’, Oc, Where the 7’ 
and o’ indicate, respectively, symmetry and anti- 
symmetry with respect to the symmetry plane 
of the CsH, group, normal to the actual plane 
of symmetry of the molecule. The N —H stretch- 
ing frequency, which may be called simply s, 
will be symmetric with respect to the symmetry 
plane. Table I shows the distribution of these 
frequencies among the irreducible representations 
of the appropriate point groups, and the infra-red 
activities of their fundamentals. It should be 
noted that Table I defines the location of the z 
axis normal to the plane of the three-membered 
ring, and of the x and y axes in the plane of the 
ring, with the former passing through the 
hetero-atom. The selection rules for the third 
harmonics and combinations of these frequencies, 
which are of greatest importance in this region 
of the spectrum, may be obtained by the usual 
group-theoretical methods.* This process yields 
the information that when only the above fre- 
quencies are considered, seven parallel and seven 
perpendicular bands are permitted for cyclo- 
TABLE I. Infra-red activity and irreducible representations of 


the high eee fundamental modes of C3He, 
Cy 40, and C:H;N. 


Du(CsHe) C2v(C2HsO, C2HaS) C,(C2HsN) 


Tey Ai, 
Tay Bi, M, 
Cs, Bo, M, 


Ga, A» 


*See King, Armstrong, and Harris, J. Am. Chem. Soc. 
58, 1580 (1936) for diagrams of symmetry coordinates 
corresponding to these cyclopropane frequencies. 

‘L. Tisza, Zeits. f. Physik 82, 48 (1933). 
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TABLE II. Observed bands, rotational structure types, 
and assignments. 


C2H.O 


8333 
8466 
8602.1 A or C 
8667.0 B, 37s 


8767.3 C 
8773.8 A, 3rq 
8815.8 B 
8835.8 C 
8881.3 A 


9021.7 C, 3e5 
9483 


C2HaS 


8637 A, 8642.3 B, 37,5’ 


8691.411 


8770.51, 
8806.01! 


8771.9 B, 37a 


8833.8 C 
8909 


8982 
9080.2 C, 3a 


8775.0 A, 3x4’ 


8921 
9065.3 C, 305’ 
9552.5 3s 


9006.81! 
9114.411, 


9704 
9807 
11333 complex 


11390 complex 
12420 4s 


11448 complex 11413 complex 


propane in the third harmonic region, while for 
ethylene oxide and sulfide fifteen bands are 
permitted, five of each of the three types A, B, 
and C. In ethylene imine, no less than twenty 
bands are permitted. While many of these bands 
will, of course, not be strong enough to be 
observed, the spectra do prove to be so complex 
that assignments can be made in only a few cases. 
In Table II are listed the observed frequencies, 
as well as a correlation of bands which appear 
to be analogous, and assignments of the active 
third harmonics in the four compounds. Wher- 
ever possible, the probable rotational structure 
types, A, B, or C, whose identification will be 
discussed in more detail in the next section, have 


also been indicated. 


BAND ENVELOPES 


As might well be expected, none of the ob- 
served bands showed a resolvable rotational 
structure, but the unresolved band envelopes do 
possess considerable significance in themselves. 
The envelopes of symmetric rotator bands have 
been calculated theoretically by Gerhard and 
Dennison,’ and the observed cyclopropane bands 
seem in satisfactory agreement with the theory, 
though overlapping of the 8770 and 8806 cm~ 
bands makes the interpretation of the micro- 
photometer traces of the perpendicular band 
rather difficult. The envelopes of the asymmetric 
rotator bands, however, are of the greatest 
interest. When the normal modes of vibration of 
an asymmetric rotator are such that the changes 
of electric moment are along single principal 


5 Gerhard and Dennison, Phys. Rev. 43, 147 (1933). 
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Fic. 1. Microphotometer traces of the three types of band envelopes of ethylene oxide. A trace of the water spectrum 
is included with the A and B type bands, which are contaminated with strong water lines. 


axes of inertia of the molecule, the resulting 
rotational structures may be said to be of pure 
types, and are conveniently characterized as of 
types A, B, or C, as the change of electric mo- 
ment is along the minor, intermediate, or major 
axis of the inertial ellipsoid. In molecules of 
symmetry Ci, Ci, C2, or C,, where the change of 
electric moment may have components along 
more than one axis, bands of a twofold (C2, Cs) 
or threefold (Ci, C;) hybrid type may appear, 
whose envelopes are linear combinations of the 
intensity distribution functions for the three 
related pure band types, with coefficients pro- 
portional to the squares of the components of 
the change of electric moment along the corre- 
sponding principal axes. In order to know which 
type of band should be associated with each of 
the C—H frequencies, it is necessary to know 
the location of the three different inertial axes in 
the molecule. For this purpose, approximate 
moments of inertia of ethylene oxide and 
ethylene sulfide were calculated on the basis of 
reasonable molecular models,® giving the follow- 
ing values, in g cm?X10-*: 
Ethylene Oxide Ethylene Sulfide 
=31.8 38.8 
I,=B=39.2 Iy=B= 74.8 
I,= C=59.0 I,= C=102.0 
The moments of inertia of ethylene imine are 
very nearly equal to those of ethylene oxide. 
These results, in conjunction with Table I, 
make possible the following correlation between 
the band types and the irreducible representa- 
tions of the normal modes: 
Ethylene Oxide A~B, BrA;, C~B, 
Ethylene Sulfide A~A,, BrB,, C~B, 
Ethylene Imine A~A", B, Chybrid~A’ 


6 Dr. H. A. Lévy has kindly informed me of preliminary 
results of his electron diffraction studies of the compounds. 


The A’ harmonics in ethylene imine may possess 
hybrid rotational structure, but this could hardly 
be observable for the C—H harmonics unless 
they also involve a considerable stretching of the 
N-—H bond. The N—H bands, on the other 
hand, would be expected to present a twofold 
hybrid structure involving the band types B 
and C. 

These predictions are confirmed in every re- 
spect by the observed band envelopes. The 
ethylene oxide bands at 8667, 8881, and 9022 
cm~'!, shown in Fig. 1, are definitely of three 
different types, of which two are indeed some- 
what similar in appearance but differ greatly in 
the relative intensities of the side and central 
branches. Since the 9022 cm band appears to 
be the analogue of the 9114 cm band of cyclo- 
propane, one may with confidence identify the 
former as of type C. The 8667 cm band, appears 
to be of type B, which will differ from the per- 
pendicular band of the symmetric disk-shaped 
rotator principally in the appearance of a mini- 
mum of intensity at the origin. The remaining 
band at 8881 cm~! may consequently be desig- 
nated as type A. This identification may seem 
a little less arbitrary after the publication of 
theoretical calculations of the three types of 
band envelopes, which are now in progress in 
this laboratory. The C—H bands of ethylene 
imine showed no evidence of hybrid structure, 
and were very similar in appearance to the 
corresponding ethylene oxide bands. The type 
band of ethylene sulfide showed scarcely 4 
trace of side branches, and even when micro- 
photometered seemed to be little more than a 
single intense central branch. The type A bands 
were weak, yet seemed rather similar in appeat- 
ance to the type A bands of ethylene oxide. 
However, the type B band, shown in Fig. 2, 
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demonstrates in a very striking manner the 
effects of the change in the moments of inertia 
from those of ethylene oxide. The band is 
obviously considerably distorted as a result of 
the difference in the moments of inertia in the 
two vibrational states, but the important feature 
is that the double maxima about the band 
center have here become more widely separated, 
and much less intense relative to the rest of 
the band. 

The ethylene imine band at 9553 cm! does 
indeed show the expected hybrid structure, as 
can be seen by comparing Fig. 3 with the corre- 
sponding pure band types, B and C, of ethylene 
oxide. The low frequency maximum seems a 
little too intense, but it may have been 
strengthened by an accidental superposition of a 
weak band of type C. This view is supported by 
the doubling observed in this maximum. The 
surprising feature of this band is the low intensity 
of the type C Q branch, which means that the 
contribution of the type C structure is very small. 
If it were feasible to determine intensities in this 
region, the contribution of the two components 
could be determined with some accuracy, since 
the intensity distribution functions for the two 
pure band types involved could then be deter- 
mined experimentally from the C—H bands of 
the very similar ethylene oxide molecule. If it 
were then possible to be certain that the change 
of electric moment took place along the line of 
the N—H bond, the angle between this bond and 
the plane of the ring could be calculated. The 
problem is a complicated one, but it surely 
appears that this angle is much smaller than the 


8750 8800 cm” 


Fic.2.M icrophotometer trace of the strong type B band 
of ethylene sulfide. A ‘trace of the water spectrum in this 
Tegion is included. 


9525 9575 cM 


Fic. 3. Microphotometer trace of the hybrid N—H band of 
ethylene imine. 


64° angle between one N—H bond in NH; and 
the plane of the other two bonds. The structure 
of N-methyl ethylene imine is being studied in 
these laboratories by Dr. H. A. Lévy, to see if 
the analogous angle in this compound can be 
determined with assurance. 

It is of interest to observe that in all these 
molecules analogues of the 8720A (11465 cm~) 
band of ethylene? appear with considerable 
intensity. These bands in cyclopropane and 
ethylene oxide were photographed under high 
dispersion, and in both cases were found to be 
unresolved bands of considerable complexity, 
apparently the result of the superposition of a 
number of simple bands. Definite assignments of 
the components is not possible, since even in 
cyclopropane, which has been most carefully 
studied, the assignments of the C—H bending 
frequencies does not seem completely satis- 
factory. It is, however, probable that in all cases 
these bands are the result of the superposition 
of combination bands involving three C—H 
stretching and two C—H bending frequencies. 

The other bands of Table II were too weak 
for high dispersion study, and many were visible 
only on the low dispersion plates. 

The author is glad, in conclusion, to éxpress 
his appreciation to Dr. E. R. Buchman and Mr. 
Herbert Sargent, Jr., of these laboratories, for 
their advice and assistance in the preparation of 
ethylene imine and ethylene sulfide, and to 
Professor R. M. Badger for much _ helpful 
assistance throughout the course of these in- 
vestigations. 


7 Badger and Binder, Phys. Rev. 38, 1442 (1931). 
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The Spectrum of Allene in the Photographic Infra-Red* 


EvuGENE H. EystTEerR 
University of Michigan, Ann Arbor, Michigan 
(Received July 18, 1938) 


The spectrum of allene between 7000 and 12,000A has been photographed under high dis- 
persion, the observed bands being assigned as harmonics and combinations of fundamental 
frequencies of the molecule. Rotational analysis of the parallel band at 11,444A has given 
97.0 X 10-* g cm? as the large moment of inertia. The perpendicular band at 11,017A has been 
found to be unexpectedly complicated, but was too weak for detailed analysis. It has been 
concluded that the dimensions in ethylene and allene are: 


C—H=1.087A, C=C=1.330A, H—-C—H=116°. 


INTRODUCTION 


LTHOUGH the rotational structure of the 
parallel type rotation-vibration bands of 
symmetric rotators is simple and well under- 
stood, that of the perpendicular bands may be 
complicated by strong rotation-vibration inter- 
actions which, although adequately understood 
in most cases, appear capable of producing com- 
plexities which, as in the 4.54 band of SiH,,' 
remain at present unexplained. Consequently it 
has seemed of interest to investigate the rota- 
tional structure of the bands of allene in the 
photographic infra-red under dispersion sufficient 
to resolve the fine structure of both the perpen- 
dicular and parallel types, comparing that of the 
former with the structure predicted by the first 
order theory of Teller,? and with the use of that 
of the latter to determine the large moment of 
inertia of the molecule. Since this moment is 
insensitive to reasonable changes in the carbon- 
hydrogen bond distance and the carbon-hydrogen 
bond angle, a knowledge of its value would de- 
termine with considerable accuracy the carbon- 
carbon double bond distance, a quantity of great 
interest in the theory of molecular structure. 


EXPERIMENTAL 


The allene employed in these investigations 
was a portion of a very pure sample sent from 
Harvard University by Dr. W. E. Vaughan to 
Professor Linus Pauling and Dr. L. O. Brockway 


* Contribution from the Gates and Crellin Laboratories 


of Chemistry, California Institute of Technology, No. 658. 
1 Steward and Nielsen, Phys. Rev. 47, 828 (1935). 
2E. Teller, Hand- und Jahrbuch der Chemischen Phystk 
Vol. 9, No. 2 (1934), p. 125. 


for use in electron diffraction studies. A second 
sample for photographs at higher pressures was 
kindly sent by Dr. J. W. Linnett of Harvard 
University. A Pyrex absorption tube six meters 
in length and having a volume of one liter gave 
adequate absorption with relatively small gas 
samples. Exploration photographs taken at a 
pressure of one atmosphere on a Bausch and 
Lomb glass Littrow spectrograph revealed a rich 
spectrum in the region between 7000 and 
12,000A, and all intense bands were then photo- 
graphed in the first order of the 21-foot grating 
spectrograph of this laboratory, at pressures 
ranging from 40 mm to 1.5 atmospheres. Wave- 
length measurements on the prism spectrograms 
were made with respect to sodium, potassium, 
barium, and calcium lines, and on the grating 
spectrograms with respect to standard iron lines 
in the second and third orders. Eastman I-Z and 
144-P plates, hypersensitized with ammonia, 
were used, and in every case visible radiation was 
removed with filters. 


VIBRATIONAL STRUCTURE 


In Fig. 1 are shown those three of the eleven 
allene symmetry coordinates which correspond 
to the modes of vibration most important in this 
region of the spectrum, namely those in which 
the motion consists chiefly of stretching of the 
C—H bonds, and also their distribution among 
the irreducible representations of the point group 
Dea, which describes the symmetry of the mole- 
cule. The fundamental frequencies may, for 
clarity, be associated with these symmetry ©0- 
ordinates, which will in this case approximate the 
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INFRA-RED SPECTRUM OF ALLENE 


true normal coordinates with considerable ac- 
curacy. The symmetry selection rules are given 
in the usual fashion* by the transformation 
properties of the electric moment and polariza- 
bility tensor components. The irreducible repre- 
sentations of these components and of the har- 
monic and combination states are given in 
Table I. Raman investigations of liquid allene by 
Kopper and Pongratz* give values for two of the 
C—H frequencies, »1=2995 and »;=3060 
Bonner and Hofstadter® report 3086 cm as the 
position of the unresolved fundamentals vz and 
v3. With the aid of these approximate frequen- 
cies, the analogous frequencies in ethylene, and 
the symmetry selection rules, the strong bands 
between 7000 and 12,000A may be assigned as 
harmonics and combinations of the three C—H 
stretching frequencies. In Table II are given the 
observed bands, their assignments, and the irre- 
ducible representations of the active components 
of the excited states, which are determinable 
from the observed band types. The unassigned 
bands, with the exception of the one at 8758A 
which will be mentioned later, were extremely 
weak and were observable only on the low dis- 
persion plates. The photographic data for allene, 
in conjunction with similar data for ethylene,® 
whose fundamental frequencies are well known, 
lead one to the conclusion that the probable 
values of the C—H stretching frequencies in 
gaseous allene are: 


vi = 3005, ve= 2981, v3= 3095. 


TABLE I. Selection rules for the group Dea. 


COMBINATIONS 
A, Aso B, 
A, Ao B, 


A, Bz 
A, 


HARMONICS 


J Ai, even 
[Ai n odd 


[B.]= {$e n even 


E3 
2Ai1+A2+Bi+B, 
M,+iM,~E 


Orr +Qyy, ly Ayzy 


sd L. Tisza, Zeits. f. Physik 82, 48 (1933). 

al and Pongratz, Weiner Berichte B141, 840 
’ Bonner and Hofstadter, Phys. Rev. 52, 248 (1937). 
‘Bonner, J. Am. Chem. Soc. 58, 34 (1936). 


Fic. 1. C-H stretching symmetry coordinates for allene 


ROTATIONAL STRUCTURE 


Two of the bands, 3v2 and 273+ 72, show rota- 
tional structure of the usual parallel type. The 
lines, however, are incompletely resolved even at 
a pressure of 40 mm, and cannot be detected on 
the microphotometer trace shown in Fig. 2, 
though on the original plates a number of them 
are distinctly visible. By reversing the plates on 
Velox No. 5 with contrast developer, about 
twenty lines in 3v2 could be measured directly 
under the comparator with respect to the iron 
comparison spectrum on the same print. The 
other band, however, was so badly overlapped by 


TABLE II. Assignment of bands in the photographic infra-red. 


FREQUENCY 
(cm7) ASSIGNMENT 


8,738.96 
11,394.0 


2vitve 
11,205 
11,139 2v3t+ve 
11,097 
11,017.2 
10,920 
9,602 
9 338 
8,758 
7,192 


WAVE-LENGTH 
(A) 


11,444.32 


A, 
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.B, 
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Fic. 2. Microphotometer trace of the 11,444 and 11,394A bands of allene; the lower 
trace shows the water lines in this region. The general trend in the background in the upper 
trace is not significant, and is absent in the lower trace only because the water plate was 


much more heavily exposed. 


the much stronger 32 as to be quite useless for a 


determination of the rotational constant. Al- 
though it was impossible to measure the first few 
lines in either branch, the relative numbering of 
the lines in the P and R branches was easily made 
because the convergence was found to be ex- 
tremely small. Furthermore, only one scheme of 
relative numbering produced a set of A’’s7’s 
which showed a reasonable dependence on J. 
The absolute numbering, determined by the 
position of the Q branch, may be in error by one 
unit because this branch is rather weak and wide, 
but the error thus induced in the value of the 
rotational constant will be insignificant owing to 
the small convergence. The broadness of the 
band lines, the presence of strong water lines, 
and the overlapping, with 271+ 72 on the high 
frequency side made it impossible to obtain long 
branches or measurements of extreme accuracy. 
The usual combinations, shown in Table III, 
give a value of 97.0X10~° g cm? for the large 
moment of inertia of the molecule. 

The perpendicular band, 3v3 at 11,017.2A, was 
found to be extremely weak, and increasing the 
pressure above one atmosphere in the six meter 
tube served only to broaden the lines. A number 


of broad, rather widely spaced band lines could 
be distinguished on the original plates, and were 
even more clearly evident in twofold enlarge- 
ments made on Velox No. 5. A microphotometer 
trace made with the microphotometer slightly 
out of focus to eliminate the effects of the coarse 
grain of the Eastman type J plates, is shown in 
Fig. 3. Careful comparison with the enlarged 
prints has shown that the maxima in the trace 
do indeed correspond to the band features ob- 
servable visually, but it must be remembered 


TABLE III. Rotational lines and combinations in 3v2 
at 11,4444. 


R(J) P(J)  A’2T(J)/2(2I +1) 


J 

5 8741.488 

6 42.046 

7 42.680 

8 43.330 9.574 
9 43.911 8733.106 — 
10 44.507 _— 11.934 
11 _ 31.977 13.117 
12 45.731 31.390 — 
13 46.274 30.841 15.484 
14 46.817 30.294 16.576 
15 29.698 17.637 
16 29.180 

17 28.677 

18 28.090 
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9050 9075 9100 


9125 9150 


Fic. 3. Microphotometer trace of the perpendicular band of allene at 11,017A, smoothed 
and corrected for the change in background. Suggested numberings of the Q branches are 
given. The pQ branch only extends to the fifth line because the band overlaps the 11,100A 


bands beyond that point. 


that many other weaker lines may have escaped 
detection. In spite of this fact, the rotational 
structure can be seen to deviate markedly from 
that predicted by Teller’s theory, in that the 
spacings are not uniform, many of the Q branches 
are doubled or tripled, and a number of extra 
strong lines have made their appearance. It is, 
however, of interest to note that the small 
moment of inertia, calculated from the average 
spacing of the lines near the center of the band 
with the constant ¢ in the Teller formula taken 
as zero, is 5.8 X 10-*°. One would expect the small 
moments of ethylene and allene to be very nearly 
equal, and this value does indeed agree very 
satisfactorily with the small moment 5.67 X 10-*°, 
determined by the fine structure of the 3.24 and 
the 11.54 bands of ethylene reported by Levin 
and Meyer.’ Hence one is led experimentally to 
the conclusion, which also follows from a con- 
sideration of symmetry coordinates in the repre- 
sentation E, that the value of ¢, which is a 
measure of the vibration-rotation interaction, is 
not large. A satisfactory explanation of the 
rotational structure of this band must quantita- 
tively account for the difference between the 
rotational levels associated with this 3v3 level of 
allene and those associated with the analogous 
3v,8 levels of the methyl halides, which appear 
to be in complete agreement with the Teller 
theory. The data are at present far too meager to 
permit a quantitative consideration of the prob- 
lem, but it seems appropriate to point out two 
alternative explanations of the complex rota- 


Levi 
16, 137 (1998) J. Opt. Soc. Am. and Rev. Sci. Inst. 


5In th : 
(1935), —_—" of H. Sponer, Molekiilspektren Vol. I 


tional structure of this band, and to suggest a 
possible method for distinguishing the correct one. 

It is easy to show* that whereas the 37; level 
of allene will be split into two degenerate E levels 
by the anharmonicity, the 3v, level in the methyl 
halides will be split into three levels, only one of 
which is of the type E and capable of giving per- 
pendicular bands in absorption from the ground 
state. The remaining twofold degeneracy of the 
E states will only be removed by terms in the 
Hamiltonian such as the vibration rotation inter- 
action, which do not possess the symmetry of the 
molecule. Since the magnitudes of these two types 
of perturbations which remove the degeneracy 
of the quadruply degenerate third harmonic state 
cannot be estimated with any certainty, the 
possibility that the rotational structure of the 
perpendicular allene band may be due to either 
of the following causes must be considered. 

1. In case the degeneracy of the E levels has 
not been removed, and if, moreover, they have 
not been widely separated by the anharmonicity, 
the 11,017A band may be the result of transitions 
to the two E states, whose rotational levels have 
been so perturbed by interaction with the vibra- 
tion that the resulting band cannot be simply 
decomposed into two ordinary perpendicular 
bands. In this case the difference between the 
allene and methyl halide third harmonics is easy 
to understand, and the fundamental perpen- 
dicular bands of the former might be expected 
to agree with the Teller theory. 

2. Alternatively, it is possible that the two E 
levels of allene have been widely separated, and 
that the 11,017A band is the result of a transition 
to only one of these states, whose degeneracy has 
been removed by a strong interaction with the 
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rotation, not present in the methyl halides, and 
not described by the first order Teller theory. In 
this case the fundamental perpendicular bands 
should show a similar structure. 

The writer proposes to investigate the funda- 
mentals of the degenerate frequencies of allene 
in the near future, in the hope of determining the 
true cause of this unusually complicated perpen- 
dicular band structure. 

The band at 11,100A, shown in Fig. 4, presents 
the appearance of a superposition of two parallel 
bands, and seems to be the result of the approxi- 
mate coincidence of the levels 2v3+»; and 
2v3+v2, which have the symmetry properties 
characterized, respectively, by the irreducible 
representations and Ai+Bi+B2. 
The only infra-red transitions to this composite 
state are to the two Bz components, which will 
produce two parallel bands. Fine structure 
similar to that observed in the 11,444A band 
could be faintly distinguished on the original 
plates, but not even the Q branches, which are, 
of course, rather weak in a molecule of this ex- 
treme spindle type, can be located with certainty 
on the microphotometer traces. 

The band at 8758A, shown in Fig. 5, is of 
special interest as the analog of the ethylene 
band at 8720A,° whose anomalous intensity and 


9 Badger and Binder, Phys. Rev. 38, 1442 (1931). 


EUGENE H. EYSTER 


complicated rotational structure have never been 
satisfactorily understood. The allene band ap.- 
pears to consist of a superposition of two paralle! 
bands of unequal intensity, and the components 
may be assigned as 3ve+2v4 and 
where the new fundamental frequencies are : 


1381 (Bo), v;=1031 


These frequencies, observed by Bonner and 
Hofstadter in the infra-red, are associated chiefly 
with the bending of C—H bonds. It is probable 
that the 8720A ethylene band is similarly a super- 
position of two or more combinations of C—H 
stretching and bending frequencies. 

The remaining bands were too weak for high 
dispersion investigation, but the band at 7192A 
seems to be 5ve, and its doublet structure is 
clearly visible on the low dispersion plates. 


STRUCTURES OF ALLENE AND ETHYLENE 


Because of the importance of the carbon- 
carbon double bond distance in establishing the 
single-double bond resonance curve,'® it seems o/ 
considerable interest to attempt to interpret the 
moments of inertia of ethylene and allene in 
terms of bond distances and angles. The C—H 


10 Pauling and Brockway, J. Am. Chem. Soc. 59, 1225 
(1937). 
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Fic. 4. Microphotometer trace of the 11,100A bands of allene. 
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Fic. 5. Microphotometer trace of the 8758A bands of allene. 


distance in these compounds can be obtained 
with some confidence by a consideration of the 
distances and force constants in a number of 
molecules containing the C—H group. When the 
C-H distances in C2H2, HCN, and the 
ground state of CH are plotted against the in- 
verse cube root of the C—H stretching force 
constant, following the procedure first suggested 
by Badger,!! one obtains a very satisfactory 
straight line. The force constants of ethylene and 
allene may then be used to fix the C—H distance, 
which is thus found to be 1.087A in both com- 
pounds. This information then makes it possible 
to calculate the C—H bond angle in ethylene 
from the small moment of inertia, giving the 


" Badger, J. Chem. Phys. 2, 128°(1934). 


angle H—C—H as 116°. Calculation of the C=C 
distance in ethylene and allene from their large 
moments of inertia, using the above parameters 
for both molecules, gives, respectively, 1.325 and 
1.330A. The difference is probably insignificant, 
since the large moment of inertia of ethylene, 
calculated by Badger’ from the rotational 
Raman band of Lewis and Houston, is probably 
not so accurate as the others, and none of the 
moments of inertia have been corrected to the 
vibrationless state. 

In conclusion, the writer is pleased to express 
his appreciation to Professor R. M. Badger for 
his valuable advice and assistance during the 
course of this investigation. 


1 Badger, Phys. Rev. 45, 642 (1934). 
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The Raman Spectrum of Liquid Diborane 
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The Raman spectrum of BeHg(/) is found to contain an isotopic triplet (v1) 793, 806, 821 


VOLUME 6 


cm~!, a broad line (v2?) 1108 cm and a very intense doublet (v3?) 2102, 2523 cm as well as 
several weaker lines. Abnormally large values for the B—B, C—C, and Si—Si bond force 
constants in BzH¢, C2He, and SisHs suggest that these bonds are compressed in these molecules 
by forces which may be associated with those which restrict internal rotation. The value of 
55.341.5 cal./deg. is calculated for the entropy of BzH¢(g) at 25°. Combined with the known 
heat of formation this value leads to a free energy of formation of —30,000+4000 cal. for 
B2He(g) at 25°. This large value does not truly indicate high stability, for diborane decomposes 


easily to form other boron hydrides. 


HE x-ray! and electron diffraction? patterns 

of diborane (BsH,) show that its atomic 
arrangement is analogous to that in ethane. 
However, the formulation of its electronic 
structure is made difficult by the deficiency of 
two electrons.’ In the hope of contributing 
toward the solution of this problem, we have 
obtained the Raman spectrum of diborane in the 
liquid condition. 


EXPERIMENTAL METHOD 


Preparation of the substance 


Diborane was prepared by the method of 
Schlesinger and Burg,’ and purified with great 
care. The vapor tensions of the first and last 
fractions corresponded to that previously re- 
corded for the pure substance: 225 mm at 
—111.8°. The quantity employed was slightly 
more than three liters of the gas at standard 
conditions. The sample was completely con- 
densed (by the aid of liquid nitrogen) in the 
Raman tube, which was sealed off and at once 
put to use. 


Optical arrangements 


The apparatus in which the diborane received 
and scattered the light from a mercury lamp, is 


1H. Mark and E. Pohland, Zeits. f. Krist. 62, 103 (1925). 

2S. H. Bauer, J. Am. Chem. Soc. 59, 1096 (1937). 

3 A review of the literature on this subject is given by 
Wiberg, Ber. 69B, 2816 (1936). 

‘Schlesinger and Burg, J. Am. Chem. Soc. 53, 4321 
(1931). 
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shown in Fig. 1. The Raman tube was made 
from a 12.5 cm length of standard wall, 10 mm 
diameter Pyrex tubing, sealed at the lower end 
to a Wood light trap of black Pyrex glass, and at 
the upper end to a wider (18 mm) tube carrying 
a horizontal plane window and the side tube for 
introducing and removing the sample. This 
Raman tube was surrounded by a 19X4 cm 
copper cylinder, out of whose wall (2 mm thick) 
a 13X3 cm window had been cut. The piece cut 
out was soldered to the inner wall to make a 
well for receiving liquid nitrogen, by whose aid 


{ 


DY 


MMMM 


Fic. 1. Schematic drawing of Raman tube assembly. The 
alcohol is kept at —100°C with liquid nitrogen to liquily 
the low boiling material in the Raman tube. 
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RAMAN SPECTRUM OF LIQUID DIBORANE 


the required low temperature (in the neighbor- 
hood of —100° or lower) could be maintained. 
The entire inner surface of the copper cylinder 
was silver plated to minimize the loss of light. 

The Raman tube and reflecting cylinder were 
immersed in an alcohol bath, contained in a 
cylindrical Dewar vessel which had been silvered 
only around two-thirds of its circumference. As 
a safety device, in case of an accidental rise of 
temperature, solid carbon dioxide was placed 
in the space below the window of the copper 
cylinder. 

The mercury vapor lamp used in this work 
was of the very convenient design shown in Fig. 
2. In the use of this relatively simple type of 
lamp, the current passes through the Torricellian 
vacuum in the water-jacketed tube. The arc is 
struck by forcing the mercury to the top of the 


TABLE I. Observed frequencies in the Raman spectrum of 
BH, (1) with displacements from the exciting Hg 
lines }3984(1), 4047 (k), \4077 (2), \4339(g), 
\4347.5(f), and \4358(e), d=diffuse. 


PROBABLE ALTERNATE 
(ARB.) (cm~!) (cm™) (cm7) 

9 22603.9 2101.5k 

3 22571.1 25271 

3 22414.2 2101.97 
10 22182.0 2523.4k 

5 22145.6 792.5e 

4 22131.9 806.2e 

3 22116.7 821.4e 

4 21993.2 2522.97 

2d 21930.0 2775.4k 1008.1e 
2 21757.7 1180.4e 2947.7k 
3 20933.1 2106g 

1 20898 2101f 
10 20836.6 2101.5e 

0 20517 2522¢ 

0 20478 2521f 

3 20449.6 2488.5e 4255.8k 
8 20416.2 2521.9e 

7 20155.3 2782.8e Hg 4960.3A 
2 20070.3 2867.8e 4635.1k 
4 19818.4 3119.7e 4887.0k 


discharge tube, where it makes contact with the 
mercury in the closed end, and then releasing the 
pressure. The discharge tube has an internal 
diameter of 7.5 mm and the connecting ‘‘barom- 
eter’’ tube has a bore of 2.5 mm. Connected to a 
source of 110-volt direct current, through a six- 
ohm resistance, this lamp draws a current of 11 
amperes, and produces mercury spectrum lines 


quite sufficiently sharp for accurate Raman 
studies. 


Fic. 2. Water-cooled mercury lamp. The arc is struck in the 
18 cm long Torricellian vacuum. 


This mercury lamp was placed against the 
window of the Dewar cylinder, and the two were 
wrapped with bright aluminum foil, in order to 
save as much of the light as was feasible. The 
scattered light was taken off through a right- 
angle prism into the slit of a three-prism 
Steinheil spectrograph, set for a dispersion of 
18.7A per mm at a wave-length of 4350A. The 
instrument was in sharp focus between 4350A 
and 5100A. Three photographs were obtained: 
the first was exposed half an hour, the second, 
two hours, and the third, five hours, on Eastman 
DC ortho plates. The spectrum of an iron arc was 
photographed upon each plate for comparison ; 
on the third plate the spectrum of the mercury 
lamp itself also was photographed, with an in- 
tensity of exposure almost the same as that of 
the scattered mercury spectrum. 


EXPERIMENTAL RESULTS 


Seven Raman lines were obtained with the 
half-hour exposure, sixteen with the two-hour 
exposure which had a low background, and 
sixteen or more with the five-hour exposure 
which gave a uniform background. In Table I are 
given the estimated intensity (I) and frequency 
(v) of each of the observed lines together with 
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Iron Arc 


- Raman Spectram 


Fic. 3. The Raman spectrum of B.H¢(/) and its microphotometer tracing. 


the displacements from the various exciting lines. 
In Fig. 3 the five-hour exposure is reproduced 
together with the iron comparison spectrum, the 
observed frequency shifts, and a microphotom- 
eter tracing of the Raman spectrum. 

The lines 2101 k,e and 2523 k,e are most 
intense and exceedingly sharp; these frequencies 
have also been scattered from the Hg i, f, and g 
and the Hg / and 7 lines, respectively. The pair 
792e and 806e are sharp and the former is con- 
siderably stronger than the latter. The line 821e 
is broad, however, and its integrated intensity is 
probably greater than that of 792e. 

As indicated in Table I, in a number of cases it 
is uncertain whether the exciting line is \4358 
or \4047. We feel, however, that \4358 is re- 
sponsible for most of these lines since if they 
were of higher frequency (i.e., scattered from 
4047) they should also be observed as scattered 
from 44358. The broad line at 21,930 cm may 
be an exception for if it were 2775k, 2775e would 
be obscured by the mercury line \4960. The 
fuzzy appearance of the line Hg 4960 causes us 


to favor this latter possibility. Note that the 
apparently sharp line marked x on the micro- 
photometer tracing is in reality only a dust 
mark on the plate. 


ASSIGNMENT OF FREQUENCIES 


Selection rules for this type of molecule have 
been given by Stitt and Yost.’ They allow 6 
frequencies for the D3a configuration (center of 
symmetry), 9 for the Ds,’ configuration (in 
which internal rotation is unrestricted) and 9 
Raman active frequencies for the D3; configura- 
tion (plane of symmetry perpendicular to the 
central bond). Howard® has worked out the 
normal coordinates using a three-constant va- 
lence-force potential function with restricted 
internal rotation. 

The symmetrical stretching vibration v1,° for 
each of the three isotopic molecules containing 


5 Stitt and Yost, J. Chem. Phys. 5, 90 (1937). 

6 J. B. Howard, J. Chem. Phys. 5, 442 (1937). We have 
used the notation used in this paper to designate the 
frequencies. We are indebted to Dr. Howard for sending us 
the unpublished formulae for the v9, v10, and »; vibrations. 
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RAMAN SPECTRUM OF LIQUID DIBORANE 


B'-—B", and B'’—B", are readily 
identified as the lines 793, 806, and 821 cm, 
respectively. The calculated isotopic shifts in 
frequency agree within 1 cm~ with the observed 
shifts, and the relative intensities are of an order 
of magnitude corresponding to the relative 
abundance, 4 to 1, of the boron isotopes except 
for the line 821 cm which is more intense than 
corresponds to this ratio. It has been noted, 
however, that the line at 821 cm™ is broad and 
therefore doubtless derives its principal intensity 
from another frequency, possibly v9. 

A comparison of the spectrum of BH. with 
those of the analogous compounds C2H¢. and 
SisHs aids in making further assignments. In 
Table II are listed our results together with 
those for CoH¢(/),7 CoHe(g),8 and SigH¢(g).* The 
frequency at 1180 cm~ apparently corresponds 
to the 1463 cm™ line for CoH¢ and the 933 cm 
band for SigeHs which have been assigned to ve, 
associated with the symmetrical bending of the 
B-—H, C—H, and Si—H bonds, respectively. 
Likewise the intense pair 2102, 2523 cm—' would 
correspond to the pair 2884, 2941 cm of 
C,H,(/) and the line 2163 in SigH¢(g) which have 
been assigned to v3, the frequency associated 
with the symmetrical stretching of the six bonds 
to the hydrogen atoms. These frequencies lie 
close to the B—H stretching frequency 2535 cm— 
in In the 57 split in the 
vibration is supposed to arise through coin- 
cidence and consequent resonance” of this mode 
with one of the overtones 2v;, or 2v9. It is dif- 
ficult at this time to imagine another explanation 
for the appearance of such an intense pair of 
lines in this region; yet a resonance splitting of 
421 cm! seems so abnormally large that we 
make this assignment only tentatively. 

The frequency 2489 cm—'! might be tentatively 
assigned to vi; in analogy with the assignment" 
of the line of CoH¢(g) at 2940 cm-; the latter 
assignment is confirmed by the appearance of 
the rotation lines associated with it which were 


*P. Daure, Ann. de physique 12, 375 (1929); Glockler 
and Renfrew, J. Chem. Phys. 6, 295 (1938). 

* Bhagavantam, Ind. J. Phys. 6, 596 (1931) ; Lewis and 
Houston, Phys. Rev. 44, 903 (1933). 

*B. L. Crawford, Jr. and J. T. Edsall, private communi- 
cation of unpublished results. 

‘°E. Fermi, Zeits. f. Physik 71, 250 (1931). 

" Bartholomé and Karweil, Zeits. f. physik. Chemie 
B39, 1 (1938). 
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observed by Lewis and Houston.’ The meaning 
of the remaining frequencies is uncertain; a 
knowledge of the infra-red frequencies would 
probably help in assigning them. 


THE ForRcE CONSTANTS 


With the assignment of the three parallel fre- 
quencies 71, v2, vs, and Howard’s formulae® for 
the normal modes of vibration derived from a> 
valence force potential function it is possible to 
calculate the force constants of the molecule. 
We obtain 3.010° dynes/cm for K; associated 
with stretching the B—H bond, 0.26105 
dynes/cm for Ke associated with bending the 
B-—H bond and 3.6X10° dynes/cm for Ks; 
associated with stretching the B—B bond and 
determined principally by 11. The two B—H 
force constants are in satisfactory agreement 
with the B—H force constants in the triborine 
triamine molecule B;N3H, in which the B—H 
stretching constant is 3.42105 dynes/cm and 
the B—H bending constants are 0.35 and 


TABLE II. Raman Frequencies of B2He, C2He, and SizHe. 


C2He(l)? C2He(g)§ SieHe(g)5 
792.5(5) 
806.2(4) 979(1) 
821.4(3) 994(10) 993(10) 435(7) 
1180.4(2) 1463(5d) 1344(0) 910 
1460 to }(2) 
955 
2733(2) 2744(2) 
2768(2) 2778(0) 
2102(10) 2884(10) 2899(10) 
2489(3) 2922(1d) 2940(1) 2163(10) 
2523(10) 2941(10) 2955(10) 
2775(3) 2963 (1d) 2969 
2868(2) 2979 
3120(2) 2989 
2999 
3009 


0.2310° dynes/cm for the planar and out-of- 
plane bending, respectively. 


INTRAMOLECULAR FoRCEs IN B2H,g, 
CoH,, AND SioHs 


A comparison of the internuclear distances for 
B.H, and analogous molecules as calculated from 
the force constants and Badger’s rule” with the 


%q, R. M. Badger, J. Chem. Phys. 2, 128 (1934); 
b, 3, 710 (1935;. 
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TABLE III. Comparison of internuclear distances calculated from force constants and those observed by electron diffraction, 


BURG 


Ks R R R(oss.) K1 r r r(OBS.) 
(DYNES/CM (CALC.) (OBs.) —R(CALC.) (DYNES/CM (CALC.) (OBs.) —r(CALC.) 
ComMPOUND X 10-5) (A) (A) (A) X 10-5) (A) (A) (A) 
BoHe, 3.37 1.48 1.86? 0.38 3.0 1.19 1.272 0.08 
CoH, 5.62 1.37 1.5533 0.18 4.79 1.06 1.0918 0.03 
1.7 2.13 2.324 0.19 — 


internuclear distances as measured by electron 
diffraction, leads to interesting indications con- 
cerning intramolecular forces. In Table III are 
given in column II the values of the force 
constants K3; for the B—B, C—C, and Si-—Si 
bonds in the compounds C2H¢, and SizH, 
listed in column I. The corresponding inter- 
nuclear distances R calculated from Badger’s 
rule are given in column III for comparison with 
the values from electron diffraction given in 
column IV; differences between observed and 
calculated values are given in column V. 
Columns VI, VII, VIII, and IX contain the 
values for the B—H and C—H force constants 
K,, the calculated distances 7, the observed dis- 
tances, and the differences between observed and 
calculated distances, respectively. 

The fact that the calculated values of R are 
much smaller than the observed values for all 
three compounds strongly suggests that the 
bonds between the central atoms are compressed, 
in analogy with Badger’s association’ of the 
apparent stretching of the bonds in tetrachlorides 
and hexafluorides by repulsion between halogen 
atoms. If the magnitude of the compressional 
forces indicated here depend upon the internal 
angle, internal rotation would be hindered. 

The results for the central atom-hydrogen 
distances indicate that these bonds, too, are 
under compression. An effective attraction 
between hydrogen atoms of one CH; group in 
ethane with the C—H bonding orbitals of the 
other group possibly arising through considera- 
tion of the terms suggested by Penney! might 
give rise to just such compressional forces. 

Any such attraction between hydrogen atoms 
would give the molecule a plane of symmetry 
perpendicular to the central bond (D3;,) rather 


13 ane and Brockway, J. Am. Chem. Soc. 59, 1223 
(1937). 

( — and Beach, J. Am. Chem. Soc. 60, 1836 
1938). 

% WW. G. Penney, Proc. Roy. Soc. Al44, 166 (1934). 


than a center of symmetry (D3a). In fact 
Bartholomé and Karweil! recently have ob- 
tained an infra-red spectrum for ethane which 
they could interpret using Ds, selection rules, 
but not with D3q selection rules. They concluded 
that ethane has a plane of symmetry perpen- 
dicular to the C—C bond rather than a center 
of symmetry. 

It is difficult at this time to make an estimate 
of the magnitude of the restricting potential," 
or the magnitude of the dependence of the com- 
pressional forces on the internal angle. 


THE STABILITY OF DIBORANE 


Since the heat of formation is known" and it 
is now possible to calculate an approximate 
value for the entropy, it seems desirable to 
discuss the stability of diborane at this time. 
From electron diffraction results? the moments 
of inertia are calculated to be A=B=55.2+2.0 
and C=14.3+0.5X10-* g cm?, on the assump- 
tion of tetrahedral distribution of bonds about 
boron atoms. The contribution of rotation and 
translation to the entropy of the gas at 25°C is 
thus 50.94+0.1 cal./deg. including the sym- 
metry term with the symmetry number equal 
to 18. The maximum contribution of internal 
rotation is 3.26 cal./deg. for no restricting poten- 
tial, Vo=0, while for Vo=3000 cal. the con- 
tribution would be about 2.0 cal./deg.'® A safe 
estimate of this contribution would probably be 
2.4+40.8 cal./deg. We have calculated the ap- 
proximate values for the missing frequencies 
using Howard’s formulae* and from them find 


16 American authors favor a value of 3000 cal. in the case 
of ethane (see Kistiakowski, Lacher, and Stitt, J. Chem. 
Phys. 6, 407 (1938) for references) while Eucken and co- 
workers estimate the height of the barrier to be about 350 
cal. (W. Hunsmann, Zeits. f. physik. Chemie B39, 23 
(1938)). 

17 Roth and Borger, Ber. 70, 48 (1937). 

18K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 
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the contribution of vibration to the entropy to 
be only 1.9+0.6 cal./deg., allowing a wide 
margin of uncertainty for the missing fre- 
quencies. We thus have the approximate value 
of 55.3+1.5 cal./deg. for the entropy of BsH¢(g) 
at 25°. 

The heat of formation of BeH¢(g) from the 
elements has recently been determined to be 
44,000+3000 cal.'7 We may make an estimate of 
the entropy of solid boron and write the following 
equations : 


2B(s)+3H2(g) = 
S* 98.1, H.2(g) =: 31.23 cal./deg., 
S*o98.1, B(s) =4-+1 cal./deg., 
S*98.1, BeHo(g) =55.341.5 cal./deg., 
AS° 998.1 = —50.443.5 cal./deg., 
AHo95= —44,000+ 3000 cal., 
AF°o9s = — 30,000+4000 cal. 


The uncertainty in the entropy of diborane is 
indeed rather large, but it is to be pointed out 
that this uncertainty can at most lead to an error 
of only 1000 cal. in the standard free energy 
while the uncertainty in the measured heat of 
formation is estimated!’ to be 3000 cal. The 
total of 4000 cal. uncertainty in the free energy 
cannot possibly invalidate the following con- 
clusions regarding the stability of diborane. 
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Namely, at room temperature diborane is 
stable with respect to decomposition into the 
elements and will not begin to decompose into 
the elements until about 500°C. Actually this 
process is not analytically complete short of 
700°. It is well known, however, that diborane is 
unstable at room temperature and above, de- 
composing to form other hydrides in a com- 
plicated set of reactions: 


25° 
+BsH.+---+H2(g) 


= white solid+yellow solid+yellow oil 
+B;Hot+BioHi4+He 


Side) 
B(s) +Ha(g). 


It will, of course, be impossible to relate the 
stability of diborane with respect to these 
products without a knowledge of their thermo- 
dynamic properties. 

We are indebted to Professor H. I. Schlesinger 
for the use of large amounts of diborane prepared 
with his material aid. We are also grateful to 
members of the staff of the California Institute 
of Technology for the use of equipment to 
evaluate the plates. 
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The Variation of the Raman Spectrum of Nitrate Solutions with Cation 


ELizABETH R. AND FRANCES M. PLEssEt* 
Department of Physics, Mount Holyoke College, South Hadley, Massachusetts 


(Received July 29, 1938) 


The strong Raman line from solutions of cadmium, silver and lead nitrates has been re- 
measured at varied concentrations, and the results of Grassmann confirmed, namely that the 
frequency shift diminishes with increasing concentration for lead and silver, and varies little 
for cadmium. A qualitative explanation of this behavior and that of sodium nitrate is sug- 
gested in the effect of the cations on the structure of water, and attention is drawn to the 
similar order of these nitrates in x-ray diffraction phenomena. 


HE Raman effect in solutions of nitrate has 
been studied by a number of observers, and 
from various standpoints. One of these has been 
the question of the dependence of the magnitude 
of the frequency shift on concentration, cation, 
and temperature. Gerlach! found for NaNO; and 
LiNQOs slight, barely measurable, increases in the 
frequency shift with concentration, but larger 
changes in various nitrate crystals. Grassmann? 
found variations in both directions, depending on 
the cation. Franklin® confirmed the increase in 
frequency shift for NaNO; solutions from 
1048.0+1 cm at 1.25 N to 1050 cm at 8 N 
and to 1068 cm in the crystal. On the other 
hand, a change of 12.5 cm™ reported by Em- 
birikos* for solutions of magnesium sulphate was 
not confirmed by Woodward and Horner® nor 
by Coon.* It seemed, hence, worth while to repeat 
some of Grassmann’s experiments, and for this 
the nitrates of silver, cadmium, and lead were 
chosen. 

The source of illumination was a glass mercury 
arc, the spectrograph a Hilger D78, the Raman 
tube of the jacketed horn form, with the horn 
and the circle around the window painted black. 
It was placed under the arc, and so that the light 
scattered from the farthest portion would just 
fill the collimator lens of the spectrograph with- 
out use of lens between the tube and slit. Water 
was ordinarily used in the jacket, but for some 
of the long exposures a solution of praseodymium 


* Skinner Fellow in Physics, 1935-36. 
1 Gerlach, Ann. d. Physik 5, 196 (1930). 
2 P. Grassmann, Zeits. f. Physik 77, 616 (1932). 
3 Franklin and Laird, Phys. Rev. 45, 738 (1934). 
* Embirikos, Zeits. f. Physik 65, 266 (1930). 
a — and Horner, Proc. Roy. Soc. 144, 129 
6 Coon and Laird, Phys. Rev. 47, 889 (1935). 


neodymium sulphate was used as filter. Times of 
exposure varied from three hours to thirty. The 
slit width was such as to give a slit image some- 
what narrower than the Raman line. An iron 
arc comparison spectrum was put on most of the 
plates at the end of the exposure. 
Measurements were made on a Hilger measur- 
ing micrometer. Earlier experience had shown 
that the use of a comparison iron arc in determin- 
ing wave-lengths was unsatisfactory, as shifts 
between this spectrum and the unmodified Hg 
lines sometimes appeared. These were suspected 
to be due to changes in temperature or baro- 
metric pressure during exposure. The room tem- 
perature was controlled by a thermostat, and 
thermograph records showed a variation of less 
than 2°F, but there was no way of controlling 
the pressure. While the effect of variations in 
pressure and temperature on the mercury lines 
would be to broaden them, but not to shift the 
relative center of gravity by a measurable 
amount, the iron spectrum put on in a few 
seconds might be measurably shifted, and by 
different amounts at different parts of a plate. 
Hence it was decided to use the mercury lines 
only, with the Hartmann formula, in computing, 
and to make a correction curve for the Hart- 
mann formula from some test plates where both 
the mercury and iron spectra were recorded with 
short exposures. Any slight error in determining 
this correction would affect all the measurements 
equally, and not alter comparative values. Dur- 
ing the progress of the work the paper of Teal 
and MacWood’ appeared giving data on the 
effect of changes in barometric pressure even 


7 Teal and Mac Wood, J. Chem. Phys. 3, 760 (1935). 
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over a period of three hours. The mercury lines 
4916.0A, 4339.23A, and either 4108.08A_ or 
4077.83A were used as standards. The deviations 
of the other mercury lines measured as checks 
were of the order of. 0.03A. 

The results reduced to vacuum are shown in 
Table I ; they have been computed from measure- 
ments on the strong Raman line excited by the 
Hg line 4358A. The corresponding line excited 
by 4047A was also measured but the accuracy 
was not so great. The dispersion in the neighbor- 
hood of 4566A is 25A/mm. The table of wave- 
number differences compiled by P. L. Bayley was 
used in checking. All but one of the silver nitrate 
plates were measured twice, and all but three of 
the lead plates, the } NV lead plates were measured 
four or more times on different days and at 
different points on the lines. The average was 
taken for each plate as the measurement of that 
plate, and the deviations shown are the devia- 
tions of these from the final average. In the case 
of lead nitrate there is also shown the average 
deviation of the separate measurements. The last 
column shows Grassmann’s values. 

It will be seen that in a general way the results 
are in agreement with those of Grassmann. His 
assumption that for silver nitrate the shift in- 
creases with dilution is confirmed, and it is shown 
that for cadmium nitrate the shift, if any, is very 
slight. For lead nitrate the variation found looks 
smaller than his, but the difference is within the 
sum of the estimated errors. Combining these 
results with those obtained earlier for sodium 
nitrate, one may arrange the four nitrates in the 
order sodium, cadmium, silver, lead, in that the 
effect of increasing concentration is with sodium 
to slightly increase the shift of the Raman line, 
with cadmium to change it very slightly if at all, 
and with silver and lead to diminish it, possibly 
more for the lead at the same concentration of 
the cation, also, the difference between the fre- 
quency shift in the crystal and in the dissolved 
State is larger for sodium nitrate than for 
cadmium nitrate, and for lead and silver is 
questionable. 

In the observations of Prins® on x-ray diffrac- 
tion by nitrate solutions a somewhat parallel 
order is found. Prins distinguishes the cases 


(193s J. Chem. Phys. 3, 72 (1935); Physica 3, 1016 


where a diffraction ring is found, the diameter of 
which increases with concentration, and inside of 
which, the diffraction diminishes as the angle 
approaches 0°, and the cases where the position 
of maximum diffraction does not change with 
concentration. He finds a mobile maximum in the 
case of cadmium nitrate but not for either lead 
or silver nitrates. He attributes the mobile ring 
to a kind of super-arrangement of the heavy ion, 
and the lack of it to an irregular gas-like distri- 
bution. Sodium is too light in comparison with 
the other atoms present for the method to be 
applied with success, and Prins found no small 
angle scattering of importance for sodium nitrate 
solution, but a pattern which he attributed to an 
arrangement that the water molecules assume in 
the electric field of the cations. As he notes that 
the phenomenon of super-arrangement decreases 
in the series K, Rb, Cs, one might infer that it 
would be large for sodium, although the experi- 
mental method does not show it directly. 

Prins attributed the lack of superarrangement 
in the case of lead and silver to nondissociated 
molecules. One finds no support for this in the 
electrolytic conductivities of these nitrates, and 
measurements of the intensity of the Raman line 
at different concentrations by several observers 
has shown that for sodium, cadmium, and lead 
nitrates it is proportional to the concentration of 
the NO; ion. It has been argued that this is not 
an evidence of dissociation as the NO; constituent 
would be present in either case. But as in all cases 
the frequency shift in the crystal is different from 
that in the dilute solution, there should also be a 
difference for light scattered from associated 


TABLE I. 
Nor- No. Av IN 
SUBSTANCE MALITY | PLATES Av IN CmM™! GRASSMANN 
4 5 1044.1+-0.4+0.6) 1043.8+0.6 
Pb(NOs)2 2.7 1046.0 
2 5 1045.7 +0.6+0.9 
1 2 1045 +1 
0.5 3 1046.7 +0.4+1.0} 1048.0+1.0 
Crystals 1 1044.3 1045 
AgNO; 7 2 1043.2 +0.1 1042.3+0.8 
3 4 1045.4+0.4 
1 2 1046.1+0.4 
Crystals 1045 (Gerlach) 
Cd(NOs)2 8 4 1047.4+0.3 1047.4+0.4 
2 3 1047.7 +0.3 
Cd(NOs)2°4H20 
Crystals 1 1052.2 
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molecules. The photometric curves of the Raman 
lines from the concentrated solutions appear how- 
ever those of single lines. If one concludes that 
dissociation is complete, one may still find in the 
ideas concerning the structure of water of Bernal 
and Fowler,’ used also by Prins, a qualitative 
explanation of the difference with respect to 
Raman scattering between the behavior of the 
nitrates of sodium and cadmium on the one hand 
and silver and lead on the other. The sodium and 
silver ions are singly charged, but the sodium 
being smaller exerts a greater influence in main- 
taining or increasing a certain regularity in 
structure in the water. The cadmium and lead 
ions are doubly charged, so that although the 
cadmium ion is as large as silver it produces an 
effect more like that of sodium, but the still 
larger lead ion again has a smaller directing 
effect, and its presence may lead to a weakening 


® Bernal and Fowler, J. Chem. Phys. 1, 515 (1933). 
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of the structure. This might be tested in a study 
of the water bands of small frequency shift 
referred to by Hibben’ as surely due to inter- 
molecular interaction. 

The NO; ion is considered to be plane, with 
the nitrogen at the center of an equilateral tri- 
angle, and the oxygen atoms at the corners. The 
Raman line studied here is ascribed to the radial 
vibration. The effect of the cations on this vibra- 
tion would be indirect according to the point of 
view suggested here. Where there is a super- 
arrangement of the cations, there are crystal-like 
groupings of molecules, and the symmetry makes 
for a higher vibration frequency, but if the water 
structure is weakened, and the cations are more at 
random, the central forces in the NO; ion may be 
weakened, and the frequency lowered. 

Our thanks are due to Dorothy Franklin, who 
assisted considerably in measuring the plates. 


10 Hibben, J. Chem. Phys. 5, 166 (1937). 


OCTOBER, 1938 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 6 


Temperature Scale and Thermionic Emission from Rhodium 
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When rhodium is heated to incandescence in an atmosphere of air or oxygen, a coating 
forms on its surface. This is shown by a change in the emissivity. When the temperature is 
increased sufficiently this coating suddenly disappears, and when the temperature is lowered it 
reappears. This coating is probably an activated adsorption of oxygen, or an oxide of rhodium. 
It is found experimentally that the temperature and partial oxygen pressure at which the 
coating disappears are related by the Clausius-Clapyron equation: 


(d log. P)/dT =L/RT?. 


P is the partial oxygen pressure, T the absolute temperature, and L the heat of vaporization 
of the surface. The electron work function of a clean rhodium surface is found to be ¢=4.80 


volts, and A =33 amp./cm?/degree’. 


EMISSIVITY AND TEMPERATURE SCALE 


HE peculiarities of rhodium surfaces de- 
scribed below were first encountered when 
measurements of spectral emissivity were being 
taken.! A specimen of rhodium in the shape of a V 
wedge 6 cm long, 4 mm wide on a side, and about 
0.04 mm thick was mounted inside a brass 


1L. V. Whitney, Phys. Rev. 48, 458 (1935). 


chamber which was evacuated to a pressure of 
1-2X10-* mm. The true temperature was 
measured by an optical pyrometer sighted inside 
the wedge, and the apparent temperature by 
another pyrometer sighted on the outside surface. 

A typical sequence of readings is shown in 
Fig. 1. When the temperature was first raised, 4 
series of values shown by the lower curve was 
found. Upon lowering the temperature there was 
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a sudden shift to the upper curve. The emissivity 
as calculated from points of the lower curve 
ranged from about 0.22 to 0.29 for \=6698A; 
that from points of the upper curve was about 
0.60. High emissivities were always found with 
this chamber at low temperatures after the 
temperature had once been raised beyond the 
transition point, the initial low emissivities at 
low temperatures never being found again after 
the first set of readings. The change from one 
curve to the other was always abrupt. If the 
temperature were being raised slowly, a critical 
point would be reached at which the emissivity 
would change and the true temperature flash up 
100 to 150 degrees. The sudden rise in tempera- 
ture was due to the fact that a low emissivity 
surface replaced the high emissivity surface. 
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It was difficult to take readings in the neighbor- 
hood of the critical point, and a shorter wedge, 
about 13 mm long, was therefore substituted. 
Because of the steep temperature gradient along 
this wedge, it was possible to have the center hot 
enough to show a low emissivity while the 
extremities showed a high emissivity. Figs. 2 to 5 
are photographs of this wedge at various 
temperatures. 

Figure 2 shows the outer surface of the wedge 
below the critical point, after it had previously 
been conditioned above it. The apparent temper- 
ature at the center is about 1500°K, and the 
surface has a high emissivity. Fig. 3 shows the 
wedge when the apparent temperature at the 
center is about 1790°K; this center part has low 
emissivity, while high emissivities obtain nearer 
the ends. Fig. 4 shows the appearance of the 
wedge when the temperature was slightly lower 
than it was in Fig. 3; the central low emissivity 
region is shorter. Fig. 5 shows the inside of the 
wedge when the outside had the appearance of 
Figs. 3 and 4. 

Figure 6 shows the same wedge after it had 
burned out. Five zones are distinguishable. The 
central zone had a blue appearance when cold, 
but at the time of burning out it had low 
emissivity characteristic of temperatures above 
the critical point. Zones adjacent to the center 
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had high emissivity, and the top and bottom 
zones low emissivity. (These last zones are not as 
clearly marked as the others and are not visible in 
Figs. 3 and 4.) These end zones had never 
become hot enough to pass the critical point and 
thus had kept their initial low emissivity through- 
out the experiment. X-ray pictures taken of the 
high and low emissivity surfaces showed no 
crystalline differences in the metal. 

Further investigation showed that prolonged 
outgassing affected the temperature of the critical 
point. This indicated a surface layer dependent 
on the temperature and pressure conditions. The 
brass chamber was therefore discarded and a 
study made of this effect in an all Pyrex system 
where better vacuum conditions could be ob- 
tained. At high vacua (2X10-* mm) the effect 
disappeared entirely, and the low emissivity 
0.242, corresponding presumably to the clean 
metal surface, was obtained at all temperatures. 
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The temperature scale for this emissivity and for 
a wave-length of 6698A is given in Table I. 

When air was introduced into the system the 
transition phenomenon immediately reappeared. 
This suggests very strongly that the high 
emissivity surface is covered by a surface layer 
which disappears at high temperature. Readings 
were taken of the temperatures of the critical 
point for various pressures of air, oxygen, 
nitrogen, and hydrogen. There was a marked 
difference between these results in Pyrex glass 
and those obtained in the brass chamber; the 
pressures necessary to obtain the transition in 
the all-Pyrex system were greater by a factor of 
10° or more. Perhaps the wax, which was used in 
sealing the brass chamber together, saturated the 
atmosphere about the hot filament and acted 
similar to a catalytic agent. When wax was 
purposely introduced into the Pyrex system the 
layer formed at lower pressures, but it is question- 
able whether conditions duplicated those inside 
the brass chamber. 

Figure 7 shows the transition temperatures 
obtained when gases were introduced into the 
Pyrex system. It is evident that the temperature 
of breakdown of what is apparently a surface 
layer is a function of the partial oxygen pressure, 
since the readings for air reduced to partial 
oxygen pressures fall on the curve found for pure 
oxygen. The readings for nitrogen indicate a 
small oxygen content, but the temperatures of 
the transition point were so low they could not be 
read accurately. No surface layer appeared when 
hydrogen was introduced into the system. 

Figure 8 is a graph of the logarithm of the 
oxygen pressure to the base e plotted against the 
reciprocal of the transition temperature. The 


TABLE I, 


TRUE TEMPERATURE APPARENT TEMPERATURE 


1200°K 1111°K 
1300 1197 
1400 1282 
1500 1365 
1600 1447 
1700 1528 
1800 1609 
1900 1689 
2000 1768 
2100 1847 
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equation of this line is 
log. P= —L/RT+48, 


and the slope gives the value of L/R. The heat of 
vaporization of the surface, (L), was found from 
this graph to be 6.15 X10‘ calories per g mole. 


THERMIONIC EMISSION 


A straight strip of rhodium 3 mm wide, 0.07 
mm thick, and 10 cm long was mounted in a 
conventional type two-guard ring tube such as 
that described by Wahlin and Reynolds.? The 
tube was evacuated and the glass baked for 700 
hours. During this baking the temperature of the 
filament was gradually raised to a value above 
the transition point. As soon as the furnaces were 
removed a thermionic determination was made. 
The Richardson plot showed a definite break at 
the transition point, the low temperature slope 
being the greater and corresponding to a work 
function higher by 3/10 volt. The high tempera- 


°H. B. Wahlin and J. A. Reynolds, Phys. Rev. 48, 751 
(1935). 
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ture portion of the plot gave a work function 
slightly greater than that of a clean surface. 
Continued heat treatment at a temperature 
above the transition point caused the break to 
disappear gradually and after 1100 hours the 
values ¢=4.80 volts, A=33 amp./cm* degrees* 
were obtained. At this time a final pressure of 
210-8 mm was reached. Additional heat treat- 
ment for 500 hours caused no further change. 
Two other specimens of rhodium gave the same 
values. 

This value of the work function is not in 
agreement with that reported by Dixon.’ He 
observed, however, an abnormal variation of the 
work function with temperature which can best 
be explained as due to incomplete outgassing. He 
also found definite changes in the thermionic 
properties of the metal at 1100°C, which he 
attributed to a change in the structure of the 
metal. In view of the x-ray findings reported 
here, such a conclusion seems untenable. 

This project was supported in part by a grant- 
in-aid from the Wisconsin Alumni Research 
Foundation. 


3 E. H. Dixon, Phys. Rev. 37, 60 (1931). 


a 
4 

for 
40- “a a 
red. 

high 
ings 20 
‘ical 

gen, 
‘ked 
lass 

d in 
the 

‘ted 
was 
the 

Ires 
the 
we 

| 
ace 
| 

tial 
ure 
of 

‘be 
en 
the 
the 


OCTOBER, 1938 


The diffusion of hydrogen through commercial and very 
pure anode nickel has been studied in a precision manner 
throughout the temperature interval 1100° to 150° Cen- 
tigrade. It was found that the diffusion rates vary in a 
perfectly regular manner throughout this region, except 
for a discontinuity at the Curie point, namely near 360°C. 
A very exact agreement with the formula R= AT =p» 
Xexp (—b/T) was found in all cases. After long heat 
treatment with hydrogen the value of y in the above 
equation for anode nickel becomes +0.50+0.01 with the 
possible exception of a region within 3° of the Curie point. 
All the nickel used had isotherms whose slope at the start 
of observations was a maximum in the Curie region and 
markedly greater than 0.5, but which approached 0.5 as 
the heat treatment, ordinarily considered as having a 
decarburizing effect, progressed. 
The magnetic transition as shown by diffusion is abrupt, 
the isobars below the Curie region being quite as straight 
_as above. There is always a sudden drop in the value of b as 
calculated from the simple equation R=A exp (—b/T) at 
the Curie point, the value changing in the purest nickel 
from 6600 to 5800 with rising temperature. No systematic 
breaks appear in the isobars at high temperatures, such as 
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found in pure iron, but nevertheless the value of z indicates 
that the diffusing hydrogen possesses appreciable energy 


_ within the lattice. 


The Curie region shows marked hysteresis, and is similar 
to the Curie region of pure iron in this respect. 

Existing theories of ferromagnetism seem to require 
modification, inasmuch as the Curie transition in nickel 
appears to take place about as abruptly as the body- 
centered-face-centered transition in iron. The loss of mag- 
netism at the Curie transition is essentially an atomic 
phenomenon and should not be associated primarily with 
conditions in aggregates of atoms or atomic domains. 
During the transformation for a finite time, there appears 
to be an additional electron in the valence shell of some 
nickel lattice points. 

A possible description of the process of the change from 
the ferromagnetic to the paramagnetic condition in the 
transition elements is that an electron is moved from the 
3d shell to one of the 4 levels, followed more or less quickly 
by a drop of another 4-level electron back to such a 3d 
level as to result in the pairing of electron spins in the 3d 
shell. 


HE diffusion of hydrogen through metals at 

high temperatures has been studied quite 
extensively during the past several years in the 
Pennsylvania State College laboratories. The 
data for iron have been published,! and in view of 
the anomalies found in the diffusion rates of 
hydrogen through iron, it was felt that a precision 
study of the diffusion of hydrogen through nickel 
would be justified. This study would have the 
double purpose of showing whether such a series 
of breaks in the diffusion isobars as found in iron 
exists in the isobars for hydrogen through nickel 
(incidentally this would also show whether the 
series was due to peculiarities of the apparatus), 
and also show the applicability of the steady- 
state method (I) to a precision study of hydrogen 
through nickel. Several other sets of data on the 
diffusion of hydrogen through nickel are available 
from the Pennsylvania State College laboratories, 


1C. B. Post and W. R. Ham, J. Chem. Phys. 5, 913 
(1937), hereafter referred to as (I). 
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namely Craninsky? and Ham.’ The technique has 
been improved from that used to obtain these 
two sets of data, especially the temperature 
control, therefore we may expect an improvement 
in the data for nickel. In this precision study by 
both the trapping and steady-state! methods, we 
have used exactly the same apparatus, with the 
exception of the metal disk used for diffusion. 
as in I. 

The diffusion rates were found to be very 
regular and to obey exactly the diffusion 
formula,‘ 


D=AP'T' éxp (—6/T) (1) 


where P is the pressure exerted by the gas on the 
ingoing side, J the temperature in absolute 
degrees, and b the so-called work function. A is 4 


2 Craninsky, Thesis for M.S. degree, Pennsylvania State 
College (1932). See also Fig. 4, this article. 

3W. R. Ham, Am. Soc. Metals, 25, 536-564 (1937). 
(1908 Richardson, Nicol, and Parnell, Phil. Mag. 8, 1-29 
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HYDROGEN DIFFUSION THROUGH NICKEL 


constant and D is the rate of hydrogen diffusion 
measured in cubic centimeters of hydrogen per 
square centimeter area per mm thickness at 
standard pressure and temperature. 

The Curie point in nickel is easily recognized 
when the steady-state method (I) is used as a 
precision measure of the diffusion rates. Moreover 
the hydrogen diffusion rates obey the same sort 
of formula below the Curie point as above, 
consequently we can interpret the offset at the 
Curie point as an indication of demagnetization. 
The fact that the law governing the diffusion 
rates is the same below the Curie point as above 
indicates that the disappearance of magnetism is 
a relatively abrupt phenomenon when no field is 
present in the specimen and the effect is brought 
about entirely by temperature effects. 

A marked hysteresis effect is noted at the 
Curie point when this transition is measured by 
hydrogen diffusion methods. This indicates that 
some of the transitions in the individual nickel 
atoms which are responsible for the loss of 
magnetism must be governed by some sort of 
“order-disorder’’ rearrangement as developed by 
Bragg and Williams® and extended by Fowler.® A 
possible mechanism for the demagnetization of 
nickel will be discussed in a later paragraph. 

Isotherms were obtained above and below the 
Curie point. The variation of the slope of the 
isotherms was found to be regular above and 
below this point, that is, for commercial nickel 
the slopes varied from a value of 0.5 about 150° 
above the Curie point to a value of 0.56 at the 
Curie point (370°C) and decreased to 0.5 again 
about 50° below this point. A similar variation 
was observed with pure anode nickel except that 
prolonged treatment of this nickel by the 
diffusing hydrogen gas gradually eliminated the 
variation in slope from 0.5, except in a possible 
region within 3° of the Curie point (359.2°C). 
Such a variation is difficult to explain in terms of 
Smithells and Ransley’s’ application of the 
Langmuir absorption isotherm to deviations of 
the pressure exponent from +3, but can be 
explained from an extension of Ham’s discussion 
of the effect of carbon in solution in iron. 


° Bragg and Willi Proc. . Soc. 151A, 540 (1935); 
145A, 699 (1934), Roy. Soc. 151A, 540 (1935) 


1 (1935), owler, Proc. Roy. Soc. 149A, 1 (1935); 151A, 


*Smithells and Ransley, Proc. Roy. Soc. 150A, 172 (1935). 
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EXPERIMENTAL PROCEDURE 


The apparatus used in this research has been 
previously reported in I, and is shown here in 
Fig. 1. Exactly the same procedure was used in 
studying the diffusion rates by use of the steady- 
state method as reported in I. 

As explained in I, the steady-state method is 
characterized by having a stopcock barely open 
to the pumps on the vacuum side of the diffusing 
disk. This enables a small back pressure to be 
built up between this stopcock and the diffusing 
disk, which can be conveniently read on a 
McLeod gauge. As shown in I, this back pressure 
is directly proportional to the rate as determined 
by the trapping method (see Fig. 2, I). We must 
also mention the fact that the success of this 
method depends to a large extent upon having a 
system of very fast pumps to maintain a high 
vacuum on the pumping side of the stopcock. 

The same temperature control was used as 
reported in I, vz., a photoelectric cell relay. 

The melting point of a pure silver bead 
furnished by the National Bureau of Standards 
was noted by means of a telescope, and the 
thermocouples standardized without disturbing 
their position. The silver bead holder was 
manipulated by magnets outside the apparatus. 


FFUSER 


pean 


| 


THERMOCOUPLE | 


Fic. 1. Apparatus for measuring gas diffusion 
through metals. 


Diffusers were prepared by welding the nickel 
diffusing disk to nickel tubes by atomic hydrogen 
welding as described in I. 

Isotherms were obtained by operating the 
pump on the hydrogen inlet side of the diffuser. 
By manipulating the opening to this pump and 
the opening into the diffuser any convenient 
partial pressure of the hydrogen can be obtained 
from atmospheric down to 0.5 cm Hg. The 
logarithm of the rate was plotted against the log 
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Fic. 2. Diffusion of hydrogen through commercial 
nickel. Sample XIII. 53 mil. AA, logio (steady-state rate) 
vs. 1/T; BB, logi (trapping rate) vs. 1/T; CC, logio 
(steady-state rate) —}logiT+1 vs. 1/T, DD, logi 
rate) —}logi T—ylogiP (in mm) +1 
vs. 


of the pressure and the slope taken to determine 
the value which we designate as the ‘‘slope of the 
isotherm.”’ 

Trapping methods were also used to establish 
the value of the diffusion rate in cubic centimeters 
of hydrogen at standard pressure and tempera- 
ture diffusing through the disk per square centi- 
meter area per millimeter thickness. These 
trapping rates also served to establish the 
proportionality factor between these rates and 
the rates as determined by the steady-state 
readings. 


EXPERIMENTAL RESULTS 


' Figures 2 and 3 show the diffusion rates of 
hydrogen through commercial nickel plotted as 


logio (steady-state rate) vs. 1/T, (2) 
logio (trapping rate) vs. 1/T, (3) 
logio (steady-state rate) —3 logio vs.1/T, (4) 


logio (steady-state rate) 
—} logio T—y logio P (in mm) +1 s.1/T. (5) 


Note that the successive corrections applied by 
subtracting }logi and logio 7+y logio P 
straighten out these curves in a surprising manner. 


Wy 


Fic. 3. Diffusion of hydrogen through commercial nickel. 
Sample XII. 38 mil. AA, logio (steady-state rate) vs. 1 7; 
BB, logo (trapping rate) vs. 1/T; CC, logio (steady-state 
rate) —}logio 7+1 vs. 1/T; DD, logio (steady-state rate) 
— logio T—y P (in mm) +1 vs. 1/T 


The various values of the constants used in 
several types of diffusion formulae are given in 
Table I for these two samples. 

Figure 5 shows the logarithm of the diffusion 
rate of hydrogen through anode nickel* plotted 
against the reciprocal of the absolute tempera- 
ture. Fig. 5 also shows the variation of the 
isotherm slope as a function of T before the 
sample has been completely purified by the 
diffusing hydrogen. Note that the maximum 
occurs at the Curie point within our experimental 
error. These slopes of the isotherms appear in the 
diffusion equation as P’, where y is the value of 
the isotherm slope as described in experimental 
procedure. Table II shows the constants in 
several types of diffusion equations for anode 
nickel. 

After several months running this sample o! 
anode nickel showed that the variation of the 
isotherm slopes from 0.5 was not as great as that 
shown in Fig. 5. It was found that the slopes 


* This anode nickel was kindly furnished by the Inter- 
national Nickel Co. According to information receiv 
from their laboratories this is the purest nickel they have 
been able to make up to the present time. Complete 
chemical analyses of this type of nickel are not available, 
but we have been assured that this metal is at leas 
99.99 percent pure. 
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HYDROGEN DIFFUSION THROUGH NICKEL 


remained 0.5 within 10° of the Curie point and 
then showed a slight variation to about 0.53 at 
the Curie point. 

A new sample of anode nickel 42 mils in 
thickness was installed and after several months 
running with hydrogen diffusing through the 
sample it was found that the deviation of the 
isotherms from the 0.5 law (i.e., R» P) ap- 
parently occurs only near the Curie point and 
there is no variation from 0.5 except possibly 
within 3° above this point. Fig. 6 shows the log of 
the rates plotted against 1/T for this sample in 
the neighborhood of the Curie point. A small 
variation is apparent in the region directly above 
the temperature 360°C. 

Avery careful attempt has been made from the 
data obtained from the last diffuser (42 mil) to 
find if the plot of log Rvs.1/T is really linear 
between 360° to 1060°C. It has been found that 
the plot is not linear and logio R—4} logio T must 
be plotted against 1/7 in order to obtain a linear 
relation. This is shown in Fig. 7. As mentioned 
before, this diffuser was in the system long 
enough so that the diffusion rates were completely 
reproducible and the variation of the isotherm 


TABLE 1. Diffusion constants for hydrogen through 
commercial nickel.* 


ABOVE CURIE BELOW CURIE 


SAMPLE No. XIII 


EQUATION 


D=kip) exp (—b/T) =6680 b1 =6780 
ki =0.93 X10 ki =1.02 X 1072 
D=k:P'T exp (—b/T) =6280 bs =6440 
ko =1.7 X10-4 ko =2.42 X 1074 
SAMPLE No. XII 
D=kP} exp (—b/T) bi =6650 bi =6590 
it Ki=2.0X1072 ki =1.9X 1072 
D=k:P'T} exp (—b/T) bz =6200 b2 =6200 
ke =4.0X 1074 ke =5.3 X1074 


*Inall of the above equations, D is the rate of diffusion of hydrogen 
measured in cubic centimeters at N.T.P. per square cm area per milli- 
meter thickness per second. P is the impressed pressure of the gas in 
millimeters of mercury. 


slopes from +0.5. were at a minimum. This 
proves experimentally the complete diffusion 
equation, Eq. (1), above the Curie point and 
doubtless below that point also, but the lower 
temperature range is not sufficient to determine 
the exponent of T with certainty. 

Figure 8 shows the experimental results of 
Craninsky? and Ham? obtained with relatively 
pure nickel obtained from Baker & Co. The 
existence of an offset at the Curie point in the 
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diffusion isobars is clearly indicated by the 
results of these investigators. 

The hysteresis at the Curie point for com- 
mercial nickel may be observed in Fig. 9. It 
should be noted that the rate obtained by the 
steady-state method enables one to follow this 
transformation under various conditions, and 
moreover, the method gives a clear insight into 
the behavior of this point when no magnetic field 
is present. The observations were taken for 
various speeds of heating through the transfor- 
mation and it should be emphasized that even 
with very slow rates of heating (or cooling) a 
hysteresis is observed which can only be mini- 
mized by slowly working back and forth through 
the region until the transformation point is 
closed in between two temperatures. If readings 
are taken too fast the lower phase can be carried 
into the upper region for about 100° while the 
upper phase can only be brought for a much 
shorter distance into the lower, say about 25°. 
This latter effect is difficult to study because of 
the slow cooling rates of the furnace at these 
temperatures. The hysteresis for pure anode 
nickel seems to be considerably smaller than that 
shown in Fig. 9. 

The possible effect of hydrogen concentration 
on the Curie point was studied by determining 
isobars at various pressures. No change in the 
position of the Curie point could be found for 
pressures from 1 atmosphere to 0.5 cm Hg. 
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Fic. 4. Variation of isotherm slopes for commercial nickel. 
Samples XIII and XII, 
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Fic. 5. Diffusion of hydrogen“through ‘anode nickel. 22 mil. AA, logio (steady-state rate) vs. 1/7; BB, logy (steady- 
rate) vs. 1/T for low temperatures. 
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DISCUSSION OF RESULTS 


The first fact to be noted from this research on 
the diffusion of hydrogen through nickel is that a 
series of anomalies in the diffusion isobars as 
found in iron! do not exist in nickel. The same 
apparatus and technique as used in I was used 
here, consequently we must assume that the 
Rydberg series found in iron is truly charac- 
teristic for hydrogen diffusion through iron. 
While the constants given in Table I for the 
diffusion of hydrogen through commercial nickel 
are quite accurate they vary considerably with 
the different samples, and it has been found that 
these constants vary considerably with the life 
history of a particular sample. 

The isotherms for hydrogen through nickel do 
not have a slope of 0.50 near the Curie point 
unless the nickel has been thoroughly decarbu- 
rized by hydrogen diffusion. Information from 
the International Nickel Company indicates that 
minute quantities of carbon were present in the 
anode nickel at the start. All evidence in the data 
for both iron and nickel is that the explanation 
given by Smithells and Ransley’ for a deviation 
from 0.5 of the slope for hydrogen diffusion 
isotherms for iron or nickel is completely errone- 
ous, their explanation being that the variation is 
due to surface conditions. On the other hand, 
there is nothing that contradicts the opinion 
expressed by Ham? that this deviation from 0.5 is 
due to the simultaneous application of two or 
more mass action effects, or in simple terms that 
some of the diffusing hydrogen may be associated 
with an alloying element (often carbon) in the 
metal in pairs or ions (or as He+). This also 
appears to be the case if the metal is alloyed 
with N or O. ‘ 

The proof of this is that isotherms for carbonyl 
irons are much above a 0.5 slope in the region of 
500°C until the carbon is removed whereupon 
they approach 0.50. In the case of pure nickel, 
isotherms that were considerably above 0.50 are 
completely reduced to 0.50 by a process recog- 


TABLE 2. Diffusion constants for hydrogen through 
anode nickel. 


EQUATION ABOVE CURIE BELow CuRIE 
D=kP} exp (—b/T) bi =5615 b1 =6600 
ki =3.3 X1072 ki =0.17 


N DIFFUSION OF HYDROGEN THROUGH 
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Fic. 6. Diffusion of hydrogen through nickel in the 
Curie region. 42 mil. Slopes based on Formula R 
=A exp (—)b/T). 


nized as a decarburizing process. The material 
upon which Smithells and Ransley base their 
theory to explain the deviation of the isotherms 
from the P} law had not been decarburized. In 
all the ranges of temperature they discuss in the 
case of iron and nickel we have found almost 
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Fic. 7. Relation between diffusion rates and temperature 
for hydrogen through pure anode nickel above the Curie 
point. 
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Fic. 8. Diffusion of hydrogen through Baker nickel, 
Craninsky? and Ham. 


perfect agreement in the variation of diffusion 
rate with the square root of pressure (P! law), 
provided the carbon is completely removed. 
Hence, the application of the Langmuir ab- 
sorption formula which Smithells and Ransley 
employed is not necessary or justified. 

The experimental fact that the exponent of T 
in the diffusion equation is +} is an important 
theoretical point in connection with the energy 
of the hydrogen in the gas phase and in the metal. 
From this the hydrogen atoms appear to possess 
very nearly the same energy within the metal as 
without, and if the hydrogen is tightly bound to 
lattice points during most of its residence within 
the metal, because of the proton’s light mass as 
compared to the metal atom of the lattice, but 
little energy of translation is allowed for the 
hydrogen. Its energy must be either vibrational 
or rotational or both. The reason that breaks 
indicating this rotation show in data on iron and 
do not in nickel may be because the latter has a 
relatively larger number of free electrons at any 
particular temperature as indicated by con- 
ductivity data, and the number of rotational 
states in which the hydrogen protons may exist 
at a particular temperature is correspondingly 
greater, thereby smoothing out the diffusion 
isobar. 

The abrupt appearance of the offset at the 
Curie point as shown by hydrogen diffusion 
studies is surprising in view of the existing 
theories of ferromagnetism. It was for this 
reason that the lower portions of the curves 


below the Curie point were observed with great 
care to ascertain whether a small bowing of the 
curves below the Curie point existed. If the 
electronic shifts necessary to cause demag- 
netization are governed by a statistical distri- 
bution, the fact that we are able to detect the 
transition at the Curie point would lead one to 
expect that we should be able to follow any trend 
or shifting due to temperature below the Curie 
point. Such a trend has not been found and if it 
does exist the magnitude is almost negligible 
because we have found the rates below this point 
to follow the same kind of diffusion law as the 
rates above. Evidentally the intrinsic magnetism 
of the nickel atoms must be governed by forces 
which break suddenly. The above arguments, of 
course, refer only to demagnetization caused by 
temperature effects in the absence of a magnetic 
field. 

To be specific, the theory of Slater® needs 
modification since this theory assumes a gradual 
rearrangement, as the metal is heated from 
absolute zero, of electrons in the third shell of the 
nickel atom, this rearrangement increasing 
rapidly as the Curie point is approached. The 
diffusion data indicate that such cannot be the 
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Fic. 9. Hysteresis at the Curie point in commercial nicke! 
determined by steady-state methods. 


case since the isobars are straight from 153° to 
360°C. 

A detailed picture of the nature of the Curie 
transformation is that certain atoms of the 
lattice must temporarily possess an additional 
valence electron during the transition, no doubt 
obtaining it from the 3d shell. This additional 
electron is necessary to interpret the behavior of 
the isotherms on pure anode nickel in the range 


8 J. C. Slater, Phys. Rev. 49, 537 (1936); 49, 931 (1936). 
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within 3° of the Curie point. However, this 
condition does not persist and evidently another 
valence electron drops back into the 3d shell and 
has such a spin that the electrons in the 3d level 
are now paired. When this has occurred the 
element may be considered as paramagnetic. 


CONCLUSIONS 


(a) The diffusion of hydrogen through nickel 
has been investigated from 1100° to 150°C with 
both commercial and pure anode nickel. Con- 
stants in the diffusion equations have been 
obtained. 

(b) The Rydberg series found in iron has not 
been found in nickel. The same apparatus as used 
in I was used here. 

(c) The isotherms are found to vary from 0.5 
in the neighborhood of the Curie point for pure 
anode nickel. When this nickel has been fully 
decarburized by exposure to diffusing hydrogen 
for long periods this variation from 0.5 only 
occurs within 3° of the Curie point. 

(d) The exponent of T in the diffusion equa- 
tion R= AT?P! exp (—b/T) has been found to be 
+} for the diffusion isobars from 360° to 
1050°C. 

(e) The magnetic transition of nickel as shown 
by diffusion is abrupt, the isobars below the 
Curie region being as straight as those above. 

(f) There is always a sudden drop in the value 
of B, as calculated from the simple equation 


R=A exp (—6/T), at the Curie point, the value 
changing in the purest nickel from 6600 to 5800 
with rising temperature. 

(g) The Curie region shows marked hysteresis 
and is similar to the Curie region of pure iron in 
this respect. 

(h) The explanation prepared by Slater to 
account for the effect of temperature on ferro- 
magnetism seems to require serious modification 
inasmuch as the Curie transition for nickel 
appears to take place about as abruptly as the 
body-centered-face-centered transition in iron. 

(i) A variation in hydrogen concentration in 
nickel near the Curie point does not change the 
position of the Curie point to a measurable 
amount. 

(j) The loss of magnetism at the Curie 
transition is essentially an atomic phenomenon 
and should not be associated with conditions in 
aggregates of atoms or atomic domains. 

(k) During the Curie transformation for a 
finite time there appears to be an additional 
electron in the valence shell of some nickel lattice 
points. 

(l) A possible description of the process of the 
change from the ferromagnetic to the para- 
magnetic condition in the transition elements is 
that an electron is moved from the 3d shell to one 
of the 4-levels followed more or less quickly by 
the drop of another 4-level electron back to such 
a 3d shell as to result in pairing of electron spins 
in the 3d shell. 
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An electron diffraction study of the oxide films formed 
on abraded iron specimens by heating in air at temperatures 
ranging from 80° to 260°C has failed to detect structural 
effects due to a ‘‘critical temperature’”’ in the neighborhood 
of 200°C. a-Fe.O; was found on specimens which were 
oxidized at any temperature above about 130°C provided 
the period of oxidation exceeded a minimum value which 
depends upon the temperature. Thus 1 hour at 180°, 
4 hr. at 165°, 13 hr. at 150° or 1000 hr. at 130° are just 
sufficient to produce an a-Fe,0; film whose diffraction 
pattern is well defined. Lower oxidation temperatures or 
shorter times tend to strengthen the diffraction pattern 
due to the primary oxide film which is assumed to be 
y-Fe2O3 (or Fe3;O4) and is detectible immediately after the 
initial abrasive treatment. The strengthening of the 
primary oxide pattern is not accompanied by appreciable 


grain growth, however, and traces of iron are present in 
the pattern after oxidation periods as long as 1700 hr. at 
120°C or below. Based on considerations of the presence 
of these iron rings it is estimated that the primary oxide 
film cannot be thicker than about 200A over some portions 
of the abraded surface. At 193°C and above the develop- 
ment of the first order blue interference color is charac- 
terized by the appearance of a new a-Fe,Q; film which is 
more coarsely crystalline and has a different type of pre- 
ferred orientation than the first-formed a-Fe,O; film. It is 
believed that this secondary film is formed by the outward 
diffusion of iron through the oxide layer. The axis of pre- 
ferred orientation of the a-Fe,O; is always symmetrical 
about the normal to the surface. This orientation cannot 
be attributed to pseudomorphism of the oxide crystals on 
the underlying iron as is the case with the primary oxide. 


INTRODUCTION 


HE nature and thickness of the air-formed 
oxide film formed on iron at various temper- 
atures has been investigated by a number of 
workers! during the last few years. Many of the 
results indicate that the film composition or the 
mechanism of film formation or both depend 
upon whether the oxidation is effected at temper- 
atures above or below 200°C. Vernon? found that 
when iron is heated in air at temperatures above 
200°C interference colors develop which are 
characteristic of the oxide film thickness whereas 
at lower temperatures the color sequence does not 
appear even under prolonged heating. From 
gravimetric measurements he concluded that 
oxide films formed below 200°C were much 
thicker than indicated by the surface color. In 
addition he found by gravimetric methods that 
the oxidation vs. time curve is a parabola above 
200°C but deviates therefrom at lower tempera- 
tures. Miley* confirmed the absence of interfer- 
ence colors below 200°C. Mehl and McCandless,‘ 


1 For recent reviews see: U. R. Evans, Metallic Corrosion, 
Passivity and Protection (London, 1937), Chapters II and 
rl Ae H. J. Vernon, Chemistry and Industry 56, 751 
(1 

2 W. H. J. Vernon, Trans. Faraday Soc. 31, 1668 (1935). 
a as} A. Miley, Carnegie Scholarship Mem. 25, 197 

*R. F. Mehl and E. L. McCandless, Trans. Am. Inst. 
Min. Met. Eng., Iron and Steel Div. 125, 531 (1937). 
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however, obtained interference colors on single 
iron crystals at temperatures as low as 140°C 
although their color vs. time curves at this 
temperature indicate that the color reaches a 
constant value which is unchanged on further 
heating. 

Several electron diffraction investigations have 
shown that oxide formed on iron from room 
temperature up to some rather poorly defined 
temperature is cubic with a lattice constant of 
about 8.4A. Unfortunately both Fe ;O, and 
y-Fe2O3 have a unit cell of this size and they 
are virtually indistinguishable by diffraction 
methods. Chemical considerations, however, have 
led Miley and others® to conclude that the cubic 
oxide is y-Fe2O3.6 Preston and Bircumshaw’ 
found the cubic oxide on specimens heated at 
135°C and 175°C but rhombohedral a-Fe203; was 
present when the oxidation temperature was 
above 200°C. Their tests were made on small 
pieces cut from larger sheets which Vernon had 
used in his gravimetric work and for this reason 
were not well suited for the purpose. 

5 |. litaka, S. Miyakeand T. limori, Nature 139, 156 (1937). 

6 The conclusion of others that this primary oxide 's 
y-Fe,O; has been accepted here for purposes of desig 
nating diffraction patterns but it must be emphasized that 
the experiments described herein do not distinguish 
between Fe;0, and y-Fe2O3;. 


7 Reported by Vernon, Trans. Faraday Soc. 31, 1677 
(1935). 
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OXIDATION 


Since these investigations indicate a critical 
temperature in the neighborhood of 200°C it 
seemed desirable to find out what electron 
diffraction methods could determine about films 
formed in air at temperatures above and below 
200°C. In some earlier work® it was demonstrated 
that the cubic crystals of the room-temperature 
air-formed film bear a definite orientation re- 
lationship with the underlying iron but no 
attempt was made to study films formed at 
higher temperatures. 


EXPERIMENTAL 


The specimens were small (3 
rectangular blocks cut from ingot iron (Armco) 
of the following analysis: 


Carbon 0.02 percent 
Manganese 0.03 
Phosphorus 0.019 
Sulphur 0.026 

Copper 0.05 


The initial surface treatment consisted of grinding 
all six surfaces on successively finer grades of 
emery paper ending with a very light abrasion on 
dry 4/0 paper. The resulting surface is in a finely 
abraded condition and gives a faint diffraction 
pattern due to y-Fe2O; overlying the stronger 
a-Fe pattern. 

The samples were heated in cylindrical electric 
furnaces equipped with automatic controls which 
maintained the temperature constant to within 
+1°C. The ends of the furnaces were closed with 
loose fitting refractory plugs which permitted free 
access of air. Ordinary air without purification 
was used for oxidation. Some of the samples were 
placed in the furnace in small alundum boats; 
others were supported by the thermocouple wires 
used to measure the temperature of the sample. 

The oxidation temperature was determined by 
means of a thermocouple the elements of which 


*H. R. Nelson, J. Chem. Phys. 5, 252 (1937). 

* Abraded surfaces were used partly because they most 
nearly represent a type of mechanically worked surface 
often met in practice. More important, however, is the fact 
that freshly abraded iron surfaces give an easily repro- 
ducible and interpretable pattern in which changes can 
readily be followed. Polished surfaces, on the other hand, 
8ive a poor pattern which, irrespective of its interpretation, 
's certainly not characteristic of metallic iron. With etched 
surfaces difficulty was experienced in removing the last 
traces of etching products which, it is reasonable to 
Suppose, might profoundly affect the oxidation process. 
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Fic. 1. Graphical representation of the results of 
heating iron specimens for various times at various 
temperatures. 


were wedged or peened into two small holes 
drilled in the sides of the specimen. Air tempera- 
tures near the test surface were sometimes 
measured by asecond thermocouple or a mercury- 
in-glass thermometer placed close to but not 
touching the specimen. The temperature close to 
the sample as measured by this method was not 
appreciably different from that of the sample 
itself as determined by the attached couple. 
The diffraction camera is similar to that 
described by Thomson and Fraser! with a hot 
cathode source of electrons. Constant accelerating 
potentials of 30 to 40 kv are provided by a 
transformer, rectifier and filter system. The 
grazing incidence or “‘reflection’’ method was 
used with a sample-to-plate distance of 40 cm. 


RESULTS 


The results on abraded iron surfaces are 
summarized graphically in Fig. 1. Each point 
represents a single experiment in which a 
specimen was heated continuously at the indi- 
cated temperature for the time given on the 


10G, P. Thomson and C. G. Fraser, Proc. Roy. Soc. 
A128, 641 (1930). 
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Fic. 2. Observed and calculated diffraction patterns. 
Heights of the curves represent estimated relative 
intensities. 


logarithmic time scale after which the surface 
was examined in the diffraction camera. The 
character of the resulting diffraction pattern is 
indicated by the type of plotted point. A circle 
means that rings due to iron were present in the 
pattern. Crosses indicate the presence of the 
primary oxide y-Fe:O;, while squares represent 
a-Fe,O3. The lines AB and CD are drawn simply 
to separate the time-at-temperature area into 
fields within which the indicated type of diffrac- 
tion pattern was observed. 

A freshly abraded iron sample gives a diffrac- 
tion pattern consisting of broad diffuse rings and 
a few fairly sharp rings. As has been shown 
previously® this pattern is due to iron and the 
primary oxide. The latter is present in exceedingly 
small crystals as evidenced by the lack Of 
résolution of the diffraction rings (cf. Fig. 2 E, F, 
and G). The iron rings are fairly sharp unless they 
happen to lie so close to oxide maxima as not to 
be resolved therefrom. 

When such an abraded surface is heated in air 
at temperatures of 120°C or lower there is very 
little change in the pattern after heating periods 
as long as 1700 hours. (See Fig. 2, H and J.) In 
general the oxide pattern becomes slightly 
stronger with a slight tendency for the diffuse 


maxima to be resolved into their component 
rings. The resolution is never good even with the 
longest oxidation periods used and there was 
little difference in the patterns obtained from 
specimens heated at 80°C and 120°C for 1700 hr. 

The presence or absence of iron was determined 
by the presence or absence of the strong (110), 
(211) and (321) reflections of iron in the pattern. 
The (110) ring was never completely resolved 
from the (400) ring of the oxide but the other 
two could be checked easily since they do not 
fall near oxide maxima. As can be seen in plots H 
and J of Fig. 2 iron was present after heating 
1700 hr. at 120°C or below. 

When an abraded iron specimen is heated at 
temperatures above about 130°C the observed 
patterns follow the sequence shown in plots £, D, 
C, B of Fig. 2. Rings due to iron disappear fairly 
rapidly in heating times which depend upon the 
temperature. The disappearance of the iron is 
accompanied by a strengthening of the primary 
oxide pattern but there is no appreciable 
sharpening of the diffuse rings before they begin 
to be replaced by new rings due to rhombohedral 
a-Fe,O;. After still longer heating periods the 
primary oxide is entirely replaced by a-Fe.0;, 
Figs. 3 and 2B. This sequence of patterns has 
been cbserved at temperatures as low as 133°C 
and as high as 210°C. At the higher temperatures 
these changes take place before the specimen has 
reached temperature equilibrium, consequently 
the time required for the various changes could 
not be determined accurately. These results show 
conclusively that abraded iron surfaces develop a 
film of a-FesO; on heating in air at temperatures 
above about 130°C. The slope of the right-hand 
end of the line AB suggests that this oxide can be 
produced at even lower temperatures if suf- 
ficiently long oxidation periods are used. One 
specimen heated for 850 hr. at 127°C gave a well- 
developed a-Fe,O; pattern. It is not known why 
this result falls so far out of line. 

At higher temperatures the surface film of 
a-Fe,O; undergoes a change characterized by a 
gradual replacement of the diffuse diffraction 
maxima by new, sharply defined rings due to the 
same oxide (see Figs. 4 and 5). There is also a 
change in the type of preferred orientation which 
is discussed later. The transition from the diffuse 
to the sharp a-Fe.O; pattern is coincident with 
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OXIDATION OF IRON 


Fic. 3. a-Fe2O; (diffuse pattern) from_iron heated in air at 
133°C for 800 hr. 


the development of the first order blue inter- 
ference color. Sharp diffraction rings begin to 
appear on the diffuse pattern when the surface 
reaches the purple stage. The diffuse maxima are 
entirely gone when the surface has become deep 
blue. 

Above and to the right of CD in Fig. 1 only the 
sharp a-Fe.O3 pattern is obtained ; below and to 
the left is a region in which the pattern due to 
this oxide is at least partly diffuse. It will be 
apparent from Fig. 1 that the data do not 
determine the shape of the curve CD. A straight 
line could be drawn to separate the fields nearly 
as well as the curve. The fact that at 193° this 
curve has not deviated sharply to the right 
indicates that the sharp a-Fe.O; pattern charac- 
teristic of higher temperatures might be obtained 
from samples oxidized sufficiently long at temper- 
atures well below 200°. Experimental data in the 
range 170° to 180° and 1000 to 10,000 hr. are 
necessary to test this point.” 

A diffraction pattern due to some unknown 
film was obtained in a number of experiments. 
(See Figs. 2/ and 6.) The pattern is characterized 
by three strong diffuse rings which are rather 
sharply defined on the inside but poorly defined 
on the outside. The innermost ring which corre- 
sponds to a spacing of about 3.5A, is much 
stronger than the others. This pattern was always 
obtained when samples were heated at tempera- 
tures of 125°C or below in one of the several 
furnaces used for these tests. It could not be 
reproduced in any of five other furnaces which 


Fic. 4. a-Fe,O; (diffuse pattern) from iron heated in air at 
240°C for 4 hr. 


Fic. 5. a-FeO; (sharp pattern) from iron heated in air at 
240°C for 14 hr. 


were used. Above 130°C this particular furnace 
gave results entirely consistent with those ob- 
tained in the other furnaces. The possibility that 
some unknown material from the walls of the 
furnace condensed on the surface and produced 
this pattern was tested by heating the test sample 
simultaneously with copper and nickel specimens 
placed beside it in the furnace. In every case the 
iron gave the diffuse ring pattern while the 
control surfaces did not. Moreover, the unknown 
pattern was obtained on a wide variety of iron 
surfaces, e.g. on electrodeposited iron films, on 
mirrors formed by vacuum evaporation and 
condensation of iron, and on etched and polished 
as well as abraded surfaces. 

Attempts to analyze and explain this pattern 
have failed and at present it must be considered 
an anomalous result. It has not been included in 
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Fic. 6. Unknown diffuse ring pattern from iron heated in 
air at 122°C for 1850 hr. 


the experimental results plotted in Fig. 1. The 
effect has been mentioned here, however, because 
it suggests the possibility that an amorphous 
oxide may be formed when iron is heated at 
temperatures below 125°C under certain special 
conditions which this experiment has failed to 
define. 

A few experiments were performed on highly 
polished surfaces and various etched surfaces 
were tried. Although the results were similar to 
those obtained on abraded specimens there was 
evidence that the initial films left by such 
methods of preparation were too variable to give 
reproducible results (see footnote’). 

In a previous paper® it was shown that the 
primary oxide crystals stand in a definite 
orientation relationship with the iron crystals on 
which they form. It follows, therefore, that on 
abraded iron surfaces in which the iron crystals 
are randomly oriented the primary oxide should 
also be randomly oriented. Such in fact was 
observed on all abraded surfaces irrespective of 
the temperature at which the specimens were 
heated. By similar reasoning we conclude that 
a-Fe,O; crystals should be randomly oriented if 
they bear a definite orientation relationship to 
the primary oxide or to the underlying iron. 
Actually, however, the a-Fe,O; patterns always 
show a low degree of preferred orientation. In the 
initial a-FesO; (open squares in Fig. 1) the [111] 
directions of the rhombohedral crystals tend to 
be perpendicular to the surface. Referred to 
hexagonal axes this means that basal planes tend 
to lie parallel to the surface of the specimen. 

The appearance of the sharp a-Fe2O; pattern is 


accompanied by the development of a new type 
of preferred orientation characterized by a slight 
tendency for (211) rhombohedral planes to lie 
parallel to the surface. As can be seen in Fig. 5 
this is a very low degree of orientation. 

In order to determine the possible effect of 
orientation in the original iron surface the oxi- 
dation of electrodeposited surfaces was studied. 
The crystals formed when iron is plated onto 
polished or abraded iron tend to stand with [111] 
directions approximately perpendicular to the 
surface. The primary oxide film forms with an 
equally high degree of orientation because of the 
definite spacial relationship between it and the 
underlying iron. The preferred orientation ex- 
hibited by a-Fe2O; films grown on such surfaces, 
however, was found to be exactly the same as 
was observed on abraded surfaces. 

No careful attempts were made to determine 
the interference color of the specimens as a 
function of time and temperature. As has already 
been mentioned, however, the curve CD in Fig. 1 
very nearly coincides with the appearance of the 
first order blue. The lowest temperature at which 
blue was obtained was 193°C. 1000 hr. at 138° 
gave a faint yellow similar to that observed after 
about 4 hr. at 180°. Specimens heated at 120° 
became perceptibly darker as the oxidation was 
prolonged. At 80° there was only a barely 
detectable darkening after 1700 hr. 


DISCUSSION 


During the course of this work Iimori" reported 
electron diffraction experiments in which he 
studied oxide films formed on etched iron surfaces 
heated in air at various pressures and tempera- 
tures. Below 340°C he found only the cubic oxide 
which he believes is y-Fe:O3. At first sight this 
seems difficult to reconcile with the present 
results. However, Iimori reports only four experi- 
ments performed at atmospheric pressure and of 
these only two fall in the temperature range 
covered in Fig. 1. He finds y-Fe2O3 after 1 hr. at 
250°C or 24 hr. at 150°C. Both of these points 
fall within the region of diffuse a-Fe,O3 patterns 
obtained on abraded surfaces. The discrepancy is 
undoubtedly due to the different surface treat- 


1 T. Iimori, Sci. Papers, Inst. Phys. Chem. Research, 
Tokyo 34, 60 (1937). 
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ment used in the two investigations and one is 
forced to the important conclusion that results on 
one type of surface cannot be applied to another. 
In fact Iimori states that on polished surfaces he 
obtained sharp a-Fe,O; patterns as low as 250°C. 

The real purpose of this investigation, viz., to 
study the nature of a critical temperature near 
200°C, has not been achieved. The results give no 
evidence of a critical temperature in this region. 
Long oxidation periods in the neighborhood of 
180° might establish a lower temperature limit 
for the formation of the sharp a-Fe2O; diffraction 
pattern. On the other hand comparatively short 
oxidation periods were sufficient for Vernon to 
show the existence of a critical temperature. It is 
apparent, therefore, that the difference between 
the oxidation mechanism above and below 200°C 
does not manifest itself in structural differences 
in the topmost layer of the oxide film. 

The preferred orientations observed in the 
a-FesO; films cannot be due to pseudomorphism 
on underlying crystals as is the case with the 
primary films. Since the axis of orientation is 
always normal to the surface with complete 
randomness in other directions it seems likely 
that the crystalline orientation is in some way 
associated with the diffusion of oxygen through 
the oxide film or, more probably, the diffusion of 
iron outward to form new layers of oxide at the 
outer surface of the film. Pfeil'® has shown that 
iron diffuses outwards through thick layers of 
oxide when ferrous materials are heated at high 
temperatures. This suggests a possible reason for 


© L. B. Pfeil, J. Iron and Steel Inst. 119, 501 (1929). 
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the appearance of the sharp a-Fe,O; pattern with 
its new orientation at the higher temperatures. 
The sharp diffraction rings mean that the new 
oxide crystals are much larger than those which 
gave the diffuse pattern. Careful examination of 
plates showing the transition from the diffuse to 
the sharp pattern leads to the conclusion that 
this transition is not due to crystal growth or toa 
process resembling recrystallization. Instead it 
appears that new crystals with the new orienta- 
tion grow on top of the old film. If this interpre- 
tation is correct, iron must diffuse through the 
oxide film to the surface. The curve CD in Fig. 1 
may, therefore, be some sort of measure of the 
rate of diffusion of iron through the oxide layer 
and the curve’s lower temperature limit, if there 
is such a limit, would determine the temperature 
at which the diffusion rate becomes negligible. 
The presence of diffraction rings due to metallic 
iron after oxidation periods of 1700 hr. at 120°C 
and less is surprising. In the ‘‘reflection’’ method 
the electrons pass through jagged ridges and 
projections left by the abrasion process. In 
passing through such a projection an electron 
would pass first through an oxide film, then 
through the unoxidized core and finally emerge 
through the oxide film. In such a path the 
electron might be scattered in either the iron or 
the oxide. Since, however, electrons cannot pene- 
trate more than several hundreds of angstroms of 
solid material without suffering multiple scat- 
tering or absorption, the presence of iron rings in 
the diffraction pattern fixes an upper. limit of 
about 200A for the thickness of the oxide film. 
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The effects of temperature and of heat treatment on 
Pyrex chemically resistant glass have been studied by 
measuring the rate of diffusion of helium and of hydrogen 
through a tube of about one millimeter wall thickness. 
Measurements were made with helium between 180°C and 
590°C, and with hydrogen at 512°C. The influence of gas 
pressure was determined and the diffusion rate was found 
to be directly proportional to the pressure for both helium 
and hydrogen. 

The diffusion data throw light on the processes of 
molecular rearrangement involved in stabilization of the 
glass. Below 548°C the rates of diffusion (R) were increased 
about 10 percent by giving the glass sufficient heat treat- 
ment to assure ‘‘stabilization,’’ i.e., equilibrium in the 
processes of dissociation and association. This increase 
occurred despite the fact that such heat treatment normally 
increases the density of the glass and therefore points to a 
chemical theory rather than a simple theory of holes 
governing the diffusion process. 

The data also showed a continuous linear dependence 
for the stabilized glass of log Ron 1/T from 440°C-590°C. 


This range extends into the region of a true viscous liquid 
and it follows that the properties of a stabilized glass are 
continuous with that of a viscous liquid and that there is 
no need to postulate a separate glassy state. The “‘trans- 
formation point’’ mentioned by some investigators is 
merely the temperature of rapid strain release of an unan- 
nealed sample. 

At or near 440°C the log R versus 1/T curve takes a 
lesser slope (for lower temperatures) and this change is 
attributed to the loss of rotational or vibrational freedom 
in the silicate complex, analogous to that frequently found 
in dielectric studies on solids. Activation energies have 
been calculated for the two temperature ranges. 

At 512°C the rate of helium diffusion is approximately 
45 times that of hydrogen, and it is believed that chemical 
forces delay the progress of the hydrogen molecule through 
the silicate network. 

The general conclusion is drawn that gas diffusion fur- 
nishes a highly sensitive method for studying the molecular 
condition of silicate glasses and similar materials. 


HE principal feature distinguishing the 
glassy state from the ordinary liquid state 

is a very high viscosity which delays the estab- 
lishment of equilibrium of molecular processes in 
glass. The result of this delay is that the proper- 
ties of a glass depend to an important degree on 
its thermal history or heat treatment. Examples 
are expansion coefficient, viscosity,” elastic prop- 
erties such as Young’s modulus,’ electrical con- 
ductivity, absorption spectrum,® and dielectric 
constant.* When the temperatures are sufficiently 
high equilibrium is very rapidly attained in a 


1M. E. Spaght and G. S. Parks, Glass VIII, ‘‘The 
Coefficient of Thermal Expansion of Boron Trioxide,”’ J. 
Phys. Chem. 38, 103-10 (1934). 

2H. R. Lillie, ‘‘Viscosity-Time-Temperature Relations in 
Glass at Annealing Temperatures,’’ J. Am. Cer. Soc. 16, 
619-631 (1933). 

3N. W. Taylor, E. P. McNamara and J. Sherman, ‘“‘A 
Study of the Elastico-Viscous Properties of a Soda-Lime- 
Silica Glass at Temperatures near the ‘Transformation 
point,’’’ J. Soc. Glass Tech. 21, 61-81 (1937). 

4J. T. Littleton and W. L. Wetmore, ‘‘The Electrical 
Conductivity of Glass in the Annealing Zone as a Function 
rl and Temperature,’’ J. Am. Cer. Soc. 19, 243-245 
(1936). 

5 W. Weyl and E. Thiimen, ‘‘The Constitution of Glass, 
II. Theoretical Basis of Glass Coloration,’’ Glastech. Ber. 
11, 113-120 (1933). 

6 W. Weyl, personal communication. 
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glass and it cannot be distinguished in its be- 
havior from an ordinary viscous liquid. 

The nature of the liquid state partakes in part 
the properties of crystals,” and in part the prop- 
erties of a condensed gas, but its study has proven 
rather difficult. One may say that at present we 
have less knowledge of the structure of liquids 
than we have of either gases or crystalline solids. 
One of the difficulties lies in the fact that dy- 
namic equilibrium is rapidly attained in ordinary 
liquids and that it is impossible to “freeze in” 
or preserve a certain molecular condition, for 
later study. In glasses this is not so. One may by 
quenching or rapid chilling preserve a high tem- 
perature structure down to low temperatures, or 
one may anneal a glass at any desired tempera- 
ture and study the properties of the molecular 
configurations so developed. These characteristics 
of glass are utilized industrially in many ways. 

There are various physical methods for in- 
vestigation of the structure of glass. Several have 


7G. W. Stewart, ‘X-Ray Diffraction in Liquids,” Rev. 
Modern Physics 2, 116-122 (1930). ‘‘ Molecular Structure as 
Interpreted by X-ray Diffraction Measurements in Liq- 
uids,’’ Chem. Reviews 6, 483-501 (1929). 
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DIFFUSION OF GASES THROUGH GLASS 


been mentioned. X-ray analysis* has demon- 
strated a certain regularity of atomic or ionic 
arrangement. This technique, however, presents 
a statistical picture of average interatomic dis- 
tances and does not show the result of heat treat- 
ments in much detail, i.e., it is not particularly 
sensitive to heat treatment of the glass. 

A method which depends on individual atomic 
properties, or on the properties of small groups 
of atoms or ions, rather than on the whole com- 
plex, is the method of gas diffusion through the 
solid. The progress of the diffusing atom is de- 
layed by its binding (electrical interaction) with 
neighboring host atoms and the rate of diffusion 
is proportional to a quantity e~4/*7 where A is 
the activation energy necessary to overcome 
these atomic forces. Diffusion thus becomes a 
property subject to study and mathematical 
treatment by the usual procedures of chemical 
reaction kinetics, and precise data can be ob- 
tained which lead to rather direct interpretation 
in terms of atoms and molecules. 

As an example of the success of the diffusion 
method in the study of solids we may mention 
the work of W. R. Ham,’ on the diffusion of 
hydrogen through iron and nickel. It was found 
that “the magnetic changes at about 360° (Ni) 
and 750° (Fe) were clearly defined by diffusion 
data. In addition, iron exhibits changes in the 
isobars at 945°, 900° and at between 200° and 
550°C. A definite difference in temperature is 
observed between Ar; and Acs, and this appears 
to be a true hysteresis. The other points are all 
ascribed to electronic changes of the iron atom 
and it is suggested that the well known maximum 
in the tensile value of iron at 200-350° is ex- 
plained in part by one of these electronic shifts. 
The general conclusion is reached that the isobars 
from hydrogen diffusion tests furnish an accurate 
means of tracing phase or electronic changes in 
iron and nickel, and that the change in slope 
of the isotherms provides a means in some in- 
stances of determining small amounts of foreign 
elements in solution in the metals.” 


*B. E. Warren and J. Biscoe, “Fourier Analysis of X-ra 
Patterns of Soda-Silica Glass,” J. Am. Cer. Soc. 20, 250-65 


(1938), 
*W. R. Ham, “The Diffusion of Hydrogen Through 


(sa and Iron,’’ Trans. Am. Soc. Metals 25, 536-64 


CHEMICAL THEORY OF DIFFUSION 


The simple physical picture of diffusion as a 
movement of atoms through open channels or 
cracks which are large compared to atomic di- 
ameters, does not offer an explanation for much 
of the most accurate and reliable experimental 
work in this field. On the other hand, several 
workers, notably Dushman and Langmuir,'® and 
later Eyring," have treated the diffusion process 
by the same methods as for chemical reaction 
rates. The fundamental assumption is made that 
the diffusing atom moves from lattice point to 
lattice point through the solid and is subjected to 
the influence of strong electrical fields of a 
periodic nature. In order for a diffusing atom to 
break away from its immediate environment and 
move to a new position of stability in the solid, 
it requires an amount of energy, £, to activate it. 
By application of the Boltzmann principle, the 
probability that a diffusing atom can cross this 
energy barrier is proportional to e~#/*7, where k 
is the Boltzmann constant. For a gram molecular 
weight of diffusing atoms, since R= Nk and 
A=NE where N is the Avogadro number, the 
expression becomes e~4/”7, 

The quantity A is sometimes called the activa- 
tion energy, the energy of transfer, or the heat of 
loosening per gram atomic weight. This is where 
the chemistry of diffusion may enter the theory, 
for obviously those atoms which can exert a 
strong attraction for the diffusing atom (large A 
value) will delay the progress of the diffusing 
atom through the solid. Thus for example at 
500°C, helium has a rate of diffusion through 
silica glass about 22 times faster than has hydro- 
gen, as shown by Williams and Ferguson” “‘Were 
it a case of effusion the ratio would be less than 
unity,” that, is approximately \/2/4=0.7, since 
the molecular weights of helium and hydrogen 
are 4 and 2 respectively. 

This large difference is to be expected on the 
theory that H, forms temporary attachments to 


10S, Dushman and I. Langmuir, ‘“‘The Diffusion Coeffi- 
cient in Solids and its Temperature Coefficient, ’’Phys. 
Rev. 20, 113 (1922). 

1H, Eyring, ‘Viscosity, Plasticity and Diffusion as 
Examples of Absolute Reaction Rates,’’ J. Chem. Phys. 4, 
283-291 (1937). 

12 G. A. Williams and J. B. Ferguson, ‘‘The Diffusion of 
Hydrogen and Helium through Silica Glass and other 
Glasses,"’ J. Am. Chem. Soc. 44, 2160-2167 (1922). 
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614 W. FTAYLOR 
the oxygen ions which are of greater stability 
than that which helium may form with these ions. 
The two gases have nearly the same diffusion 
velocities through rubber (CioHis) which has a 
low oxygen content." In fact, because of the 
chemical inertness of the helium atom, it is 
probable that the only forces which tend to bind 
it to the ions of the glass are the so-called van 
der Waals forces which arise from the interaction 
of dipoles between the two kinds of atoms. 

Our own experiments show that at 512°C the 
rate of helium diffusion through Pyrex glass is 
approximately 45 times as fast as that of hydro- 
gen, and this despite the fact that the diameters 
of He and of He are nearly the same, namely 4.00 
and 3.94A, respectively. These values are com- 
puted from the densities of the two liquids, at 
—271.5° for He and at — 258° for He. 

The theory that weak bonds between diffusing 
atoms and host atoms lead to relatively rapid 
diffusion has been also emphasized in another 
paper by one of the authors," in discussing the 
diffusion of gold, silver bismuth, thallium, and 
tin into lead, and the self-diffusion of lead. 
Further evidence that diffusion is essentially 
controlled by chemical factors has been given by 
W. R. Ham! who has shown that a mixture of Ne 
and He diffuses through metallic iron very much 
more rapidly than either gas alone. This effect 
may be explained by assuming that Nz and He 
together form a complex molecule, possibly of an 
amide type, which has less attraction for the 
iron atoms than have the two H atoms which 
form when N¢ is not present. 

In connection with the problem of the glassy 
state there are several questions which might be 
clarified by gas-diffusion studies. The first ques- 
tion was concerned with the nature of the 
so-called ‘‘transformation point: is it a true 
thermodynamic property of a glass or is it simply 
interpreted as such because of a confusion of rate 
effects with equilibrium properties? Most expan- 
sion measurements are dynamic, in the sense that 
the temperature is rising or falling at a finite rate. 
It was considered desirable to measure some prop- 


13 J. D. Edwards and S. F. Pickering, ‘‘Permeability of 
Rubber to Gases,’’ U. S. Bureau of Standards, Scientific 
Paper 387 (1920). 

4N. W. Taylor, ‘‘Reactions Between Solids in the 
OTe of a Liquid Phase,”’ J. Am. Cer. Soc. 17, 155-163 
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erty such as diffusion under constant temperature 
conditions, above and below the ‘‘transformation 
point.’’ The second question is related to the first 
one and may be stated: ‘“‘What influence does 
thorough annealing have on the permeability 
of glass to gas, and how shall the results be 
interpreted?” 


MATERIALS AND METHODS 


Earlier work on gas diffusion through glass has 
shown that hydrogen and helium diffuse more 
easily than do other gases. For this reason, and 
because of the relative simplicity of atomic 
structure of these gases, they have been chosen 
for this study. For the glass, Pyrex chemically- 
resistant glass was chosen, primarily because it 
was easily made into a suitable shape for the 
diffusion sample and also because it could readily 
be sealed to the rest of the vacuum system. From 
the standpoint of theoretical interpretation a 
glass of simpler chemical composition, such as 
fused silica or boric oxide would have been 
better, but from the practical standpoint Pyrex 
is of greater interest. 


EXPERIMENTAL METHOD 


The diffusion sample was a 41.0 cc Pyrex test- 
tube, having 94.6 cm? effective diffusing area and 
1.1 mm wall thickness, sealed into a larger Pyrex 
tube as in Fig. 1. It was placed in a wire wound 
cylindrical muffle furnace, whose temperature 
was maintained uniform to +0.25°C by three 
thermocouples, at ends and center, and by auto- 
matic control to +0.5°C. The temperature of the 
“sample” was read from a thermocouple placed 
inside the inner tube. The high pressure side con- 
nected with a manometer and gas supply, the 
low pressure side with a calibrated McLeod gauge 
and mercury diffusion pump and oil pump. 
Vacuum stopcocks (Hg seal) were used through- 
out. The apparatus was thoroughly tested for 
leaks, using air, before introduction of helium or 
hydrogen. These gases were from commercial 
sources and were purified before use, by contact 
with charcoal at liquid-air temperatures. This 
treatment probably removed water vapor 
hydrocarbons, but would have no effect on any 
nitrogen which might be present. The speed of 
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Fic. 1. Diffusion sample. Direction of diffusion indicated 
by arrows. 


nitrogen diffusion through glass is extremely 
slow in comparison to that of helium, however. 

The procedure consisted of trapping the dif- 
fused gas for definite intervals from 15 seconds 
to 30 minutes depending on the temperature, 
compressing this gas in the McLeod gauge and 
reading its volume at barometric pressure and 
room temperature which were of course known. 
At least 3 hours were devoted to each tempera- 
ture, so that the glass had an opportunity to 
reach a steady state during that period. The 
diffusion data therefore represent essentially 
static conditions as compared with ordinary 
expansion studies. 


PRESENTATION OF DATA 


In view of the fact that the properties of a 
glass depend in part on the thermal history of 
the sample, the schedule of heat treatment will 
be given. 

The rate measurements fall into five groups: 


I. August 4~August 11, 1936. Temperature rising from 
406°C to 590°C. 

II. August 11-August 27. Temperatures falling from 
590°C to 180°C. 

III. August 28-September 11. Temperatures rising and 
falling between 352° and 570°C for 23 cycles. This 
series will be described as “‘the isotherms.”’ 

IV. September 12-September 14. Temperatures rising 
from 187° to 367°C. 

V. September 14-September 27. Glass held at 350° to 
400°C. September 27-October 1. Temperatures 
rising from 242° to 544°C. 


Group L August 4-August 11. Temperature 
rising from 406° to 590°C (curve A) 


After the whole vacuum system had been 
demonstrated to be free of leaks, measurements 
of diffusion rate for helium were begun at a 
temperature of 406°C. Since the glass sample had 
not been made strain-free by annealing at a 
temperature in the annealing range (approxi- 
mately 545°C or above), the glass was really in a 
condition of strain. For this reason the rates of 
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diffusion are less than those obtained later. The 
experimental points are plotted in Fig. 2, Curve 
A. It will be seen that at temperatures of 548°C 
or above the rates obtained on the strained glass 
are almost the same as those obtained later in 
the Group II series (Curve B) which presumably 
represents strainfree glass. Below 548°, however, 
i.e., between 406° and 548°C, the helium diffusion 


rates are about 10 percent less than those of 


Group II. Previous workers in the field of gas 
diffusion through glass have nearly always failed 
to remove the condition of strain in the glass 
before taking precise measurements, and it is 
therefore doubtful whether much of the pub- 
lished work on gas diffusion should be accepted 
in a quantitative sense. 


Group II. August 11—-August 27. Temperature 
falling from 590° to 180°C (curve B) 


The data of Group II were taken over a period 
of 17 days, measurements made at 590, 584, 
574, 567°C, etc., until a low point of 180°C had 
been reached. As in Group I and all other groups, 
each temperature was held constant for approxi- 
mately 3 hours before the rate data were ob- 
tained. At all temperatures down to about 525° 
ample time was available to permit complete or 
nearly complete stabilization (molecular re- 
arrangement) of the glass to a condition charac- 
teristic of that temperature. At or near 525°C 
there was a very slight change in direction of the 
log R versus 1/T curve which thereafter continued 
as a straight line down to about 440°C. Below 
440°C the log R versus 1/T line was again straight, 
but had a somewhat lesser slope. Note that on 
Fig. 2 it diverges from the broken line which is a 
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continuation of the upper slope. At about 225°C 
the line again changed direction, but as the rate 
measurements are less certain in this low tem- 
perature region, no emphasis will be laid on this 
effect, at present. The significance of the 525°C 
and 440°C breaks will be discussed later. Data 
below 337°C are not shown in Fig. 2. 


Group III. The isotherms, August 28-Sep- 
tember 11 


The relation between diffusion rate and gas 
pressure was a subject of special interest, not 
only to permit reducing all rate data to a common 
pressure basis of 1 atmosphere, but also in con- 
nection with the more general problem as to the 
“molecular” state of the helium atom as it 
worked its way through the glass. For these 
reasons a series of isotherms were made at 5 
different pressures, 16, 32, 48, 64 and 80 cm Hg. 
respectively, and at 4 different temperatures, 
355°, 460°, 525° and 570°C. It was convenient to 
work at one helium pressure over the four tem- 
peratures, then to reduce its pressure and work 
back over the same four temperatures in reverse 
direction, etc. This resulted in heating the glass 
from 352 to 568, cooling to 355, heating to 570, 
cooling to 305 and finally heating to 571°C. The 
rate values are plotted in Fig. 3, and in Fig. 2, 
Curve C. 

These cycles of alternate heating and cooling 
apparently produced no major effect on the 
diffusion rates except a slight increase, ca. 1 per- 
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cent above 450°C and about 23 percent at or 
near 355°C. The “transformation point” at about 
525 was wiped out, however, and its disappear. 
ance was confirmed by later tests, Group V. 

The accuracy and reproducibility of the diffu. 
sion rates of Group III are shown by Fig. 2 where 
the plotted points fall nicely on the straight lines 
at each temperature. 

Table I exhibits values of dm/dt for various 
values of p, and also shows values of R, the 
characteristic rate per atmosphere pressure gradi- 
ent (g helium per cm? per sec. per mm wall 
thickness), calculated from the data of the two 
columns by means of the relation R= 76/p dm, di. 
The temperature is also shown. It will be seen 
that the values of R are practically constant, 
except for small shifts due to temperature. 

In the range 440°-600°C, the data of Curve C, 
Fig. 2, may be satisfactorily expressed by the 
equation 

—log R=9.075+1562/T, 


where R has the units already stated, and T is 
the absolute temperature, degrees Kelvin. This 
equation leads to rate values approximately |! 
percent higher than those given in Table I, but 
is representative of the later data. 

The linear relation between rate of diffusion 
and gas pressure is not surprising, since helium 


TABLE I. Rate of helium diffusion through stabilized 
Pyrex glass. 


PRESSURE 


Temp. (°C) (cm He) dm/dt R 
355 80.11 3.05 X 10-8 2.82 
355.5 63.71 2.42107” 2.89 
357 47.96 1.83 x 10-2 2.94 x10" 
355 31.69 1.23 X 107? 2.88 x 107" 
359 15.96 0.59 x 10-2 2.90 x 
465 80.97 7.62 X10-” 6.29 x 
455.5 63.87 5.18 5.96 x 10 
461 48.32 3.86 X 6.15 X10" 
459 31.59 2.56107” 6.1210 
460.5 16.03 1.28 6.10 X10 
521.5 80.48 9.65 X 10-2 8.92 x10 
526 64.08 7.69X 10-2 9.15 x10 
525 47.28 5.73 9.21 x10™ 
526 32.01 3.86 X 10-2 9.20 
529 16.16 1.93 107” 9.27 
568 80.76 12.50 10-2 11.6210 
573 64.15 9.84 10-2 11.84 X10 
569 47.65 7.39 X 107? 11.70x10™ 
570 32.18 4.93 10-” 11.66 
571 15.86 2.46 X 10-2 11.83 x10 
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is normally monatomic. It may be mentioned 
that the solubility of helium in silica glass is in 
direct proportion to the pressure.’ 


Group IV. Temperatures rising from 187 to 367°C 


Following the series of isotherms, a group of 5 
determinations were made at temperatures be- 
tween 187 and 367°C. They are consistent with 
Curve C, Fig. 2. This confirms the small increase 
in rate due to the heat treatment received by the 
glass sample in the 23 weeks subsequent to 
August 27. 


Group V. Temperature rising from 242 to 544°C 
following 2 weeks heat treatment at 400- 
350°C 

After the series known as Group IV had been 
completed the sample stood one week at 400°C 
and one week at 350°C. Rate measurements were 
then made at approximately 244°, 440°, and 
542°C. They indicate a further small increase, 

about 1-1.5 percent, over the rates in Group III 

and IV taken two weeks previously. They sup- 

port the conclusion that further stabilization 
took place in the glass during this period. 

Secondly, they confirm the disappearance of the 

“transformation point.’’ On the other hand, they 

appear to confirm the reality of a change in slope 

of the log R versus 1/T curve at or near 440°C. 


DISCUSSION 


Interpretation of the energies of activation 


The heats of activation computed from the 
log R versus 1/T curves for the temperature 
ranges 350°-440°C and 440°-600°C are 6480 and 
7150 calories per mole respectively. For the un- 
annealed glass it is 6920 calories (Temp. 400- 
500°C). These values may be compared with Van 
Voorhis"'® value of 8700 cal. for helium through 
Pyrex (200° to 500°C), and W. D. Urry’s!? value 
of 5340 cal. for helium through unannealed Pyrex 
in the temperature range 81°-283°C. 


*G. A, Williams and J. B. Ferguson, ‘‘The Solubility of 
Helium and Hydrogen in Heated Silica Glass and the 
Relation of this to the Permeability,’’ J. Am. Chem. Soc. 
46, 635-639 (1924), 

'* C. C. Van Voorhis, ‘The Diffusion of Helium Through 
(mae Widely Different Glasses," Phys. Rev. A23, 557 


rh. W. D. Urry, “Further Studies in the Rare Gases. I. 
fn Permeability of Various Glasses to Helium,” J. Am. 
hem. Soc. 54, 3887-3901 (1932). 


617 


Two possible factors may determine the mag- 
nitude of the activation energy barriers. One is 
the attraction of helium for the ions of the glass, 
which we believe to be small compared to the 
second which is the repulsive force acting on the 
helium atom as it moves from “‘hole’’ to “‘hole’”’ 
in the molecular network. This quantity might 
be quite large in magnitude if the passages 
through which the helium has to move are very 
small. The term “‘hole’’ as used here indicates a 
space of about atomic dimensions which is not 
occupied by any of the ions constituting the 
glass. Just as a silicate crystal is formed of SiO, 
tetrahedra linked, by sharing of corners, into 
chains and rings and sheets, these larger units 
being joined through positive ions, so we believe 
a silicate glass to have a similar structure, the 
essential difference being that there is no repeat- 
ing “unit cell’ in the liquid, but only a random 
arrangement of SiO, tetrahedra in space, and 
secondly, that the centers of gravity about which 
the atoms oscillate or rotate, themselves move 
around in a kind of Brownian motion. This sort 
of structure provides small ‘‘open spaces’ or 
“holes” of atomic size in which the diffusing 
helium atoms may easily stay. On the other 
hand, considerable energy may be required by a 
helium atom to allow it to squeeze past its nearest 
oxygen neighbors and find a similar ‘‘hole.’’ The 
magnitude of the repulsive energy developed as 
the helium atom passes through a small molecular 
orifice will depend on an inverse 6th or 7th power 


TABLE II. 


DIAMETER OF SQUARE ENERGY (CAL/MOLE) 


3.6A 13,920 
3.8 8,393 
4.0 4,250 
4.2 2,183 
4.4 1,078 


law integrated over all the ions which act on a 
given helium atom as it moves from the center of 
its “hole” to the tightest point in the orifice. At 
present exact calculation of this energy from 
fundamental electrostatic or wave-mechanical 
theory cannot be made. Analogous calculations 
have, however, been made by Barrer,'* who has 

18 R. M. Barrer, ‘‘The Mechanism of Activated Diffusion 


use Silica Glasses,’ J. Chem. Soc. (London) 378-86 
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considered how the size of a hole between four 
helium atoms arranged at the corners of a square 
affects the energy required to bring in a fifth 
helium atom from an infinite distance normal to 
the square so that it finally rests in the center of 
the plane. Table II is taken from Barrer’s article, 
page 385. 

It will be seen that not only are the repulsive 
energies several thousand calories in magnitude 
and therefore comparable to those actually ob- 
served in the present work, but also they are 
quite sensitive to small changes in the size of the 
hole when the size is nearly the same as that of 
the atom which is being squeezed. This may 
partly explain why heat treatment of the glass 
changes the activation energy for helium diffusion. 


Changes resulting from heat treatment 


1. The “transformation point.’’—One question 
to which the diffusion data have given a clear 
answer is the nature of the “transformation 
point.’’ The log R versus 1/T plot for the stabi- 
lized glass (Curve C) is exactly linear from about 
440°C to 600°C and shows no discontinuity or 
change in direction at or near 550°C which is the 
approximate temperature at which the expansion 
coefficient shows large increase. Since each diffu- 
sion measurement was of several hours duration 
the data may be regarded as static, in contrast to 
the dynamic expansion measurements. The con- 
clusions must be drawn that a viscous liquid and 
a stabilized glass are continuous with one another 
in properties and that the ‘‘transformation point”’ 
is purely a property of the unstabilized glass. It 
is merely a temperature at which the viscosity is 
low enough to permit rapid stabilization by 
molecular rearrangement, thus leading to equi- 
librium between the forces of association and 
dissociation. These views are in harmony with 
those of Lillie,? Littleton,? and Weyl'® but con- 
trary to those of Berger®® and others. 

2. Strained and annealed glass.—Various work- 
ers have noted that heat treatment affects the 


19 W. Weyl, ‘“‘On The Constitution of Glasses,’’ Glastech. 
Berichte 10, 541 (1932). 

20 A. Klemm and E. Berger, ‘‘The Change in Trans- 
mission of Glasses with Temperature Above and Below the 
Transformation Point,’’ Glastech. Berichte 14, 194-206 
(1936); E. Berger, ‘‘Contributions to the Theory of Glass 
Formation and the Glassy State,’’ J. Am. Cer. Soc. 15, 
647-677 (1932). 


W. TAYLOR AND W. RAST 


permeability of glass to gases. Tsai and Hogness”! 
state ‘‘We find that silica glass when heated toa 
high temperature undergoes a permanent change 
which results in a noticeably greater permeability 
for helium. The permeability depends somewhat 
on the previous heat treatment of the glass.” 
Other workers have recorded that various samples 
of the same glass, e.g. fused silica, gave different 
permeabilities, without ascribing any reason for 
it. An examination of Fig. 2 shows that the rate 
of diffusion of helium through the Pyrex glass 
after annealing (Group II) was about 10 percent 
greater than that through the unannealed Pyrex 
(Group I). This large increase occurs despite the 
fact that Pyrex glass becomes denser (of higher 
specific gravity), as a result of annealing. This is 
shown by the density values 2.236 (annealed) 
and 2.233 (chilled).” 

These densities do not show the maximum 
spread which may be obtained under special 
conditions. On the other hand, they are in line 
with the usual experience that the chilled glass 
is less dense, and they fail to explain the lower 
permeability of this glass. The conclusion must 
follow that the chilled or unannealed glass is in a 
dissociated and unsaturated state which for some 
reason, either purely geometric or because of 
stronger chemical binding forces, delays the 
progress of the diffusing atoms. The idea of 
chemical unsaturation of the unannealed glass is 
in agreement with the views of W. Weyl based 
on absorption spectra. It has also been known 
for many years that unannealed glass is more 
susceptible to chemical attack then the annealed 
product. 

The 440°C effect.—Figure 2, Curve C, shows 
two straight lines which meet at approximately 
440°C. Above this temperature the rate of diffu- 
sion is expressed by 


—log R=9.075+1562/T, 
while below it, at least down to 350°C, the rate 


is given by 
—log R=9.285+1415/T. 


This change in slope occurs whether one ap- 
proaches 440°C from above or below and with 


2L. S. Tsai and T. R. Hogness, ‘Diffusion of (Gases 
“oe Fused Quartz,” J. Phys. Chem. 36, 2595-2600 
1 
22H. R. Lillie, personal communication. 
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essentially static measurements (three hours per 
point). 

Since there is no known pecularity in the 
properties of helium at 440°C it is believed that 
the effect is due to the Pyrex glass. The following 
hypothesis is tentatively advanced. The thermal 
energy of the glass at any temperature consists in 
atomic and molecular motion of translation, os- 
cillation and rotation. It is possible that at or 
near 440°C a new degree of freedom is gained, 
i.e., above this temperature some molecule like 
SiO, may rotate about one axis with respect to 
the rest of the complex, while below 440°C it 
may not. 

A recent paper by Smyth** presents evidence 
pertaining to rotation in crystal lattices. In 
NH,NO; and NaNO; the NO; ion begins to 
rotate in the plane of its ions at 84°C and at 275° 
respectively. Similarly in NH,Cl, NH,Br and 
NH,I rotation of the NH, ion sets in at —30°, 
—38° and —42°, respectively. Many other ex- 
amples are discussed by Smyth, who points out 
that the shape of the molecule is an important 
factor, and that symmetrical groups like NHs,, 
PH, (or SiO,) rotate more easily than groups of 
lower symmetry. 

No claim is made for great sharpness (small 
temperature range) of the ‘‘break’’ at 440°. As 
far as the experimental evidence is concerned it 
may actually extend over several degrees. In 
crystals the effect may be sudden or gradual. The 
writers have no data on dielectric constants for 
Pyrex in this temperature range, but such infor- 
mation would be useful to test this hypothesis. 
Since we regard glass as a highly viscous liquid, 
the following statement by Smyth is of interest. 
“From a molecular point of view, the essential 
difference between a solid with rotating molecules 
and a liquid is that the molecules of a solid are 
vibrating about points fixed in a lattice, while 
molecules of the liquid are vibrating about points 


* C. P. Smyth, ‘The Dielectric Constants of Solids and 
Molecular Rotation,” Chem. Review 19, 329-361 (1936). 
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which migrate slowly and irregularly through the 
liquid.” 

It is obvious that more glasses and in particular 
glasses of simple composition should be studied 
in order to get more light on this hypothesis. It 
seems highly probable that the increased loosen- 
ing of the lattice because of a rotation (or an 
extra degree of vibrational freedom) would per- 
mit more rapid diffusion of helium, as observed. 

It is also clear that exact specific heat or 
similar information on crystalline silicates on 
simple chemical composition would be of value in 
connection with the above mentioned hypothesis. 


DIFFUSION OF HYDROGEN 


Two determinations were made at 787.5°K 
(512.5°C) on the rate of diffusion of hydrogen. 
The first, with the hydrogen under a pressure of 
78.04 cm Hg, led to a value dm/dt=1.908 x 10-*, 
from which R=dm/dt=1.858 X 10-" when calcu- 
lated for a pressure drop of 1 atmosphere 
(76 cm Hg). The second, with a measured ai 
hydrogen pressure of 41.44 cm, gave a value of 
dm/dt=1.017X10-" from which R=dm/dt= 
1.865 x 10-" when calculated for a pressure drop 
of 1 atmosphere. 

The agreement between these two results 
shows that the rate of diffusion of hydrogen is 
directly proportional to the pressure, and not to 
the square root of the pressure as is the case when 
hydrogen diffuses through metals. It appears 
that the two hydrogen atoms diffuse through 
Pyrex glass in pairs, rather than as separate 
atoms. 

A comparison of the rate of diffusion of hydro- 
gen and helium at 787°K (512°C) indicates that 
the helium diffuses about 45 or 46 times as fast 
as does hydrogen at this temperature. 

Even at 120°C the diffusion of hydrogen into 
glass is perceptible. W. Weyl** has produced 
fluorescence by reduction of dissolved AgsO to 
the free atoms by hydrogen at this temperature. 


2W. Weyl, “Contribution to the Fluorescence of 
Glasses,”’ Sprechsaal 70 (1937), No. 46. 
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The free volume is defined as the total integral over that 


part of the potential energy of the molecule in the liquid 
which is due to thermal displacements of the center of 
gravity of the molecule from its equilibrium position. The 
free angle is the corresponding integral over angular dis- 
placements of the molecule. The free volume, V; is related 
to the velocity of sound, u, in the liquid by 


u (liquid) =u (gas) (V/ Vy). 


This equation is used to derive a formula connecting the 
sound velocity in the liquid with its thermal conductivity. 
The quantity, RT/p exp (AH/RT), gives the product of the 


free angle ratio and the free volume, rather than the free 
volume itself. Here p is the vapor pressure and AH is the 
heat of vaporization. The free volumes from sound veloc- 
ities agree with those determined by independent methods, 
Specific heats, entropies of vaporization, the spectroscopic 


(Received July 21, 1938) 


sound velocity, free volumes and the function RT p 
Xexp (AH/RT) are examined from the point of view of 
restricted rotation of the molecules in liquids. Lack of 
free rotation of molecules suffices to explain the ab- 
normalities of the liquids examined. The dielectric con- 
stants of polar liquids are interpreted assuming restricted 
rotation and the following formula is derived: 


=1—((1—cos /2)*, 
where yw; is the apparent dipole moment of the molecule 
in the liquid, yu, is the dipole moment as determined from 
measurements in the gas phase, and 4, is the polar angle 
related to the free angle ratio, 52, by 
62=(1—cos 0;)/2. 
The concept of restricted rotation of the molecules in the 


liquid accounts satisfactorily for the observed dielectric 
polarizations of water and methyl alcohol. 


INTRODUCTION 


N the present paper certain kinetic theory 
formulae for gases are modified in such a way 
as to make them applicable to liquids. The point 
of view employed is to treat the molecule moving 
in its free volume in the liquid as equivalent to 
the molecule moving in the total volume in the 
gas phase. Using this notion it is possible to de- 
rive a new formula for the velocity of sound in 
liquids in terms of the velocity in the gas phase, 
the total volume and the free volume of the 
liquid. An equation very similar to that given by 
Bridgman? for the thermal conductivity of liquids 
also follows from such considerations. 

The term, free volume, as employed here de- 
notes the total integral over that part of the 
potential energy which is due to thermal dis- 
placements of the center of gravity of the 
molecule from its equilibrium position. The 
name, free angle, is introduced to refer to the 
corresponding integral over angular displacement 
of the molecule around its center of gravity. For 
linear molecules the free angle ratio is the free 


1 Present address, University of Rochester, Rochester, 
New York. 


2P. W. Bridgman, The Physics of High Pressures 


(Macmillan Co., 1931). 


angle divided by 47, and for nonlinear ones an 
analogous definition holds. The free angle ratio 
is a convenient quantity to employ since it is 
given by the ratio of the partition functions for 
the restricted rotation of the molecule in the 
liquid, and its free rotation in the gas phase. The 
free angle ratio thus extends from a lower value 
(greater than zero) to unity depending on 
whether the molecule is held tightly in its pos'- 
tion or is allowed to turn freely about its axes. 
For linear molecules the free angle ratio includes 
two degrees of rotational freedom, and is desig- 
nated by the symbol, 52. For more complicatec 
molecules, three degrees of rotational freedom 
must be taken into account and the correspond- 
ing symbol is 63. 

To a first approximation both the free volume 
and the free angle ratio may be thought of 
regions in which completely free displacement‘ 
occurs, the regions being surrounded by infinite 
potential walls. This is the hard spheres a> 
proximation previously employed by Eyring 
Hirschfelder and Stevenson’: * to yield the 


3H. Eyring and J. O. Hirschfelder, J. Phys. Chem. #! 
249 (1937). 

4J. O. Hirschfelder, D. P. Stevenson and H. Eyriss 
J. Chem. Phys. 5, 896 (1937). 
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pression for the free volume, 


1 2RT 
v= ’ (1) 


where V is the molal volume of the liquid, R is 
the gas constant per mole, 7 is the absolute 
temperature, and FE is the thermodynamic 
energy. It is, however, not necessary to consider 
the molecules as hard spheres, and Kincaid and 
Eyring® have given a treatment of liquid mercury 
in which the size of the molecules is presumed to 
vary both with temperature and volume. The 
dynamic potential curve’ in this case is a parabola 
with a flat bottom. 

For monatomic liquids, the free volume is 


defined by, 
V;=RT/p exp (AH/RT), (2) 


where p is the vapor pressure of the liquid, and 
AH is the heat of vaporization. For liquids com- 
posed of polyatomic, nonspherical molecules, the 
expression on the right of Eq. (2) is equal to the 
product of the free volume, and the free angle 
ratio. When the molecules in the liquid do not 
rotate freely, this quantity will necessarily be 
smaller than the free volume. 


THE Use or KINETIC THEORY EQUATIONS 
FOR LIQUIDS 


Velocity of sound in liquids 


The velocity of sound, u, in any homogeneous 
medium is given by the general hydrodynamic 
formula: 


u=(v/B.)}, (3) 


where v is the specific volume and 8, is the 
adiabatic compressibility defined by 8,=—1/v 
X(dv/0P),. For the special case of an ideal gas, 
(3) becomes 


u=(RTy/M)}, 


where y is the ratio of the specific heat at con- 
stant pressure to that at constant volume, and 
M is the molecular weight of the gas. 


’John F. Kincaid and H. Eyring, J. Chem. Phys. 5, 
587 (1937). 

* The dynamic potential curve is defined as the curve for 
which the corresponding partition function, when raised 
to the 3Nth power, is the translational partition function 
for the liquid as a whole aside from an exponential potential 
energy factor. 
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Equation (4) may be compared to that for @, 
the average kinetic theory velocity of the 
molecules. 


(5) 


It is seen that u is proportional to é but is slightly 
smaller since the factor (8/7)! is always greater 
than y. This is what one might suppose since a 
wave propagated by matter would hardly be ex- 
pected to travel faster than the molecules which 
carry it. However, when the velocity of sound in 
liquids is compared to the kinetic theory velocity 
of the molecules it is found that u for most 
liquids is greater than ¢ by factors ranging from 
five to ten. For example for benzene at 25°C, 
is 2.83 X 104, u(gas) is 1.88 X 104, while w(liquid) 
is 13.010‘ cm/sec. Fig. I illustrates the mecha- 
nism responsible for the fact that the velocity of 
sound in liquids can be so much greater than the 
kinetic theory velocity of the molecules. The 
wave front is assumed to travel from the edge 
of molecule A to the adjacent edge of molecule B 
with the velocity given by Eq. (4). As A collides 
with B, however, the signal is transmitted almost 
instantaneously to the opposite edge of molecule 
B. Since the ratio of the total distance to the free 
space between two points in a liquid is given by 


nr’ 


Fic. 1. Illustration of the mechanism which explains 
the observation that the velocity of sound in liquids is 
greater than the kinetic theory velocity of the molecules. 


the ratio (V/V,;)', this leads at once to the 


equation,’ 


7 Equation (1) of this paper is Eq. (13) of the paper by 
Eyring and Hirschfelder and was there derived from the 
expression Vs!=2(V!—d) assuming (dd/V)7 equals zero, 
i.e., the hard sphere model. Eq. (1) is obtained by the 
same procedure from the expression V;=(V!—d) provided 
(dd/8V)r has the appropriate negative value. This corre- 
sponds to getting more free volume from increasing ( V!—d) 
> ae one would get if the sides of the free volume were 
flat faces. Since the sides of the free volume are spheres 
this sign of (dd/dV)r is obvious a priori. The necessity of 
going beyond the simplest hard sphere model has been 
encountered in two previous papers. See Eq. (22) of 
Kincaid and Eyring‘ and Eq. (21) of Hirschfelder, Steven- 
son and Eyring.’ However, it is possible that Eq. (6) should 
be employed in the more general form, u(liquid) =cu(gas) 
(V/V;)4 when c isa proportionality constant not necessarily 
equal to unity. 
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u(liquid) = u(gas)(V/V,;)} 
=(RTy/M)(V/V;)*. (6) 


Equation (6) may be employed either to find 
the velocity of sound in liquids if V; is known 
from some other source, or to determine free 
volumes if sound velocity measurements have 
been made. A detailed comparison of free vol- 
umes computed by this and various other meth- 
ods will be made in a later section. 


Thermal conductivity of liquids 


A straightforward application of the kinetic 
theory of gases gives the relation for the thermal 
conductivity K, 


K(gas) =3(V/V;)'2Lc,, (7) 


where (V/ V;)} is the number of molecules per cc, 
L is the mean free path, and c, is the specific heat 
per molecule. It is found, however, that Eq. (7) 
gives values which are too low for the thermal 
conductivity of gases, and Loeb® gives the cor- 
rection factor }(9y’—5), where y’ is the ratio of 
Cp to cy. Eq. (7) may now be rewritten, 


K (gas) =43[4(97'—5) (8) 
a relation which may be tested in the form, 
K(gas) (9) 


since the viscosity, 7, is given by 3(N/V)méL. 
Eq. (9) has been observed to give excellent 
agreement with experiment.® 

In order to convert Eq. (8) into a form ap- 
plicable to liquids the following identifications 
may be made. (N/V) becomes the number of 
molecules per cc of liquid rather than per cc of 
gas. The average velocity ¢ must be multiplied 
by the ratio, (V/V,)', for the same reasons as are 
given above for justifying Eq. (6). The distance 
that the energy is carried is now (V/N)! rather 
than L, the mean free path for the gas. c, for the 
gas is replaced by acy, where a is the accommoda- 
tion coefficient which takes account of the num- 
ber of degrees of freedom which come into equi- 
librium in the thermal conduction process. 
Making all these substitutions we have: 


K (liquid) =1/12(9y’ —5)(N/V)(8RT/xM)} 
X(V/V)V/N) ac»), (10) 


8 Leonard B. Loeb, Kinetic Theory of Gases (McGraw- 
Hill Book Co., 1927). 
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TABLE I. Comparison of the observed thermal conductivities 
of liquids at 30°C with those computed from Eq. (12). 


K (oss.) 

SUBSTANCE K(oss.)9 K(CALC.)!0 K(CALC.) 
Methyl Alcohol 21.1 X 10% c.g.s. Units 22.3 0.95 
Ethyl Alcohol 18. 19.2 94 
Propyl Alcohol 15.4 17.3 89 
Butyl Alcohol 16.7 13.3 1.26 
Iso-Amyl] Alcohol 14.8 14.1 1.05 
Ether 13.7 10.9 1.26 
Acetone 17.9 16.3 1.10 
Carbon Bisulfide 15.9 16.1 0.99 
Ethyl Bromide 12.0 12.8 94 
Ethyl Iodide 11.1 10.6 1.05 
Water 60.1 51.7 1.16 


a formula which reduces to 


K (liquid) = 1/12(9y’—5)(8/7y)! 
X(N V)iu(acy). (11) 


Here u is the velocity of sound in the liquid, y’ is 
the effective value of the ratio of c, to c, for 
thermal conduction in the gas, and y is the same 
quantity for the transmission of sound. 

In order to have formula (11) reduce to the 
one successfully employed by Brdigman,? ac, 
must be assigned the value 3k per molecule. For 
the polyatomic liquids which Bridgman studied, 
this is understandable if only the kinetic energy 
of translation and rotation transfer appreciable 
energy in the process of thermal conduction. 
Using a value of c, equal to 3k, y’ becomes 4/3 
and (11) reduces to, 


K (liquid) =0.931(1/7)!3k(N/V)3u. (12) 


This formula has the same form as that employed 
by Bridgman,? but the numerical coefficient is 
about 12 percent smaller, as a result of the factor 
0.931/y! replacing unity in his expression. 
Table I shows values for K at 30° for a number 
of liquids computed by (12) for comparison with 
Bridgman’s experimental values. The agreement 
is somewhat better than Bridgman’s earlier com- 
putations. Since carbon bisulfide is a linear 
molecule, the computation for this substance 
was made with a value of ac, equal to 2.5k, and 
a corresponding value of y’ of 1.40. The excellent 
agreement of the observed value with the one 
calculated in this way suggests that the liquid 
rare gases would be interesting substances with 
which to check the theory. A value for ac, of 


9P. W. Bridgman, Proc. Am. Acad. 59, 109 (1923). 

10 y values secured from K. Bennewitz, and W. Rossner, 
Zeits. f. physik. Chemie B39, 126 (1938). Sound velocity 
measurements from Bridgman® or Freyer, Hubbard and 
Andrews.” 
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about 3k would also be obtained from the trans- 
lational contribution alone, if potential energy 
were included. Thus from data for the condensed 
rare gases, one will be able to determine whether 
or not potential as well as kinetic energy comes 
into equilibrium in the thermal conduction 
process in liquids. 

Equation (12) not only works remarkably well 
for the variation of K from liquid to liquid at 
atmospheric pressure, but, as Bridgman has 
pointed out, it also gives approximately the 
temperature variation of K. It is therefore quite 
surprising that it does not predict the pressure 
effect correctly. Whereas most liquids increase 
their thermal conductivity by approximately a 
factor of 2 on going from atmospheric pressure to 
12,000 atmospheres,” the formula predicts an in- 
crease by about a factor of 4. This is plausibly 
explained by a decrease in the value of a due to 
quantization of the mass motions of the molecules 
at high pressures. The effective heat capacity, 
ac,, would have to become 1.5% in order to ex- 
plain the observations at the highest pressures 
to which the experiments extend. 

Debye" has treated the problem of the thermal 
conductivity of solids, with special reference to 
effects at low temperatures. The essential differ- 
ence between his treatment and the one given 
above is that the identification of L with (V/N)} 
in the derivation of Eq. (12) is equivalent to 
assuming that each molecule acts independently 
while the essence of Debye’s treatment lies in 
his assumption that the molecules are strongly 
coupled in such a way as to produce elastic waves 
of long wave-length. The effect of coupling, 
enabling the heat impulses to be transmitted by 
distances greater than the’intermolecular dis- 
tances without scattering is, however, important 
only at low temperatures. For most liquids it 
appears to be more satisfactory to treat the 
problem in essentially the manner indicated. 


A comparison of values for free volumes of 
liquids computed by various methods 


Since more complete data exist for benzene 
than for any other nonmetallic and nonpolar 
liquid, its free volume was computed by all avail- 


"P. Debye, Wolfskehl-Vortrige zu Géttingen (Leipzig, 
1914), Anhang p. 46. See also, Max Born, Atomtheorie des 
Festen Zustandes (B. G. Teubner, Leipzig, 1923), p. 707. 
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TABLE II. The free volume of benzene as a function of 
temperature at atmospheric pressure. 


Vs FROM! 13 63V¢ FROM" Vs FRoM!5 
Ea. (6) Eg. (2) Ea. (1) 
(SouND (VAPOR (ENERGY VOLUME 
7°C VELOocITyY) PRESSURE) COEFFICIENT) 
2 0.193 cc 0.214 cc 
.272 0.238 cc .288 
50° 419 423 
127 .633 
100° 1.24 .99 


able methods. The agreement between the vari- 
ous methods is seen by inspection of Table II to 
be very satisfactory. C, for liquid benzene exceeds 
the gaseous value by 5 calories per mole per de- 
gree. This suggests lack of free rotation about 
two axes for the molecules in the liquid. On the 
other hand the sound velocity method gives a 
free volume only slightly greater than the product 
of free volume with free angle ratio as determined 
from vapor pressures, indicating a free angle ratio 
not greatly less than unity. 

In general reliable values for heats of vaporiza- 
tion are not available over a sufficiently large 
temperature range to make possible a significant 
comparison of the absolute magnitude and tem- 
perature coefficients of the free volumes from 
vapor pressure with free volumes from sound 
velocity. However, although Eq. (1) is an ap- 
proximate relation based on a very simple model, 
inspection of Table II reveals good agreement 
between values based on this relation and those 
calculated from vapor pressures. On the strength 
of this agreement, free volumes were computed 
by Eqs. (1) and (6) for a number of liquids for 
which the necessary data are given by Freyer, 
Hubbard and Andrews.” The results are shown 
in Table III. Although the computation of free 
volumes from vapor pressures and heats of 
vaporization is not a practical procedure for 
many liquids because of uncertainties in experi- 
mental heats of vaporization, the converse calcu- 

12 Sound velocity data from E. B. Freyer, J. C. Hubbard 
and D. H. Andrews, J. Am. Chem. Soc. 51, 759 (1929). 

13 Values of y for computing yw for the gas were secured 
from the frequencies given by R. C. Lord and D. H. An- 
drews, J. Phys. Chem. 41, 141 (1937). 

144E. F. Frock, D. C. Ginnings and W. B. Holton, Nat. 
Bur. Stand. J. Research 6, 881 (1931). 53 is the free angle 
as Data from 25° to 65°C from R. E. Gibson and John F. 
Kincaid, J. Am. Chem. Soc. 60, 511 (1938). Data at other 
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TABLE III. Values of the free volume for various liquids. 


Sus- 
STANCE ACETONE ETHER 


CHLOROFORM 


CARBON CARBON 
TOLUENE TETRACHLORIDE DISULFIDE 


Vs FROM 
Ea. (6) 


Vs FROM 


Vs FROM 
Eg. (1) 


Ea. (6) 


Vs FROM 
Eg. (1) 


Vs FROM 
Ea. (1) 


Vs FROM 


0.54 
63 


0.45 


lation of heats of vaporization from vapor pres- 
sures and “‘sound velocity” free volumes appears 
to be a useful method. 

In order to test the applicability of the _" 
velocity method of computing free volumes for a 
metallic liquid, free volumes of mercury as com- 
puted from Eq. (6) were compared with those 
given by the partition function previously pub- 
lished.> The expression for Vy is, 


Here 0 is the characteristic temperature of the 
liquid and Vp is identified as the volume of the 
liquid at the melting point. Taking © equal to 
79°, a value consistent with the experimental 
entropy of fusion, the free volumes shown in 
Table IV are obtained. It may be noted that the 
agreement between the values at both low and 
high pressures is quite satisfactory. This is of 
particular interest since the free volumes for 
mercury are of a different order of magnitude 
than those of the nonmetallic liquids discussed 
above. However, only a more exact investigation 
of mercury can establish how satisfactory these 
values really are. 

In order to test the theory for another mon- 
atomic liquid, an attempt was made to utilize 
the somewhat inadequate data for liquid argon. 
The entropy of liquid argon at its melting point 
at atmospheric pressure, as computed from the 
data of Euken"® is 12.5 E.U. The general formula 
for the entropy of a liquid is,'” 


S=4R+R In (14) 
16 A. Euken, Verh. D. phys. Ges. 18, 4 (1916). 


1 This formula in the exact form given depends on a 
contribution of R to the entropy of fusion resulting from 


and this leads to a value of 0.15 cc for the free 
volume of liquid argon at its melting point. 
Taking the specific heats at constant volume and 
constant pressure given by Euken and Hauck,'* 
and the volumes determined by Baly and 
Donnan,!® we estimate 8, at the melting point to 
be 1.29X10-"° c.g.s. unit. This corresponds to a 
value of u equal to 0.74105 cm/sec. and a free 
volume of 0.35 cc. The agreement to within a 
factor of two is considered satisfactory in view of 
the uncertainties in the data. The free volume 
for argon at its boiling point was determined 
from its heat of vaporization to be 1.25 cc. Since 
other methods yield much smaller values, it ap- 
pears to us that the accepted figure for the heat 
of vaporization must be too low. A value of 1.65 


TABLE IV. Free volumes for liquid mercury at different 
temperatures and pressures. 


P(aTM.) Vs FROM Ea. (6) Vy FROM Ea. (13) 
— 39° 1 0.0081 cc 
0 1 0.0120 cc .0120 

100° 1 .0220 .0256 
200° 1 .0434 
300° 1 .0718 
357° 1 .0904 

0° 1 .0120 .0120 
0° 1000 .0117 .0116 
0° 2000 0114 .0114 
0° 3000 0111 O111 
0° 4000 .0107 .0109 
0° 5000 .0104 .0107 
0° 6000 .0101 .0105 
0° 7000 .0098 .0103 


the communal sharing of the free volume. It also is de- 
pendent on the assumption that each molecule moves on 
a parabolic dynamic energy curve, and on the vanishing 
of the derivative (0E,/dT)7, where E, is the static energy 
previously defined (Ref. 5). It is believed that these sim- 
plifications do not lead to serious error in this case. 

18 A, Euken, and F. Hauck, Zeits. f. physik. Chemie 
134, 161 (1928). 

WE.C.C. Baly, and F. G. Donnan, J. Chem. Soc. 81, 
907 (1902). 
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kcal./mole, 10 percent greater than that given 
in the International Critical Tables, is indicated 
by the theory. 


Determination of free volumes from spectroscopy 


The frequency of vibration of molecules in 
liquids is related to the free volume by the 
approximate relationship, 


V3 
(2amkT)! 


h 


(15) 


which connects the velocity of sound in the liquid 
with the energy hy for a translational degree of 
freedom. Thus one would expect the light 
scattered from liquids to show the familiar Stokes 
and anti-Stokes lines of the Raman spectrum 
corresponding to energies from 10 to 100 cm. 
Such lines have been observed”® and as Debye”! 
has remarked they show the same correlation 
between frequency and sound velocity to be 
expected from Eq. (15). This method of arriving 
at free volumes will no doubt become of in- 
creasing importance. 


TaB_e V. Difference between the specific heats at constant volume for the ideal gas and the liquid for various substances at 
different temperatures.*® 
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RESTRICTED ROTATION OF MOLECULES 
IN LIQUIDS 


Heat capacity as a criterion of restricted rotation 


Most of the available data for the quantity 
C,(liquid) minus C,(ideal gas) has been sum- 
marized in Table V. It is seen that in all cases 
the heat capacity at constant volume for the 
liquid is greater than the ideal gas value and 
that this inequality persists up to the critical 
point. Considering first the two examples of 
monatomic liquids which are given in the table, 
we note differences of 2.6 and 3.1 calories per 
mole per degree for argon and mercury respec- 
tively at their melting points. These are just the 
differences to be expected if three degrees of 
translational freedom for the molecule in the 
vapor have become vibrations in the liquid. 

The data for liquids composed of more complex 
molecules can be interpreted in terms of hindered 
rotation of the molecules. If rotational degrees of 
freedom of the molecules in the vapor become 
librations” in the liquid we would expect a 
specific heat difference between liquid and ideal 
gas of about one calorie per mole per degree for 
each degree of freedom so affected. This may, of 
course, be modified by anharmonic terms in the 


Mercury® 
Cy(L1Q) Cy(L1Q) 

—C,(IDEAL —C,(IDEAL 

T°K GAS) 


ARGON*3 


T°K 


CARBON DISFULIDE!® 


CARBON TETRACHLORIDE!® 26 
C,(LIQ) C,(LIQ) 
—C,(IDEAL —C,(IDEAL 


84° (2.6) 3.1 160° 4 Meiting Point 
370° 2.6 210° (3.9) 280° 3.8 

- 500° 24 265° 3.4 310° 3 

87° 2.5 630° 2.3 319° 2.3 349° (3.0) Boiling Point 
156° 1.1 Critical Point 
T CHLOROFORM!2+ 26 BENZENE! 27 ETHYLENE: 25, 28 

Cy(L1Q)  Cy(L1Q) Cy(L1Q) 

— Cy(IDEAL —C,(IDEAL —Cy(IDEAL 
T°’K GAS) T°K GAS) GAS) 
203° (6.7) 278° 5.3 104° 5.3 Melting Point 
250° (6.3) 300° 5.0 125° 5.3 
290° 5.3 325° 4.9 150° 5.3 
334° 4.7 353° 4.9 170° 5.3 Boiling Point 

283° (4.2) Critical Point 


x E. Gross, and M. Vuks, Nature 135, 998 (1935). 
"?. Debye, Chem. Rev. 19, 171 (1936). 
* Libration is used here to mean restricted rotation 


* A. Euken and Hauck, Zeits. f. physik. Chemie 134, 161 (1938). 


* Paul C. Cross, J. Chem. Phys. 3, 825 (1935). 


* G. B. Kistiakowsky and Falih Nazmi, J. Chem. Phys. 6, 18 (1938). 


* Robert D. Vold, J. Am. Chem. Soc. 59, 1192 (1935). 


4 R. C. Lord and D. H. Andrews, J. Phys. Chem. 41, 149 (1937). 
a Mathies E. Haas and Gebhart Stegeman, J. Phys. Chem. 36, 2125 (1932). 
The data for specific heats at constant volume for the vapors are, in most cases, calculations from spectroscopic data. 
The references give either the computations or the necessary data for making them. Extrapolated values are given in 
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potential functions. On this basis the difference 
of 5.3 calories per mole per degree shown by 
benzene and ethylene near their melting points 
is interpreted as being due to vibrations for five 
of the degrees of freedom and free rotation about 
the axis perpendicular to the plane of the mole- 
cule for the other. It is significant that carbon 
tetrachloride, a symmetrical molecule, has a 
smaller difference between the specific heats of 
liquid and vapor than chloroform, a relatively 
unsymmetrical one. The data indicate that even 
carbon disulfide does not rotate freely in the 
liquid, but the figures at the two lowest tempera- 
tures recorded must be interpreted with some 
caution since they represent long extrapolations. 
The same is true of the chloroform data. The 
high value of 4.2 calories per mole per degree for 
the difference between the observed and the ideal 
gas value for ethylene at the critical point may 
be due to a change with temperature of associ- 
ation into double and triple molecules, assuming, 
of course, that the data are correct. Any hypoth- 
esis of hindered rotation at such large volumes 
appears unlikely. 

Since the data indicate that the effect of rising 
temperature, and thus of increasing volume, is to 
reduce the potential restricting rotation, it is of 
interest to examine the effect of volume changes 
alone on this potential. This can be done by 
examining the effect of pressure on the specific 
heat at constant volume of those liquids for 


TaBLE VI. The heat capacity at constant volume of liquid 
benzene under pressure.*® 


°K 1 BAR | 250 BARS | 500 BARS | 750 BARS | 1000 BARS 
298.1° | 11.0 11.5 11.7 
308.1° | 10.9 11.0 11.2 11.5 $1.7 
318.1° | 10.8 11.1 11.3 11.6 11.8 
328.1° | 10.8 11.2 11.4 11.7 11.9 
338.1° | 10.7 1 11.4 11.8 12.1 


30 The units are calories per mole per degree. The changes 
of specific heat with pressure were computed from the 
data of Gibson and Kincaid® using the thermodynamic 
equations, 


dC,/dP=—Td?V/dT? and 


where C, is the specific heat at constant pressure, and a 
and £ are defined as, 
a=1/V(dV/dT),, B=—1/V(AV/dP)r. 

The values for C, at one bar are the weighted mean of the 
data given in Landolt-Bérnstein, Tabellen. The contribu- 
tions to the specific heat made by the internal vibrations 
were computed from the frequencies given by Lord and 
Andrews.?? 
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which adequate data are available. Table \] 
shows values of C, for liquid benzene at various 
pressures and temperatures. The contribution of 
the internal vibrations to the specific heat have 
been subtracted from the thermodynamic (,, 
It is necessary to assume that this contribution 
is independent of pressure since no data are 
available on the effect of high pressures on the 
frequencies corresponding to the normal modes 
of vibration of molecules. Assuming that this 
effect will not change the qualitative aspects of 
the data we have the result that the heat 
capacity at constant volume is considerably 
greater at 1000 bars*' than at low pressures. The 
fact that C, rises to as high a value as 12.1 
calories per mole per degree appears to indicate 
that at high pressures there is restricted rotation 
even about the axis perpendicular to the plane 
of the ring. 


The free angle ratio from spectroscopy 


When the velocity of sound in liquid water is 
used to compute V; at room temperature 2 
value of about 0.4 cc is secured. On the other 
hand, the product 63;V; as determined from 
heats of vaporization and vapor pressures at the 
same temperature is about 0.02 cc. A reasonable 
interpretation is to assign the discrepancy of 1 
factor of about 20 to a small free angle ratio. A 


TABLE VII. Free volumes for liquid water. 


Vy FROM? 63V FROM 
Ea. (6) Eg. (2) 1/63°3 
(SOUND (VAPOR (FROM SPEC- 
FC VELOcITY) PRESSURE) TROSCOPY) 
0 0.44 cc 0.0095 cc 26 0.25 cc 
10 43 .012 25 36 
20 42 .014 25 ao 
30 42 .017 24 Al 
40 43 .021 24 50 
50 .025 24 .60 
75 .040 24 .96 
100 .062 24 


3t The bar, or metric atmosphere, is equal to 10° dyne 
per square centimeter. This is approximately one norm 
atmosphere. 

% Sound velocity data from J. C. Hubbard and A. | 
Loomis, Phil. Mag. 5, 1177 (1928). 

33 1/5; computed from the relation 


x ( 
1—exp (—/c500/kT) 


with a value of ABC, the products of the principal ™~ 
ments of inertia, equal to 6.1 x 107°, 
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MOLECULES 


check of this hypothesis is possible. Cartwright* 
has observed frequencies of 167 and 500 cm in 
liquid water at room temperature. These fre- 
quencies were interpreted by him as arising from 
restricted rotation of the molecules. Since three 
degrees of rotational freedom must be accounted 
for, at least one of the rather broad absorption 
regions probably arises from two frequencies, and 
Stearn and Eyring** have employed a doubly 


degenerate frequency at 167 cm™ to represent - 


approximately this absorption band. Using this 
approximation we have the values for the free 
angle ratio shown in Table VII. These were 
secured by taking the ratio of the partition func- 
tion for free rotation of the water molecule and 
the product of the Einstein function correspond- 
ing to frequencies of (2) 167 cm and 500 cm. 

It may be seen that the free volume from the 
vapor pressure, the heat of vaporization, and the 
observed frequencies is almost exactly equal to 
the free volume from sound velocity at room 
temperature, the point at which the spectro- 
scopic observations were made. The discrepancies 
at other temperatures indicate that the fre- 
quencies corresponding to restricted rotation 
should be greater than those observed at lower 
temperatures, and smaller at higher tempera- 
tures. Such a temperature drift will probably be 
difficult to observe spectroscopically because 
rising temperature simply makes the absorption 
bands broader and more diffuse, as has been ob- 
served.*® 


The dielectric constant as a criterion for re- 
stricted rotation 


In order to have available a new tool for the 
determination and interpretation of restricted 
rotation of molecules in liquids, formulas con- 
necting the free angle ratio with the dielectric 


f f exp | —— Cos cos O10 sin 
0 Yo ART 


IN LIQUIDS 


7. 
5° 100° 150° 780° 
Angle 


Fic. 2. wey, plotted against the polar angle 6,. This 
ratio is determined from equation (17) and the relation, 
m=? F/3kT which holds for the gas. 
constant are given. The importance of the free 
angle ratio in determining vapor pressures and 
boiling points of liquids may be seen by inspec- 
tion of Fig. 2 where 62 is plotted as a function of 
the polar angle, @;:, through which the molecule 
is free to turn. Here 4; is analogous to the angle 
of colatitude as measured on a globe from the 
north pole down to a given “‘parallel,’’ and é¢ is 
the corresponding solid angle divided by 47. 
The two are connected by the equation, 


§2=1—cos 0,/2. (15) 


Here 62 is analogous to 6: previously defined for 
linear molecules, but differs from 6; for nonlinear 
ones since it does not take account of the re- 
stricted rotation of the molecule about the axis 
in which the dipole lies. If the potential curves 
for the rotation of the molecule about the three 
axes are symmetrical, then 52 is equal to the two- 
thirds power of 53. In some actual cases this rela- 
tion is approximately correct as we shall see later. 
We next consider dipoles free to orient under 
the influence of the field in the region between 
the pole and colatitude @;, assuming the poles 
are distributed at random with respect to the 
field direction. The expression for the average 
moment m, using classical statistics is then: 


Sin ¢o. (16) 


1 27 
f 
Jo 


“C. H. Cartwright, Phys. Rev. 49, 470 (1936). 
3) Stearn and Henry Eyring, J. Chem. Phys. 5, 


* Paul C, Cross, John Burnham, and Philip A. Leighton, 
59, 1134 (1937). 


22 uF 
f f exp cos sin 
0 Yo kT 


Here 612 is the angle between the dipole and the 
applied field, F, and yz is the dipole moment. 4, 
and yg; are the polar coordinates of the dipole 
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with respect to an axis chosen in such a way that 
the dipole can move with a potential energy 
which depends only on @. This potential energy 
arises from interaction with neighbors. 62 and ¢e 
give the direction of the applied electric field 
with respect to the same axis. The evaluation of 
m follows without difficulty using the expansion 
exp (x) =1+x+2°/2+--- and the identity cos 612 
=cos 6; cos Oe+sin 6; sin @2(sin sin 
+cos ¢1 cos ¢2). The result is 


1+ cos (17) 
2 
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Higher powers of /kT are readily calculable, 
but are in general negligible. 

The values for 62 and wu, as a function of 
6, are tabulated in Table VIII and shown in 
Fig. 2. Here uw, is the dipole moment which is 
calculated from measurements on the gas, while 
ui is the value calculated in the same way if the 
rotation is restricted to lie between the polar 
angles 0 to 

If in addition to moving freely between 0 and 
6, an extra potential energy E permits the dipole 
to move between 6; and 180°, the equation cor- 
responding to (16) gives, 


sin? 6,(1—exp (—E/RT)) 


R= 


Finally if the potential energy for rotation between 0° and @,° is zero, between 6; and 4, is 


3k 2(1—cos 6,+(1+cos exp (— E/RT)) 


E, and between 62 and 180° is again zero, then the analog to (16) gives 


m= 


Other cases can be obtained in a similar way but 
these will suffice for our purpose. 

Consider first the application of the theory to 
a solution of a polar substance in a nonpolar 
solvent. If the solution is sufficiently dilute for 
interactions between the dipoles to be negligible, 
any restriction on the free rotation of the polar 
molecule is imposed by the nonpolar molecules 
around it. Under these conditions, the polar 
molecules can orient in the solution with almost 
exactly the same ease that they can do so in the 
gas phase, provided the frequency of the applied 
field is well below the natural frequency of rotation 
of the solute molecules in the solution. For small 
solute molecules in relatively nonviscous sol- 
vents, the frequencies at which derivations due 
to dispersion can be expected are well above the 
present experimental range. 

In pure polar liquids, however, an entirely dif- 
ferent situation exists. Under these conditions all 
the molecules are not free to orient in the field 
both because of dipole coupling and because of 
the size and shape of the molecules. As we have 
already seen for the case of water, spectroscopic 


( (sin? 6;—sin? 62)(1—exp (— E/RT)) )] 
3k 


2+(cos 62—cos 6,)(1—exp (— E/RT)) 


considerations and the velocity of sound, in con- 
junction with the vapor pressures, agree in 
giving a value for 6; of about 0.04, corresponding 
to a d2 of 0.12(=0.04!). We use the familiar 
equations, 


e—1 we 
(20) 
e+2 3 3kT 
TABLE VIII 
1 0° 10° 20° 30° 40° 50° | 60° 
| .0153 | .0594| .1296 | .2202 | .3252 4374 
bo 0 | .0076 | .0301| .0665 | .1170 | .1786 | .2500 
| 
01 70° 80° 90° 100° 110° 120° 
uu? | 5508 | .6555 | .7500 | .8292 | .8919 | .9375 
62 3290 | .4131 | .5000 | .5869 | .6710 | .7500 
1 130° 140° 150° 160° 170° 180° 
rug? | .9681 | .9864 | .9954 | .9990 | 1.0000 1.0000 
be 8214 | .8830 | .9335 | .9699 | .9924 | 1.0000 
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n?>—1 4n 
V=—WNa. 
n?>+2 3 


and (21) 


At 20° and one atmosphere water has an index 
of refraction for sodium light of equal to 1.3330, 
a molal volume, V, of 18.05 cc and a dielectric 
constant, «, of 80. This leads to a value for 
(42/3)Na of 3.7 cc and = 13.7 cc 
for the liquid. For the gas (47N/3)y?/3kT =71 
cc“? From Fig. 2 we see that this value of 
urn, of 0.194 leads to a value of 62 of 0.10 to 
be compared to 0.12 found above. 

The same calculations for methyl alcohol at 
0°C yield the following results. V;63 from vapor 
pressures and heats of vaporization is 0.026 cc, 
and V; from sound velocity is 0.54 cc. The value 
of 63 from these results is 0.048, corresponding 
to 52 equal to 0.13. Utilizing the data of Akerlaf* 
we obtain the value of 0.25 for 62 at the same 
temperature from the dielectric constant and 
dipole moment of the molecule. It is of interest 
to observe that 62, as determined from dielectric 
constant, increases with the temperature being 
0.34 to 60°C. Such a result is, of course, to be 
expected. 

There is a tendency for 62 from the dielectric 
constant to be greater than (63)! determined 
from other sources. This might be interpreted as 
due to a restriction on rotation about the axis in 
which the electric moment of the molecule lies 
which is greater than that about the other two 
axes. This is probably a valid argument in the 
case of water since the higher of the two fre- 
quencies which were employed to account for 
restricted rotation undoubtedly corresponds to 
a mode of motion involving only the hydrogens, 
and thus is a libration about the axis in which 
the dipole lies. Computing this correction from 
the Einstein functions corresponding to these 
frequencies gives a value of 62 of 0.09 to be 
compared to 0.10 from the dielectric constant. 


* P. Debye, Chem. Rev. 19, 171 (1936). 
% Akerléf, J. Am. Chem. Soc. 54, 4129 (1932). 
*® Trans. Faraday Soc. 30, 905 (1934). 
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It is not unreasonable that the discrepancy in 
the case of methyl alcohol might be accounted 
for on a similar basis, but the data for making 
such a computation are not at present available. 

However, various factors which might cause 
deviations from the theory have been neglected. 
Thus, the coupling of dipoles in the liquid may 
operate in such a way as to either increase or 
reduce the apparent moment of the molecule in 
the liquid. Also the adequacy of equation (20) 
for pure polar liquids is by no means always con- 
ceded, even after account is taken of the 
restricted rotation of the molecules.” 

However, it appears that any explanation of 
the deviations from (20) which is not compatible 
with restricted rotation of the molecules in many 
liquids is not in agreement with other lines of 
evidence. Fowler,*! and Debye and Ramm*® have 
also interpreted the dielectric properties of con- 
densed phases from the point of view of restricted 
rotation. 


CONCLUSION 


The results of this work lead us to suggest 
that the abnormalities of the so-called associated 
liquids can be explained entirely from the stand- 
point of restricted rotation of the molecules. 
Thus although water and methyl! alcohol are 
ordinarily classed as highly abnormal, their free 
volumes from sound velocity are not markedly 
different from those of other liquids. It is only 
the product, V;é3, rather than V; itself, which 
is a different order of magnitude for “hydrogen 
bonded” as compared to normal liquids. The 
values recorded here for 6; for water and methyl 
alcohol thus account for their abnormally high 
entropies of vaporization. 


4° Onsager, J. Am. Chem. Soc. 58, 1486 (1936). Van Vleck, 
J. Chem. Phys. 5, 556 (1937). J. Wyman, J. Am. Chem. 
Soc. 58, 1482 (1936). R. H. Cole, J. Chem. Phys. 6, 385 
(1938). S. S. Kurtz, Jr., and A. L. Ward, J. Frank. Inst. 
224, 583 (1937). 

41R. H. Fowler, Proc. Roy. Soc. 149A, 1 (1935). 

4 P, Debye, Phys. Zeits. 36, 193 (1935); P. Debye, and 
W. Ramn, Ann. d. Physik 28, 28 (1937). 


ib] 
n of 
nin 
h i 
hi 
the 4 
olar 
and 
pole 
cor- 
(18) 
: 
is 
(19) 
— 
con- 
ding 
iliar 
(20) 
pe: 
60° 
1373 
500 
20° 
00° 
000 


OCTOBER, 1938 


JOURNAL OF CHEMICAL PHYSICS 


Order Versus Disorder in Ternary Structures Including Certain Spinels' 


F. C. BLAKE 


VOLUME 6 


Ohio State University, Columbus, Ohio 


(Received June 27, 1938) 


The doubling process as the second alternative suggested by Professor M. v. Laue for mix- 


crystals is followed through in the case of lithium ferrite and it is shown in that case that all 
superlattice lines disappear when the lattice edge is eight times as large as the small disordered 
lattice found by Posnjak and Barth. It is shown that for this large (33A) lattice the simple 
space group, O,", gives only face-centered lines on the powder photograph. The disappearance 
of the superlattice lines for this lattice may be said to be due to the coherence of the x-rays 


that cause the destructive interference involved. 


1. INTRODUCTION 


ROFESSOR W. L. Bragg in his Bakerian 

lecture, June 28, 1934, in conjunction with 
Dr. E. J. Williams,’ discusses order versus dis- 
order in alloys, and they call attention to some 
of the earlier work of others in this field. In 
general the effect of increasing temperature on 
alloys is to change an ordered lattice into a dis- 
ordered one, and one may expect a sharp dis- 
continuity in the specific heat of an alloy at the 
critical temperature as the temperature is 
lowered across it. With further lowering of the 
temperature the change in specific heat is more 
gradual as the degree of order is increased. 

In his paper on the ‘Statistical Theory of 
Superlattices” Professor H. A. Bethe® extended 
Bragg and Williams’ definition of order so as to 
include local order; thus he defines long range 
and short range order, and calculates these for 
the simple cubic lattice; but interaction between 
nearest neighbors is alone considered. Chang* 
has extended Bethe’s theory to include inter- 
action between next nearest neighbors and gets 
a better agreement with experiment.5 He shows 
that with increasing interaction between second 
neighbors, the value of the specific heat at the 
critical temperature is very sensitive to the exact 


1 Read before the American Physical Society, Indian- 
apolis, Dec. 1937. 

? Bragg and Williams, Proc. Roy. Soc. A145, 699 (1934). 
See also their second and third papers, Proc. Roy. Soc. 
A151, 540 and A152, 231 (1935). Many more references to 
earlier work are to be found in their second paper; the 
third paper is by Williams alone. 

3 Bethe, Proc. Roy. Soc. A150, 552 (1935). 

‘Chang, Proc. Roy. Soc. Al61, 546 (1937); see also 
Bethe, J. App. Phys. 9, 244 (1938). 

5 Sykes, Proc. Roy. Soc. A148, 443 (1935). 


value chosen for this second neighbor interaction. 
He states, however, that it is not expected that 
his treatment should agree with experiment too 
well, for it leaves the higher interactions and the 
heat motions out of account. More recently 
Kirkwood® has presented “‘an alternative theory 
of order and disorder, essentially equivalent to 
Bethe’s in its degree of approximation.”’ Bethe 
in a later paper’ adopts the method of Kirkwood. 

Order-disorder problems in ternary alloys, or 
crystal structures, have, so far as the author is 
aware, received very little attention in the 
literature. In studying Posnjak and Barth's‘ so- 
called principle of ‘‘variate atom equipoints” 
the writer was aware of the two points of view 
concerning “‘mix-crystals” discussed by M. von 
Laue? in explanation of the results of Vegard and 
Schjelderup’® for some mix-crystals of alkali 
halides. The two views are these. In a mix-crystal 
the two components, A and B say, are either 
distributed statistically (in a disordered fashion) 
or else there is perfect order, but the lattice is 
much larger. Vegard and Schjelderup chose the 
first alternative. It is well known that certain of 
the alkali halides form solid solutions in all pro- 
portions."" These alkali halide mix-crystals are 
ternary, not binary, structures, yet with regard 
to mixing they are binary. In KBr-KCI there is 
a common cation, K. In KI-NHg,I there is 4 
common anion, I. Both of these mixes form solid 
solutions in all proportions." 

6 Kirkwood, J. Chem. Phys. 6, 70 (1938). 

7 Bethe, J. App. Phys. 9, 244 (1938). 

8 Posnjak and Barth, Phys. Rev. 38, 2234 (1931). 

9 Laue, Ann. d. Physik 56, 497 (1918). - 

10 Vegard and Schjelderup, Physik. Zeits. 18, 93 ( 1917) 


1 E. g. Havighurst, Mack and Blake, J. Am. Chem. Soc. 
47, 29 (1925). 
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ORDER VS. DISORDER 


2. SUPERSTRUCTURES 


In binary alloys, when the degree of order is 
high, there is in general a tendency for super- 
lattice lines to appear, because of the two super- 
structures present. The difference in scattering 
power for x-rays of the two components present 
results in extra reflections that are not present 
when the disorder is complete. This is well illus- 
trated in the copper-gold alloy, Cu;Au, studied 
by Sykes and Evans.” For a copper-gold wire 
heated to 425°C and quenched in water, only 
eight face-centered lines appeared on the powder 
photograph. After ordering by various methods 
of cooling, twelve extra superlattice lines appear, 
all of which are simple cubic. Techniques were 
developed by them for changing the character 
of the superlattice lines from fuzzy to sharp. 


3. SECOND ALTERNATIVE FOR MIXED CRYSTALS 


Since M. v. Laue® derived his expression for 
the structure factor of the mix-crystal K2BrCl, 
most people have chosen the first (disordered) 
alternative in studying binary or ternary alloys, 
or structures, for oftentimes superlattice lines are 
absent, and such lines are at once eliminated by 
assuming disorder. Where possible, it is much 
more scientific to study specific heats of strugtures 
as a function of temperature, as well as to make an 
x-ray analysis of the structures. But super- 
lattice lines can be eliminated by the second 
alternative, viz., that of choosing a larger lattice. 
When properly done, it is possible, in some cases 
at least, to obtain the exact equivalent of the 


TABLE I.* 


OBSERVED CALCULATED INTENSITIES (33A LATTICE) 
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*For an explanation of the symbols ring in the intensit 
formulas see Blake, Rev Mod Ph 5 “169 (19.33 ant con 


* Sykes and Evans, J. Inst. Metals 58, 255 ( 1936). 
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disordered lattice; and even though crystal- 
lographers and metallurgists might say that this 
large equivalent of the small disordered lattice 
could never be produced owing to an infinite 
time being needed for its formation, it is worthy 
of having attention called to it. 


4. Low TEMPERATURE MODIFICATION OF 
LITHIUM FERRITE 


A couple of years ago I secured from Dr. Barth 
of the Geophysical Laboratory some of the 
powdered lithium ferrite and zinc stannate that 
he and Dr. Posnjak had used in their spinel 
studies® and that led them to announce their 
principle of “variate atom equipoints.’’ Powder 
photographs of these two compounds were taken 
with molybdenum Ka-rays, and the photographs 
were densitometered. Proper account was taken 
of all of the intensity factors and of the laws of 
photographic blackening for x-rays. The results 
for lithium ferrite are shown in Table I, together 
with the observed results of Barth and Posnjak. 
The chief differences between my results and 
theirs is due to their overestimation of the in- 
tensity for faces having large values of @ for 
fairly long times of exposure. When the exposure 
time is long (i.e., S>0.6 for the strongest line on 
the film) there is an exhaustion of silver for lines 
having @ small as against those having @ large, 
thus leading to errors. 

In studying this second alternative a frontal 
attack was made on the problem by doubling the 
lithium ferrite lattice from 4.141A to-8.282A, 
care being used to increase the symmetry of the 
crystal if possible as the size was doubled. It is 
always possible in a crystal lattice to move the 
origin of coordinates without affecting inten- 
sity relations, though symmetry properties are 
thereby affected. In all cases I have chosen an 
oxygen atom as origin. Allowing no exchange 
between iron and lithium atoms, one is, of course, 
conscious in the 8.28A cell that neighboring 
4.14A parts are not identical. Indeed, if one 
alternates lithium and iron in between the oxygen 
atoms in all directions parallel to the axes, he is 
conscious that when he has progressed along the 


13 (a) Barth and Posnjak, J. Wash. Acad. Sci. 21, 255 
(1931); (b) Posnjak and Barth, Phys. Rev. 38, 2234 (1931); 
(c) Zeits. f. Krist. 82, 325 (1932). 

4 Blake, Rev. Mod. Phys. 5, 169 (1933). 
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Fic. 1. Calculated intensities for 8.28A lattice. 


cubical axes 8.28A the structure he has thus 
built up is only a piece of a cell, for at one corner 
he finds an oxygen atom with three iron atoms as 
nearest neighbors while at the opposite corner 
(that furthest removed) an oxygen atom is sur- 
rounded by three lithium atoms as nearest 
neighbors, the oxygen atoms at the other four 
corners having two lithium and one iron as 
nearest neighbors, or vice versa. Thus this 8.28A 
lattice is only one-eighth of a proper unit cell, 
and it is not surprising that when one plots 
intensity calculations against sin 6 he gets a great 
overstressing of the superlattice (dotted) lines 
(faces), the full lines in Fig. 1 being the calcu- 
lated face-centered lines that are observed, where 
for the calculations all intensity factors are taken 
into account except the temperature factor and 
the absorption factor. The first legitimate doub- 
ling of the lattice enlarges it at once to 16.56A. 
For, by inserting axial planes of symmetry, each 
oxygen atom at the corners of the enlarged lattice 
is surrounded by three iron atoms as nearest 
neighbors, while the oxygen atom at the center 
is surrounded by six lithium atoms as nearest 
neighbors. This fact determines the space-group 
of the lattice, viz., simple cubic O',. When the 
intensity relations are calculated and plotted 
(Fig. 2), it is noticed that the superlattice 
(dotted) lines are reduced in intensity, being 
largely crowded down into the region where @ is 
small. But since the powder photographs show 
their absence, another lattice doubling (the 
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Fic. 2. Calculated intensities for 16.56A lattice. 


second legitimate one) is suggested. Thus the 
correct, ordered lattice, that is, the exact equiva- 
lent of the statistical averaging of Posnjak and 
Barth’s 4.141A lithium ferrite lattice, has eight 
times the edge of that lattice, viz., 33.13A. The 
result is shown in Fig. 3. The only superlattice 
line residue is for face 333, and it has an intensity 
less than 4 percent and hence could not be de- 
tected, for it occurs in the region of maximum 
general blackening. Moreover in Posnjak and 
Barth's observations no line is reported whose 
intensity is less than ten percent. In what I call 
the second legitimate doubling it should be 
borne in mind that care was used to preserve the 
axial planes of symmetry. This second doubling 
again gives the space-group as O!,. Thus an 
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TABLE II. Coordinates of oxygen atoms. 


TABLE III. Coordinates of iron atoms. 
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containing 512 molecules of LigFexO,, with ca- 
tions differing in scattering power as 23 to 2 (or 
as 26 to 3 if neutral atoms) gives only face- 
centered reflection intensities. I invite the atten- 
tion of the crystallographers to this fact. I think 
this gives the answer also to the question implied 
in this sentence of Posnjak and Barth’s paper on 
lithium ferrite, “Taking into consideration not 
too small a portion of the crystal, an equal 
number of points is occupied by iron and by 
lithium.” The “‘not too small a portion” is a cube 
of 33A edge. 

In order that others who care to may study 
this 33A lattice for themselves, I give in Tables 
II, II and IV the atomic positions of the 
oxygen, iron and lithium atoms, respectively, 
with certain necessary explanations. In Table II 
the figures 0 and 1 occur while in Tables III and 
IV the figures 0, 1, 2, 3 occur. Let any figure in 
any of the three tables represent the numerator 
of a fraction of which the denominator is 16. To 
illustrate, choose the figure 3, then 3/16 is the y 
coordinate of an atom whose x and z coordinates 
are given (also as numerators of fractions having 
the denominator 16) at the top and side of each 


table, where it is understood that 7/16 is the 
same as 9/16, because the unit angle chosen for 
the revolving radius vector!’ is 27/16=22.5 de- 
grees. Thus the least count in stating the exact 
position of an atom along any of the three axes 
is 1/16 of the lattice edge. One interprets the 
tables correctly thus. 

Since there are 16 units along the x and z axes, 
there must be 16 units along the y axis. In Table 
II only two such units are given by the figures 0 
and 1 as numerators of y for positions occupied 
by 256 oxygen atoms. To obtain positions of 
another 256 of the total 2048 oxygen atoms in the 
unit cell, it is necessary to substitute 2 for 0, and 
3 for 1, in the y values of Table II. Similar sub- 
stitutions of 4, 6, 8, 6, 4, and 2 for 0, and of 5, 
7, 7, 5, 3, and I, for 1, give the coordinates of 
the remaining oxygen atoms. Table III gives the 
coordinates of 256 iron atoms. The remaining 
iron atom positions, which total 1024 in the unit 
cell, are obtained by substitutions of 4, 8, and 
4 for 0, of 5, 5, and I for 1, of 6, 6, and 2 for 2, 
and of 7, 7 and 3 for 3, in the y values of Table 
III. The same substitutions made in Table IV 


% See Blake, J. Chem. Phys. 2, 320 (1934). 
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TABLE IV. Coordinates of lithium atoms. 
fo] 7/8] 7/6] 5] 4/3 
1 4 
2 2 
3 
4 4 
s 
© 6 
7 7 
8 
7 
6 6 
5 5 
4 
3 3 
2 3 
1 1 
0 


obviously give the positions of the 1024 lithium 
atoms in the unit cell. To obtain the positions 
of all atoms in any one plane of the unit cell, e.g. 
y=1/16, take the positions of the figures 1 in 
Table II as the positions of oxygen atoms in the 
plane y=1/16, the positions of figures 1 in 
Table III as the positions of iron atoms in that 
plane, and the positions of the figures 1 in Table 
IV as the positions of lithium atoms in the same 
plane. 

A few words as to the space group, O'; or 
Pm3m. In calculating the structure factor the 
general formulas give A= cos 2zhx[cos 2rky 
-cos 27lz+cos 2zly-cos 2rkz], and B=0. But I 
have found it no easier to use these formulas 
than to summarize by the revolving vector 
method.” As already stated, assuming that for 
the spinels the absorption and temperature fac- 
tors pretty well cancel each other, the correct 
intensity relations are given in Fig. 3 and these 
agree well enough with the lithium ferrite obser- 
vations of Table I. The atomic form factors for 
lithium, iron and oxygen used are those of the 
International Tables for the Determination of 
Crystal Structures, the crystal being taken as 
ionic or atomic as indicated in Table I. 


BLAKE 


In the case of zine stannate, I have tried the 
doubling process on the 8.61A disordered lattice 
given by Barth and Posnjak' to see if I could 
get the large lattice equivalent; but I haven't 
completed the study for the reason that if it is 
also true in this case that the lattice edge should 
be enlarged eight times, the labor involved be- 
comes quite burdensome, for the cell would then 
have an edge 68.88A. The 34.44A cell containing 
512 molecules of Zn2SnO, goes to zero for the 
fourth order of 311, no matter what value is used 
for the oxygen parameter. This indicates another 
doubling, or some sort of reorganization of the 
cell. I intend to study the doubling process 
further in the case of zinc stannate. 


5. LARGE-SIZED LATTICES AND MOLECULES 


It is becoming more fashionable to speak of 
large-sized lattices, even if in some cases there 
is apparently but one molecule per unit cell. The 
recent series of papers in WNature,'® ‘Crystal 
Structure of the Proteins,”’ gives a tetragonal 
lattice, a=b=63.5A, c=145A, and an ortho- 
rhombic lattice, a=49.6A, b=67.8A, c=66.5A, 
Bernal, Fankuchen and Perutz stating that ‘‘the 
proteins fall into two groups, pseudo-hexagonal 
and pseudo-cubic, in both cases with a side of 
60-75A.”’ Long needle-like structures up to a 
length of 6000A seem to have been prepared by 
Astbury and Bell.'7 What is the distinction in 
building blocks between the spinels and these 
organic compounds that would prevent 33A 
cubes for the lithium ferrite? In Sykes and Evans 
paper, ‘The Transformation in the Copper-Gold 
Alloy, Cu3Au,’’” the authors conclude that when 
this alloy was cooled at 30°C per hour it consisted 
of small ‘“‘nuclei’’ within the individual crystals: 
so they distinguish between two different states 
of the alloy, one with nuclei, the other without. 
Each nucleus could be said to have 100 percent 
local order but an aggregation of nuclei could 
lack long range order. It seems to me that the 
33A lithium ferrite lattice discussed in this paper 
shows 100 percent local order and 100 percent 
long range order, that is, the nuclei are properly 


fitted together. A recent remark of Professor H. 


16Crowfoot and Riley; Crowfoot and Fankuchen: 
Sane, Fankuchen and Perutz; Nature 141, 521-524 
(1938). 

17 Astbury and Bell, Nature 141, 747 (1938). 
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A. Bethe to the writer, that there ought to be 
more than one way in which all superlattice lines 
could be eliminated for a lattice as large as 33A, 
leads me to admit the possibility, but not the 
probability, of this unless the local order is in a 
nucleus as small as 4.14A, when the long range 
disorder is 100 percent. It should be borne in 
mind that if one turns the 33A lattice inside out, 
as it were, that means nothing but an interchange 
throughout the lattice of iron and lithium atoms, 
thus giving the same, not a different, solution. 
Granted that this 33A, ordered lithium ferrite, 
according to present notions concerning alloys, 
could not be produced in the laboratory because 
of the inordinate amount of time required to set 
it up, it nevertheless remains as the low tempera- 
ture form of the disordered 4.14A lattice of 
lithium ferrite. Further studies on rates of cooling 
of alloys with specific heat measurements and 
with the possible use of “‘stabilizing agents’’ may 
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help to clarify the knotty problems involved in 
order versus disorder studies. Bragg and Williams? 
in their first study in this field remark as their 
concluding sentence, ‘‘Maxima and minima in 
physical properties at certain relative propor- 
tions (e.g. Fe;Al and Cu;Au) are statistical 
effects, and do not imply the existence of cor- 
responding compounds.”’ It seems conceivable to 
the writer that an ordered solution of such binary 
alloys could be found, in which case it would be 
a fair question to ask whether a single ordered- 
lattice alloy is not a definite compound. It has 
been shown that proper ordering causes super- 
lattice lines to disappear, that a large (33A) 
lattice of lithium ferrite is the exact equivalent of 
the small (4.14A) disordered lattice, that a simple 
cubic lattice of edge 8 times the disordered face- 
centered cubic lattice gives only face-centered 
lines on the powder photograph in the case of 
lithium ferrite. 
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The average dipole moment of a molecule containing two movable dipoles has been calculated 


as a function of temperature and restriction to internal rotation for three different potential 
functions. The potential functions used are: (V»/2)(cos ¢+1), (Vo/m)(#—¢) and Vo/x?(r—@)?, 
where V» is the difference in potential energy between the cis and trans forms of the molecule 
and ¢ is the angle between the planes containing the dipoles and the axis of rotation. Values of 
the integral Jo” cos x exp (—bx*)dx are tabulated for values of b ranging from 0 to 2 at intervals 
of 0.1. The dipole moments of s-dichloro-ethane, 1,2-bromo-chloro-ethane and s-dibromo-ethane 
as a function of temperature are used with the three potential functions to determine values of 


INTRODUCTION 


"[TRERE have recently been several attempts 

to determine the magnitude of the po- 
tential barrier restricting internal rotation around 
single carbon-carbon bonds. Most of these at- 
tempts have been made on ethane.! In attacking 


National Research Fellow in Chemistry. 

Kemp and Pitzer, J. Chem. Phys. 4, 749 (1936); J. Am. 
Chem. Soc. 59, 276 (1937). Howard, Phys. Rev. 51, 53 
(1937). Kistiakowsky and Nazmi, J. Chem. Phys. 6, 18 
(1935). Bartholomé and Karweil, Zeits. f. physik. Chemie 
ao Hunsman, Zeits. f. physik. Chemie 39, 23 


V». The values thus obtained are compared with the values obtained by other methods. 


the problem from the standpoint of dipole 
moments it is necessary to substitute for sym- 


metrically located hydrogen atoms, polar groups 
whose relative orientations will have a large 
effect on the resultant dipole moment of the 
molecule. In order to treat successfully the dipole 
moment data of such molecules as the 1,2,di- 
substituted ethanes it is necessary to derive ex- 
pressions relating the observed dipole moment 
to the height and shape of the potential barrier 


and the individual bond moments. We have de- 
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rived such expressions for the cases of cosine, 
linear and parabolic potential functions. To find 
the average dipole moment it is necessary to 
calculate the average value of the square of the 
dipole moment, averaging over all possible con- 
figurations of the molecule and weighting each 
configuration with a Boltzmann factor taken from 
the assumed potential function. The square of 
the moment rather than the moment itself, is 
averaged because the Debye? equation applied 
to each molecular species gives the square of the 
moment, rather than the first power of the 
moment, as a function of the dielectric constant 
and the refractivity. 


CALCULATIONS 


The square of the dipole moment of a molecule 
containing two dipoles capable of relative motion 
to the extent of rotating about an axis in the 
molecule is: 


u(¢)=A+B cos ¢, (1) 


where COS a; COS a2, and 
B=2mymz sin a sin aa. m, and mz are the mag- 
nitudes of the two movable dipoles. a; and ae 
are the angles made by the dipoles with the axis 
of rotation. ¢ is the angle between the plane con- 
taining one dipole and the axis of rotation and 
the plane containing the other dipole and the 
axis of rotation. For a; and ae tetrahedral and 
m,=m2=m, A and B each become 16/9 m?. 

The observed, or average dipole moment for 
the molecule is 


exp (— V(¢)/RT)de 


exp (— V($)/RT)d¢ 


0 
TABLE I. p for linear and cosine potential functions. 


a (LINEAR) (COSINE)| @ (LINEAR) (COSINE) 
0.0 0.000 0.000 | 4.0 0.868 0.863 
0.5 0.199 0.242 4.5 0.892 0.880 
1.0 0.378 0.446 5.0 0.910 0.892 
1.5 0.527 0.596 6.0 0.935 0.912 
2.0 0.640 0.698 7.0 0.952 0.927 
2.5 0.727 0.765 8.0 0.963 0.935 
3.0 0.788 0.808 9.0 0.970 0.944 
3.5 0.835 0.841 | 10.0 0.977 0.954 


2 Debye, Polar Molecules. 
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The effect of the angle ¢ on the moments of 
inertia and thus the entropy of the various con- 
figurations is negligible except for small ¢,* and 
will be neglected. Integrating ¢ from 0 to z is 
equivalent to integrating from 0 to 27. V(4¢) is 
only required to represent the potential function 
when ¢ is between 0 and z. Substituting Eq. (1) 
into Eq. (2) we obtain 


=A+Bp (3) 


cos exp (— V(¢)/RT)d¢ 


0 


where p= 
f exp (— V(¢)/RT)dé 
0 
We shall first evaluate p for the linear po- 
tential function 
V(¢) 


Vo is the difference in potential energy between 
the cis and trans configurations of the molecule. 
In this case 


cos ¢ exp (— do 


p= 
[exp 
0 
= coth a, 
where a=V,/2RT. p=0 when a=0 and ap- 
3 Altar, J. Chem. Phys. 3, 460 (1935). 
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proaches —1 as a becomes large. When a is 
large and p= —1, the dipole moment (A —B)', 
which is the moment of the ‘rans molecule. When 
a is small, the moment approaches A? the value 
for free rotation. p is tabulated in Table I and 
plotted in Fig. 1, curve A, as a function of a. 
The integrals arising when the potential func- 
tion V=(Vo/2)(cos ¢+1) is used have been 
evaluated by Mizushima and Higasi.‘ In this case 


| cos exp (— Vo(cos ¢+1)/2RT)d@ 
0 


{ exp (— Vo(cos 
0 


f cos ¢ exp (—a cos ¢)d¢ 
0 


exp (—a cos 
0 


=1J,(1a)/Jo(ta). a=Vo/2RT. 


Numerator and denominator are Bessel functions 
which are tabulated.* In this case too p=0 for 
a=(0and —1 for a large. p is tabulated in Table I 
and plotted in Fig. 1, curve B. 

In calculating the average moment for the 
parabolic potential function V =(V 


cos ¢ exp (— do 
0 


p= 


exp (— do 


Let b= Vo/?RT and r—¢=x, then 


f cos x exp (—bx?)dx 
0 


p= — 


exp (—bx?)dx 
0 


The integral B(b) can easily be transformed 
into an error integral. The result is B(b) =(3) 


‘ Mizushima and Higasi, Proc. Imp. Acad. 8, 482 (1932). 
Jahnke and Einde, Tables. 


X is the normalized error 
integral and is tabulated.°® 

The integral S(b) does not become integrable 
by simple transformations, nor is it tabulated. 
We found it necessary to expand both cos x and 
exp (—bx*) in powers of x and integrate term by 
term. The result is the following double series: 


ton (21-4+1)(21—2n)! 


5S(b) => cub", Ca 
n=0 


This series converges fairly rapidly for small 
values of 6 and was used for calculating S(d) 
when 6=0.5 or less. When 3 is large S(b) can be 
approximated to by the integral 


cos x exp (—bx?)dx = }(2/b)' exp (— 4d). 
0 


We can test the error made by this approxima- 
tion in the following way. The error is less than 
the error which would be made by making the 
same approximation for B(b) because cos x is 
always less than 1. The error in B(d) is found, 
from the error function table, to be 0.1 percent 
at b=0.55, and to decrease rapidly as } increases. 
For calculating S(b) when b=0.6 or larger the 
approximate formula was used. S(b) and p are 
tabulated in Table II. In Fig. 1, curve C p is 
plotted as a function of a. p varies from 0 to —1 
as b (or a) varies from 0 to ~. 


TABLE II. S(b) and p for the parabolic potential function. 


S(b) 


0.0000 0.0000 
0.2606 0.0968 
0.4373 0.1858 
0.6367 0.3371 
0.7238 0.4541 
0.7563 0.5425 
0.7622 0.6097 
0.754 0.659 
0.741 0.700 
0.731 
0.757 
0.779 
0.796 
0.813 
0.825 
0.836 
0.846 
0.855 
0.863 
0.870 
0.876 
0.881 
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TABLE III. Potential barrier, Vo. 


POTENTIAL FUNCTION 
COSINE LINEAR PARABOLIC 


T°K X 10718 KCAL. / KCAL./MOLE KCAL. / MOLE 
C.H,Cl. 
305 1.12 4.02 4.17 7.97 
341 1.24 3.78 4.15 711 
376 1.32 3.71 4.46 6.85 
419 1.40 3.71 4.26 6.63 
457 1.45 3.80 4.41 6.66 
480 1.48 3.81 4.52 6.66 
485 1.48 3.85 4.57 6.73 
544 1.54 4.02 4.80 6.85 
C.H,BrCl 
339 1.09 4.65 4.77 9.32 
368 1.14 4.67 4.90 9.22 
405 1.20 4.74 5.12 9.10 
436 1.28 4.54 5.05 8.46 
C.H aBre 

339 0.94 6.00 5.66 12.87 
368 0.99 5.92 5.74 12.20 
405 1.03 6.12 6.04 12.56 
436 1.10 5.85 5.92 11.69 


Moment of CH.Cl=2.037 

Moment of CH.Br =2.027 

™ of C2H4BrCl and from 
ahn® 


DISCUSSION 


If one of the potential functions described 
above is assumed to be the actual one, V» can 
be found from the value of the dipole moment at 
a single temperature. Substituting in Eq. (3) 
values for the group moments, bond angles and 
the observed moment of the molecule, a value of 
p is obtained. From the graph of the relation 
between p and a, a value of a is obtained which, 
with the temperature at which the dipole moment 
was measured gives a value of Vo. If values of 
the dipole moment have been determined for 
several temperatures several values of Vo can be 
obtained. These values should be consistent if 
the assumed potential function is a good one. 

We have taken the dipole moment data for 
s-dichloro-ethane, 1,2,bromo-chloro-ethane and 
s-dibromo-ethane® and calculated values for Vo 
for each of the potential functions. We assumed 
that a1=a2= 109° 28’ for all three molecules. For 
the group moments we have used the gas values 
of the moments of ethyl chloride and bromide.’ 
The ethyl rather than the methyl compounds 


6 Zahn, Phys. Rev. 40, 291 (1932); 38, 521 (1931). 
7 Smyth and McAlpine, J. Chem. Phys. 2, 499 (1934). 
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were used to obtain the group moments in order 
that the inductive effects should be approxi- 
mately the same as in the disubstituted ethanes, 
Gas values of the dipole moments were used for 
the three disubstituted ethanes. The values of /, 
calculated by the above method are tabulated in 
Table III. Professor C. P. Smyth has suggested 
to us that perhaps the group moments should be 
obtained from the methyl rather than the ethyl 
compounds. If this is done, the values in Table 
III are all lowered by about 1 kcal./mole. The 
arguments in favor of either procedure are not 
conclusive. 

The values of V» calculated from the cosine 
potential function are more consistent than the 
values obtained from the linear or parabolic 
potential functions. They are also generally 
lower. In estimating a probable value for Vy we 
shall accordingly give them more weight. From 
these results the probable values of Vo for 
s-dichloro-ethane, 1,2,bromo-chloro-ethane and 
s-dibromo-ethane are 4, 5, and 6 kcal./mole, 
respectively. When values of the dipole moments 
obtained in solution are used somewhat lower 
values of Vp» are obtained. Lennard-Jones and 
Pike® have treated dichloro-ethane using the 
cosine potential function. They find that lower 
values of Vo are obtained, 2.0 to 3.2 kcal./mole, 
if a variation of the group moments with ¢ is 
considered. Altar* has treated the dichloroethane 
dipole moment data using a two-term Fourier 
series to represent the potential function. This 
potential curve extended smoothly to small 
values of ¢ gives a value of Vo of about 5 kcal./ 
mole. An extrapolation to ¢=0 is necessary as 
his potential function does not represent the 
actual potential function at small values of ¢. 
The determination of the potential curve near 
¢=0 by any experimental method is doubtful 
because of the small numbers of molecules it 
those configurations. The shape of the curve 
when ¢ is greater than 90° (say) is of much more 
significance. 

The values of V» obtained in this paper are in 
fairly good agreement with electron diffraction 
results on these compounds.°® 


and Pike, Trans. Faraday Soc. 30, 830 
1934). 

9 Beach and Palmer, in publication; Beach and Turke- 
vich, to be published. 
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Internal Rotation in Ethylene Chloride* 


J. Y. Beacut anv K. J. PALMER 
Princeton University and California Institute of Technology 


(Received July 1, 1938) 


The electron diffraction method of investigation has been used to determine the nature of the 
potential barrier restricting the rotation about the carbon-carbon single bond in ethylene 
chloride (CH2Cl—CH,.Cl). Assuming that the shape of the barrier is given by a cosine function, 
the height of the barrier is between 5 kcal./mole and infinity. This potential barrier is compared 
with the barrier resulting from the exchange repulsions and the electrostatic repulsions of the 
two ends of the molecule. Eyring’s method was used for calculating the exchange repulsions. 
The repulsions between the chlorine atoms must be calculated differently from the repulsions 
between hydrogen and chlorine atoms in order to obtain agreement between theory and experi- 
ment. Duplicate experiments were done at Princeton and at Pasadena. 


HE electron diffraction photographs of 

ethylene chloride vapor from which these 
results are deduced were taken at Pasadena and 
at Princeton. Both compounds were Eastman 
products, fractionally distilled in three-foot 
columns on which the reflux ratio could be con- 
trolled. The first photographs were taken by 
both of us at Pasadena. The Pasadena apparatus 
has been described by Brockway.! The Princeton 
apparatus is similar to the Pasadena apparatus 
and has been used previously.2 The Pasadena 
and Princeton photographs were made with 
electrons having wave-lengths of 0.0611A and 
0.0590A, respectively. The two sets of photo- 
graphs, which are identical except for the effect 
of electron wave-length, and lead to identical 
results, will not be discussed separately. 

We did not attempt to determine by electron 
diffraction all the bond distances and bond 
angles in the molecule in addition to the shape 
of the potential barrier. The carbon-hydrogen 
distance was assumed to be 1.09A, the carbon- 
hydrogen distance in methane.? The carbon- 
carbon distance was assumed to be 1.54A, the 
carbon-carbon distance in diamond‘ and a large 
number of saturated hydrocarbons.* All the bond 
angles on the carbon atoms were assumed to be 


*Part of these results were presented to the American 
Physical Society at New York City, February, 1938. 
I National Research Fellow in Chemistry. 
: Brockway, Rev. Mod. Phys. 8, 231 (1936). 
Beach and Stevenson, J. Am. Chem. Soc. 60, 475 
(1938) ; J. Chem. Phys. 6, 75 (1938). 
, Ginsburg and Barker, J. Chem. Phys. 3, 668 (1935). 
Strukturbericht. 


og puling and Brockway, J. Am. Chem. Soc. 59, 1223 


tetrahedral. This assumption is made reasonable 
by the fact that in methylene chloride and 
chloroform, where the chlorine-chlorine distances 
are only about 2.90A, the bond angles are but 
slightly (about 3°) distorted from the tetra- 
hedral value.® 
The observed values of s(sons= (4m sin 0/2)/X) 
for the maxima and minima and the intensities 
of the maxima of the electron diffraction photo- 
graphs are tabulated in Table I. The heavier 
photographs show up to nine maxima. The last 
two cannot be measured with much accuracy. 
The first three maxima are of approximately 
equal intensity. The fourth and fifth are much 
weaker. The sixth is strong. The seventh, 
eighth and ninth are progressively weaker. 
The data have been treated by the radial 
distribution method? and by the modified radial 
distribution method.* Both radial distribution 
curves are shown in Fig. 1. The modified radial 
distribution curve differs from the ordinary 
radial distribution curve in that the peaks are 
sharper. The positions of the maxima are prac- 
tically identical. The average positions of the 
maxima on the two curves are 1.73A, 2.70A and 
4.29A. The first two peaks represent carbon- 
chlorine distances, and the third peak represents 
the chlorine-chlorine distance. 4.29A is just the 
value of the chlorine-chlorine distance for the 
trans-configuration of the molecule if the carbon- 


— and Brockway, J. Am. Chem. Soc. 57, 473 
(1935). 

7 a and Brockway, J. Am. Chem. Soc. 57, 2684 
(1935). 

8 Schomaker, to be published. 
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chlorine distance is assumed to be 1.76A, the 
radius table value® for the carbon-chlorine single 
bond distance. The dissymmetry in the chlorine- 
chlorine peak on the modified radial distribution 
curve indicates that some of the molecules are 
oscillating from the ¢rans-configuration to inter- 
mediate configurations having shorter chlorine- 


TABLE I. Ethylene chloride.* 


Max.| Min. I Sobs Sa Sa/Sobs Sb Sb/Sobs 


1 12 | 3.160} 3.26 | 1.032 | 3.31 | 1.047 
4.136| 3.97 | 0.960 | 4.01 | 0.970 
15 | 4.936| 4.84 | 0.981 | 4.87 | 0.987 
6.668} 6.68 | 1.002 
13 | 7.679) 7.72 | 1.005 |} 7.79 | 1.014 
8.594} 8.58 | 0.998 | 8.70 | 1.012 
4 | 9.194; 9.14 | 0.994 | 9.23 | 1.004 
9.964} 9.98 | 1.002 | 10.09 | 1.013 
3 {10.591 10.78 | 1.018 
8 
4 


11.302 | 11.12 | 0.984 
12.041 | 12.02 | 0.998 | 12.10 | 1.005 
13.183 | 13.02 | 0.988 | 13.14 | 0.997 
14.815 | 14.95 | 1.009 | 14.90 | 1.006 
Av. 0.998 Av. 1.005 
C-—Cl=1.76 *C—Cl=1.77 


Final results C—Cl=1.76+0.02A. 
5 kcal./mole == VoSinfinity. 


_ * (a) From trans curve. (b) From 6 kcal./mole curve. The first max- 
imum and the second minimum were omitted in the averaging. 


® Pauling and Huggins, Zeits. f. Krist. 87A, 205 (1934). 


Fic. 1. Curve A, modified radial distribution function. Curve B, ordinary radial 
distribution function. 
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chlorine distances. While doubtless such motion 
occurs, the radial distribution function is not 
sufficiently accurate to make a quantitative 
estimate of the extent of the motion possible. 

To arrive at more definite results regarding the 
oscillation about the ¢rans-configuration the elec- 
tron diffraction data have been treated by the 
visual method.!® The relation V= V (cos ¢+1) 
is assumed to represent the potential barrier to 
internal rotation. ¢ is the angle between the plane 
containing one chlorine atom and the two carbon 
atoms and the plane containing the other 
chlorine atom and the two carbon atoms. V» is 
the height of. the barrier and limits will be 
placed on its value from the electron diffraction 
data. This form for the potential function has 
been shown to be a fairly satisfactory approxi- 
mation to the actual potential function by 4 
study of the consistency of the Vo values ob- 
tained from dipole moments using different 
potential functions." 

The formula for the intensity of scattered 
electrons is 

sin sli 


sl ij 


10 Pauling and Brockway, J. Chem. Phys. 2, 867 (1934): 
1 Beach and Stevenson, J. Chem. Phys. 6, 635 (1938). 
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Fic. 2. Theoretical intensity curves. Curve A, cis model. Curve B, half-cis and 
half-trans. Curve C, free rotation. Curve D, Vo=0.6 kcal./mole. Curve E, Vo=6.0 
kcal./mole. Curve F, Vo=infinity, trans model. Curve G, segment of curve for 
potential function shown in Fig. 4, Curve A. 


l;; is the distance from the ith to the jth atom; 
Z is atomic number. Intensity curves were calcu- 
lated for ethylene chloride for values of Vs 
ranging from zero to infinity. When Vo equals 
infinity all the molecules are in the trans- 
Position. When Vo equals zero all the molecules 
rotate freely, and there are equal numbers of 
molecules in all configurations. For intermediate 
values of Vo the number of molecules possessing 
a given configuration is found from the Boltz- 


2 Altar, J. Chem. Phys. 3, 460 (1935). 


— 


mann factor, exp (—V(¢)/RT). The electron 
diffraction curves are calculated by adding the 
intensity formulas for molecules having values of 
o=0°, 30°, --- 180° weighted by the appro- 
priate Boltzmann factor. The Boltzmann factors 
and the final curves are determined by the value 
given to Vo. 

It has been shown by Altar” that this pro- 
cedure is not rigorously correct. The weighting 
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factors are not exactly equal to the Boltzmann 
factors. The entropy of the molecule varies with 
@ and the statistical weights are not constant 
when Vo equals zero. This is true for this type 
of molecule even if the product of inertia is 
independent of ¢. However, these effects are 
‘small at all angles near 180°, which is the region 


occupied by most of the molecules, and will be - 


neglected. 

The electron diffraction intensity curves are 
plotted in Fig. 2. The trans-curve corresponds 
most closely to the photographs. The 6 kcal./ 
mole curve is less satisfactory because the fifth 
maximum on it is not sufficiently prominent. 
The other curves are unsatisfactory from a 
qualitative inspection, none of them showing the 
fourth and fifth maxima correctly. Quantitative 
comparisons of the photographs with curves E 
and F are made in Table I. The curves in Fig. 2 
are all calculated for molecular models having a 
carbon-chlorine distance of 1.76A. This does not 
predetermine the carbon-chlorine distance be- 
cause if the value had been incorrect a value for 
the carbon-carbon distance different from the 
assumed 1.54A would have been obtained when 
quantitative comparisons of curves and photo- 
graphs were made. The final results are: C—Cl 
= 1.76+0.02A and 5 kcal./mole = VoSinfinity. 

Wierl," from electron diffraction photographs 
showing only three maxima, reported that the 
chlorine-chlorine distance in ethylene chloride is 
4.5A. He also stated that ethylene chloride could 
be an approximately equimolal mixture of cis 
and trans-molecules. The theoretical intensity 
curve for this mixture is curve B, Fig. 2. It is 
unsatisfactory because it possesses only one 


TABLE II. Repulsions in ethylene chloride, kcal./mole. 


¢ 0° | 30° | 60° | 90° | 120° | 150° | 180° 


Electrostatic 
interactions 1.80} 1.46] 0.85} 0.18] —0.29 | —0.59 | —0.65 


Exchange repulsions 14.03) 13.30) 12.83]13.83] 14.58 | 14.38 | 13.75 
Total 15.83) 14.76) 13.68}14.01} 14.29 | 13.79] 13.10 
E(¢) —E(x) 2.73) 1.66] 0.58} 0.91) 1.77] 0.69] 0.0 


Corrected exchange re- 
pulsions 14.87 | 14.56] 13.89 


Corrected exchange 
plus electrostatic 19.17/17.28)15.07/14.63| 14.58 | 13.97 | 13.24 


E(¢) —E(7) 5.93) 4.04) 1.83) 1.39} 1.34] 0.73] 0.0 


18 Wierl, Ann. d. Physik 13, 453 (1932). 
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maximum between the prominent third and sixth 
maxima, Wierl assumed a carbon-chlorine bond 
distance of 1.82A. Debye" reported a chlorine- 
chlorine distance of 4.4A from x-ray photographs. 


DISCUSSION 


There are several contributions to the potential 
barrier to internal rotation, the more important 
ones being the exchange steric repulsion between 
all atoms not bonded to each other, the electro- 
static interaction of the bond dipoles, and the 
interaction between. the bond orbitals on the 
carbon atoms. The magnitude of this last effect 
in H2O2 has been calculated by Penney and 
Sutherland." They find that the potential energy 
of the molecule due to this interaction has 
maxima of about 2 kcal./mole in the cis- and 
trans-positions and minima when ¢ is +90°. 
This effect in ethane or ethylene chloride would 
lead to three approximately equal potential 
minima at angles of 120°. The depth of the 
minima should be less for several reasons. The 
carbon atoms are separated by a greater distance 
than the oxygen atoms. The orbitals on the 
carbon atoms make a larger angle (109° 28’) with 
the carbon-carbon bond than the oxygen orbitals 
make with the oxygen-oxygen bond (90°). The 
potential humps are separated by only 120° 
instead of 180° and the nine repulsions will tend 
to cancel. To show how this might occur we 
have calculated the exchange repulsions'’ be- 
tween hydrogen atoms for the hypothetical 
molecule HxOOH:. The groups are as- 
sumed to be planar and the H—O—O angles 
tetrahedral. The hydrogen atoms are located 
approximately on the loops of the oxygen 
orbitals, the interaction between which has been 
calculated by Penney and Sutherland. Taking 
0.94A and 1.32A for the hydrogen-oxygen and 
oxygen-oxygen distances, the steric repulsions 
between hydrogen atoms is greater by 1.9 kcal./ 
mole when ¢=0° than for the configuration with 
¢=90°. Comparing this with the steric repulsion 
potential barrier of 0.3 kcal./mole, obtained by 
Eyring'® for ethane, it is seen that the nine 
interactions tend to cancel out. If an analogy 
exists between orbital interactions and steric 


4 Debye, Physik. Zeits. 31, 142 (1930). 


16 Penney and Sutherland, J. Chem. Phys. 2, 492 (1934). 
16 Eyring, J. Am. Chem. Soc. 54, 3191 (1932). 
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repulsions, which does not seem unreasonable, 
the orbital interaction in ethane, ethylene 
chloride, etc., would be expected to give rise to a 
small potential barrier, 0.3 or 0.4 kcal./mole. 

The dipole interaction of the two CH,Cl 
groups has been calculated by Smyth, Dornte 
and Wilson.!7 They used the formula for the 
interaction of dipoles at large distances, placing 
the dipole at a point along the carbon-chlorine 
bond. Changing the position of the dipole on the 
carbon-chlorine bond altered the shape of the 
potential curve but did not change the difference 
in energy between the cis- and trans-configura- 
tions. We have repeated the calculations as- 
suming that the dipoles are due to charges 
located at the center of the atoms. This is true 
insofar as the charge distributions on the atoms 
are spherically symmetric. The magnitudes of 
the charges are found by dividing the bond 
moments!® by the interatomic distances. Sum- 
ming the nine electrostatic interactions gives a 
slightly higher potential barrier, 2.4 kcal./mole. 
This potential curve is tabulated in Table II 
and is shown in Fig. 3, curve A. 

The van der Waals attractions have not been 
calculated as the present methods of calculating 
them yield abnormally large interaction at small 
interatomic distances. 

The largest interactions between the two 
CH,CI groups are the exchange repulsions. These 
can be calculated by the method of Eyring.'® 
We have calculated and summed the nine re- 
pulsions as a function of the angle ¢. The Morse 
curve was assumed to be 20 percent coulombic. 
The sum is tabulated in Table II and is shown 
in Fig. 3, curve B. The exchange plus the electro- 
static interactions are tabulatéd in Table II and 
shown in Fig. 4, curve A. 

To see if this potential curve is compatible 
with the electron diffraction data we have cal- 
culated the theoretical intensity formula for it. 
A segment of it is shown in Fig. 2, curve G. It is 
unsatisfactory in the region shown. It does not 
show the fourth and fifth maxima. If the potential 
curve A, Fig. 4, were correct the radial distri- 
bution function would show a peak at 3.2A 
roughly equal in height to the one at 4.29A. 


"Smyth, Dornte and Wilson, J. Am. Chem. Soc. 53, 
4242 (1931). 


“Smyth, J. Phys. Chem. 41, 209 (1937). 
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Fic. 3. Curve A, electrostatic interaction. Curve B, 
exchange repulsion. Curve C, corrected exchange re- 
pulsion. 


That this peak is not present is easily seen in 
Fig. 1. If the repulsions between the chlorine 
atoms were greater at the cis-position by about 
2.5 kcal./mole, the theoretical potential curve 
would be brought into agreement with the elec- 
tron diffraction data. 

Professor Eyring has suggested to us that 
because of the large number of nonbonding 
electrons in Cl, a deeper Morse curve than the 
actual one should be used in calculating the 
exchange repulsions. We have adopted the fol- 
lowing procedure for doing this. The Morse curve 
for Cl, is multiplied by 14/8, the ratio of the 
total number of valence electrons to the number 
of electrons on the valence shell of one atom. 
This will not affect the Morse curve of a bond 
involving a hydrogen atom. Assuming that 20 
percent of this Morse curve is coulombic and 
calculating the chlorine-chlorine repulsions leads 
to the corrected repulsions. When these repul- 
sions are used the potential barrier becomes 
higher. It is tabulated in Table II and plotted 
in Fig. 4, curve B. The electron scattering curve 
calculated for this potential barrier is identical 
with the one calculated for Vo=6.0 kcal./mole, 
curve E, Fig. 2. The agreement between theo- 
retical and experimental potential barriers 
requires that the chlorine-chlorine repulsions be 
calculated by this or some analogous method. 
The same results can be reached by assuming 
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Fic. 4. Curve A, exchange repulsion plus electrostatic 
interaction. Curve B, corrected exchange repulsion plus 
electrostatic interaction. 


that the coulombic part of the Cle Morse curve 
is different (smaller) from the 20 percent used 
for the other calculations of exchange repulsions. 
The procedure above seems to us more justifiable. 

The strongest justification for the above pro- 
cedure is that it leads to results in agreement with 
the electron diffraction data on ethylene chloride. 
This calculation does not change the values 
previously calculated for the hydrogen-hydrogen 
repulsion in ethane. Our estimate of the orbital 
interaction when added to the hydrogen- 
hydrogen repulsions gives a value of about 0.6 


J. Y. BEACH AND K. J. PALMER 


kcal./mole for the potential barrier in ethane, 
The methods we find satisfactory for calculating 
the barrier in ethylene chloride do not lead to 
the high barrier (3 kcal./mole) which has been 
reported for ethane. It is possible that the analogy 
on which our estimate of the orbital interaction 
is based is not sufficiently good. However, it is 
difficult to understand why it should be as bad 
as is indicated by the existing data on ethane. 

Altar” has obtained a potential barrier for 
ethylene chloride from the temperature variation 
of the dipole moment. The theoretical intensity 
curve from Altar’s potential function is identical 
with curve E£, Fig. 2, for which Vo is 6 kcal./mole. 
The shape of his curve near the frans-position is 
the same as the shape of a cosine curve having 
an amplitude of 5 kcal./mole. The difference 
between the two potential curves at small values 
of ¢ cannot be detected in the electron diffraction 
curves because of the small fraction of the 
molecules having small values of ¢. It must be 
emphasized that in putting limits on Vo, the 
height of the barrier, we have assumed a certain 
shape for the potential curve. We have really 
only determined the shape of the potential curve 
near the trans-position of the molecule. At smal! 
values of ¢ the shape of the barrier may be quite 
different from that of a cosine function. If the 
actual barrier is steeper than the cosine function 
at small values of ¢, our Vo is small. If the actua! 
barrier is flatter than the cosine function at 
small values of ¢, our Vo is too high. In any 
event, the barrier is high at small angles. 

We wish to thank Professor Henry Eyring. 
Professor Linus Pauling and Professor Charles P. 
Smyth for many valuable discussions on this 
subject. 
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Calculations of the Lower Excited Levels of Benzene 


M. GoeprEert-MAyvER, Department of Physics, The Johns Hopkins University, Baltimore, Maryland 


AND 


A. L. SkLAR,* Department of Physics, The Catholic University of America, Washington, D. C. 
(Received August 3, 1938) 


The energy of the first excited levels of benzene is calculated by the method of antisym- 
metrized molecular orbitals. The results predict two weak bands, due to forbidden electronic 
transitions, at \= 2500 and \=2100 and a strong band at \=1500. No empirical data except 
the carbon-carbon distance in benzene were used. 


N a previous paper,! one of the authors calcu- 
lated the excited electronic levels of some 
organic molecules, among them benzene, on the 
basis of the Heitler-London method of approxi- 
mation. The same calculation is undertaken here 
for benzene with the method of antisymmetrized 
products of molecular orbitals.’ 


I. THE ORBITALS 


Of the four valence electrons of each carbon 
atom those in the 2s-eigenfunctions and in the 
2p-eigenfunctions whose node is normal to the 
plane of the benzene molecule give rise to sp 
hybridized single bonds to neighboring carbon 
atoms and to the hydrogens. The angles between 
the three sp? hybridized wave functions deter- 
mine the geometry of the benzene molecule. 
The remaining 6 electrons, one for each carbon, 
are in the 2pm atomic eigenfunctions whose node 
coincides with the plane of the benzene ring. 
These electrons are responsible for the aromatic 
properties of benzene and for the lower absorp- 
tion bands.! We are here concerned with the 
different molecular states of these 6 electrons 
only. 

The electrons shall be denoted by Greek 
letters v, uw, 1, --+, 6. The different carbon atoms 
will be denoted by k. The 2p7 eigenfunction of 
the kth carbon atom will be written K(v); 
specifically I(v), +--+, VI(v). They are all taken 


*This work was supported by a grant from the Penrose 
Fund of the American Philosophical Society. 
vA L. Sklar, J. Chem. Phys. 5, 669 (1937). 
R. S. Mulliken, J. Chem. Phys. 1, 492 (1933); Phys. 
Rev, 32, 186, 761 (1928); 33, 730 (1929); 40, 55, 41, 49, 
151; 43, 279 (1932-3). J. E. Lennard-Jones, Trans. 


Faraday Soc. 2 i i 
37, 601 (1999) 668 (1929). G. Herzberg, Zeits. f. Physik 


to be positive on the same side of the benzene 
plane. 

The interaction between the 6 electrons is at 
first neglected. The potential 7) acting on the 
electrons has then the symmetry of the benzene 
molecule and may be written as a sum of the 
contributions from each atom. 


Ho= Dx. (1) 
K 


Of the 6 K-functions the following 6 mutually 
orthogonal orbitals have to be formed? 


6 
k=1 


with /=0, +1, +2, 3. The coefficients o;~! take 
care of the normalization. If the overlapping 
integral fJ(v)II(v)dr, is neglected, all o’s be- 
come unity. 

With this approximation, and neglecting of 
integrals over products of K’s which are not 
neighbors, the energies of the orbitals are: 


f I(T +H,)Idr 
+2 cos (271/6) f I(T+H)IIdr, (3) 


where T is the kinetic energy. In section IV the 
energies will be evaluated more exactly. 

It is seen that the energies of the orbitals 
increase with |/| since /J(7+1°)IIdr is nega- 
tive. The energy depends only on the absolute 
value |/|, showing that the levels with |/| =1 
and |/| =2 are doubly degenerate. 


3 E. Hiickel, Zeits. f. Physik 70, 204 (1931). 
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The energy of the molecule is lowest if the 
orbitals $0, ¢1, ¢-1 contain two electrons each. 
The eigenfunction of the ground state in this 
order 01 approximation is the product 


(4) 


The first excited level arises if one electron 
from the state /=+1 or /=—1 is lifted to /=2 
or 1= —2. Its energy above the ground state is 
é2—€;. To this level belong four different eigen- 
functions: 


¥1 = 

b-2(6), 

(5) 
¥4= 


II. Group THEORETICAL CONSIDERATIONS 


The way in which this 4-fold level of the 
molecule as a whole splits up, due to electronic 
interaction, and the behavior of the resultant 
molecular states can be determined without 
calculation by the application of group theory. 
The potential Hy has the symmetry of the 
benzene ring, which is transformed into itself by 
rotations and reflections forming the group Dgn. 
The group operations transform any orbital ¢, 
into a linear combination of orbitals of the same 
energy, with coefficients which are the matrix 
elements of an irreducible representation of the 
group. The same holds true for a state of the 
molecule as a whole. The interaction between the 
electrons cannot alter this; it is only capable of 
splitting up ‘‘accidental” degeneracies, namely 
to separate the energies of states belonging to 
different irreducible representations which hap- 
pened to be energetically equal in the zeroth 
approximation. 

We wish therefore to investigate which irre- 
ducible representations of the benzene group are 
contained in the reducible fourth order repre- 
sentation determined by the transformation 
properties of the functions (5). 

From the group tablet of De, by the usual 
method of comparison of the characters, the 


‘ Bright Wilson, Phys. Rev. 45, 706 (1934). 


irreducible representations of the one electron 
orbitals can be found immediately. They are 


1, 6(Eig); 
T4(Bag), 


where the symbols in brackets are the notation of 
Mulliken and Placzek.> Every closed shell trans- 
forms like the identical representation T;(A,,), 
that is, remains unchanged under the operations 
of the group. We find therefore for the eigen- 
function of the ground state 


WoC 


The four eigenfunctions of the excited level 
transform like the products namely: 


XPu =(To + Piz) 
Eig X 


The excited level will therefore split up into two 
single states belonging to I'y and Tio, and one 
doubly degenerate one belonging to the two- 
rowed representation 

Transitions between an excited state and the 
ground state with emission or absorption of 
light occur only if the integral of at least one 
component of the electric moment P multiplied 
by the eigenfunctions of ground- and excited 
state does not vanish. This is the case if the 
integrand is completely symmetrical, or at least 
contains an additive part which has this prop- 
erty. In other words, an excited level will be 
active in combination with the ground state 
only if the direct product of its representation 
and the representation of P contains the identical 
representation I',. The transformation properties 
of P are the following: 


P,, (normal to benzene plane) ¢ I's(Ao.) ; 
P,, P, (in the plane) ¢ 


None of the products of I's with T's, Tio or T'n 
contains T,; nor do the products of Tj. with 
Ty and Tyo. Only T'wXT2 contains a totally 
symmetrical part, so that T'12 is the only active 
level, and is capable of emitting light polarized 
in the benzene plane, that is, perpendicular 
bands. 


5 Mulliken, Phys. Rev. 43, 279 (1933). 
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Transitions from the states I'y and Tyo to the 
ground state are, however, possible with simul- 
taneous emission or absorption of vibrational 
energy, for instance of one quantum of vibration 
frequency of symmetry I;(f2,) that is of a 
degenerate, Raman active vibration (frequencies 
vg, 7, Ys, ¥9 in the notation of E. Bright Wilson) ; 
or in combination with emission or absorption 
cf two quanta cf any degenerate vibration. 
It is to be expected that these transitions are 
weak, since they would vanish completely if the 
eigenfunction of the total molecule could be 
expressed as a product of vibrational and elec- 
tronic functions. 

The linear combinations of the four product 
eigenfunctions y (5) which have the symmetries 
Ts, Tig respectively are: 


Bou) ¥ ys), 


V3=y3, 
T 
= 


III. ELEctRonNic INTERACTION 


The change of energy of the states due to the 
electronic interaction 


v=1 (7) 


will now be calculated. 

The molecular eigenfunctions of the ground 
state (4) and the excited states (6) must be 
multiplied by an appropriate spin function and 
made completely antisymmetrical with respect 
to the exchange of electrons. The two possible 
spin functions for the vth electron will be denoted 
by a(v), B(v). If two electrons are in the same 
orbital, they must necessarily have opposite 
spins. For the ground state, therefore, the spin 
function 


Xo = (8) 


gives the most general result. 

For the excited states two cases must be dis- 
tinguished. The spins of electrons 5 and 6, the 
only ones which are unpaired, may be either 
Opposite or parallel. In the first case, the spin 


function must be made antisymmetric in 5 and 
6, namely 


Xs = 
X [a(5)8(6) —B(5)a(6) J. (9) 


In the second case it must be symmetric in 5 and 
6, which can be done in three different ways, 
corresponding to the three orientations of the 
sum of the spins in space, one of which is 


(10) 
while the others contain 8(5)8(6) and 
2-!La(5)B(6) +8(5)a(6) J. 


Since the perturbing potential //' (7) is spin 
independent, spin states (9) and (10) do not 
interact. We shall refer to the states described 
by eigenfunctions obtained with (9) as singlet 
states, to those obtained with the three spin 
functions symmetric in 5 and 6, of which (10) 
is one, as triplet states. The energy of the latter 
three is the same in this approximation. The 
magnetic interaction between spin and orbit may 
however, bring about a fine structure splitting. 

The totally antisymmetric eigenfunctions are 
then obtained by performing all 6! permutations 
P of electrons on the product of spin and co- 
ordinate function, multiplying the result by —1 
to the order (odd or even) of the permutation 


and adding. For instance, for the ground state 


(6 1)? PWoxo. 


Since the molecular orbitals ¢; are. siricily 
orthogonal, the average value of H' for any such 
function can be expressed without any approxi- 
mation, since all integrals over products of 
functions in which more than two electrons are 
in different molecular orbitals are zero. The 
energy will contain Coulomb terms, of the form 


e 


and exchange terms 


e 
f f (12) 


Of the product eigenfunctions (4), ~3 and y,4 
interact with none of the permutations of any 
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other one, bearing out the fact that the eigen- 
functions (6) are the correct ones. The inter- 
action between y; and ye brings in two more 
integrals, 


f (13) 


e 


All integrals y, 6, £, 7 are real. 
The energies, including electronic interaction, 
are then for the ground state 


Eo=2€0+4 261-1, (15) 


where €o and ¢; are the energies of the orbitals (3). 
For the excited singlet states, the energies are 


'Pio(Beu), (16) 
E3,=E® —6_:2+ 612, 
and for the triplet states 
Exn=E™ (17) 
—6_12— 412, iw), 


where the term E! which all excited levels have 
in common is 


6yoit27o2 
(18) 


The subscripts 1, 2, 3 correspond to those on the 
eigenfunctions (6). 


IV. THE ENERGIES IN TERMS OF ATOMIC 
INTEGRALS 


The next task is to express all the integrals 
y, 6, & » over molecular orbitals in terms of 
integrals over the atomic eigenfunctions K(v) 
through the use of (2). 

The normalization coefficients o; of the orbitals 
are first determined by the relation: 
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+2 cos 


=1+2 cos (2nl/6) To) 


+2 cos (4nl/6) f 
+(=1)'f 


The functions K(v) are assumed to be normalized. 
If the “‘overlap integrals’ between neighbors are 
denoted by si, between next to neighbors se and 
between opposite places in the ring s3, the above 
expression can be written 


4rl 
o,;=1+2s; cos cos ae 1)'s3. (19) 


The integrals y,. and 6, contain the terms 
in the denominator, the integrals and 7 
the term o,:o2. The expansions of the orbital 
integrals contain molecular integrals of Coulomb 
type between electrons on the same and on 
different atoms, 


2 
f (20) 


Tou 


secondly a mixed exchange and Coulomb integral 
the interaction between an electron on one 
atom with one that is shared between the same 
and the neighboring atom, 


B= (21) 


and finally the exchange integral between elec- 
trons on neighboring atoms 


Ca f f (22) 


Interactions between electrons which are pat- 
tially on atoms which are not nearest neighbors 
shall be neglected, except for the purely Coulomb 
terms (20). All integrals A, B, C, s, are positive. 
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In this nomenclature, formulae( 16) and (17) for the energies of the excited states take the form 


E,,=E™ + (1/6012) { 


E3,=E™ + (1/6010) { 


+(1/60102){ —3A94+2A1 
Ex —A 
—2Ao+ Ai- Aot2A3—3C}, 


The order of the singlet terms will be shown 
to be E:,<FEi1,<E3,, that of the triplet states 
Ey. <E3:< Et, with the triplet states lower than 
the singlets. 

Thecommon term E™ (18) of the excited states, 
as well as the energy of the ground state Eo (15) 
are too complicated, due to the different normali- 
zation factors, to be written down in the general 
form. We are, here, however, interested only in 
the height of the excited levels above the ground 
state, and therefore in E“—Ep,. Using the 
numerical values of 51, Sz, 53, and therefore of the 
normalization coefficients oo, 71, 72 (see Table I), 
we obtain 


—Ey= €2— €:— 2701+ 
+3712+ 601 — 61-1 
= €:+0.4,A o+0.8,A 1+0.64A » 


(25) 


+0.23A,—1.15B—1.26C. 


Lastly, we are confronted with the calculation 
of the energies «, of the orbitals due to the 
potential of the benzene rings. An approximate 
form for these was given in (3). The complete 


TABLE I. 


SI(v)II(v)dr, 
fI(v) ITI (v)dr, 
Ao= J 
IP(u)drydr 
f 
A;= SSE /ryT*(v)I V*(u)dr dry 
IT (v)I(u) LT (u)drydry 
— SH i(v)IP(v)dr, 
— JH i(v)I(v)IT(v)dry 


02 
16.82 volts 
8.78 


I(v), II(v) designate 2x carbon functions (32) around nuclei J, JJ, 
respectively, H is the potential of a neutral carbon atom (33). 


Ae—6A;4+8As—3As+ C}, 
Eo, = E™ +(1/60102) { —Ay—2A2 


0 —2A2+3A;3—7C}, 


'To(Biu), 
'Ti0(Beu), 
'Tie(E,u), 

*T's(Biw), 
*Ti0( Bou), 


+3A3s+5C}, (23) 
~ 


+ A3t+ 


3A1—3A2 


expression is 


1 
f I(v)(T-+Hp){ I(v) +2 cos (2x1 /6)IT(v) 


+2 cos (26) 


The effect of the hydrogen atoms will be neglected 
completely. In the expression 


Ho= 
K 


the potential Hx is then simply the attraction of 
the carbon nucleus and the repulsion of the other 
five carbon electrons. Since J(v) is supposed to be 
an atomic eigenfunction, the relation 


(T+H1(v))I(v) = Weyl (v) (27) 


holds, where W2, is the energy of a 29 electron in 
a carbon atom in the valence state. This term 
occurs once in every one of the 4 terms of €; (26), 
with coefficients that add up to a, precisely 
canceling the o; in the denominator. W3, there- 
fore drops out of the difference €2:—€; since its 
coefficient in each ¢; is independent of /. 

In the potential Hx the two electrons of the 
inner shell of the carbon atom shall be treated as 
if they were located in the nucleus. Hx(v) is then 
the electrostatic potential at the place v of a 
nucleus of charge 4 plus one electron in each of 
the three 2spo hybridized wave functions, or, 
which is the same, of the nucleus and one electron 
each in a 2s and the two 20 eigenfunctions. The 
calculation of the new integrals can be simplified 
by noticing that Hx is the potential of a neutral 
carbon atom with one electron in the 2s and each 
of the three 2 functions, minus the potential of 
a 2pm electron. The potential of the neutral 
carbon atom, Hx is a spherically symmetrical 
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attraction and we may write 


Hex(v)=Hx(v)— (Ku). (28) 


If (27) and (28) are inserted into (26) the second 
part gives rise to integrals of the types A and B 
(20), (21), the first part to new integrals, namely 
the effect of one neutral carbon atom on an 
electron at another nucleus, 


Q= (29) 


and the effect of the potential on an electron 
partially on the same and on the neighboring 
atom 


The minus sign is introduced in order to make 
both Q and R positive. All terms involving 
Hx which arise from interactions of non-neighbors 
are small enough to be neglected. The difference 
of and is then 


01-02 


( As} 
oitos 
= —1.12042.22R—1.12A; 
—1.1242—0.56A3+2.22B. 


+ {R+B} (31) 


If (31) is put into (25) we obtain the final result 
for the part of the energy common to all excited 
states above the ground state 


—Ey= —1120+2.22R+0.41A0 
—0.31A:—0.4842—0.33A; 
+1.07B—1.26C. (32) 


V. NUMERICAL VALUES 


The eigenfunctions used for the calculation of 
the numerical integrals are 27 hydrogen like 
functions at the different nuclei, 


5 \ 3 
K(r, 0,?) = rexp (—2r/2)sin@cos¢. (33) 


M. GOEPPERT-MAYER AND A. L. SKLAR 


The screening constant z was chosen to be 
z=3.18, a value which was obtained by Zener® on 
the basis of a variational calculation minimizing 
the energy of the ground state of one carbon 
atom. This screening constant is not necessarily 
the most advantageous one for our molecular 
calculation. Minimizing the energy of one carbon 
atom in its valence state would presumably give a 
more correct value for z. 

The distance between carbon atoms is 7 
= 1.39A,7 which leads to 279 = 8.37 in atomic units. 

The numerical values of the integrals used are 
collected in Table I. 

The overlap integrals s can be calculated in 
closed form.’ Ao, the repulsion between two 
electrons on the same atom, was calculated for 
this paper. The Coulomb repulsions between 
electrons on different atoms, A, Ao, A3, as well as 
the exchange integral C were taken from the 
tables of Bartlett.’ 

The integral B, the repulsion between one 
electron on one atom and one shared between the 
same and a neighboring atom was not to be found 
in the literature. The evaluation was complicated 
and details of the calculation will be published 
shortly by one of the authors and R. H. Lyddane. 

The attraction of a neutral carbon atom, in the 
valence state for an electron is: 


4e? 1 Ty 
H,;= e’a(r)rdr, (34) 
ly Ty 0 Ty 


where 


o(r) =4r{ | 
+ | (35) 
r, signifies the distance of the electron from the 


nucleus of the carbon atom I. 
Since 


f a(r)r'dr=4 


(34) may be simplified to 


H,;= (36 
Jry Vy 


6 Zener, Phys. Rev. 36, 51 (1930). 

7 Handbuch der Physik, Vol. 24, No. 2, p. 143. 

8 Handbuch der Physik, Vol. 24, No. 1, p. 643. ' 

9 J. H. Bartlett, Phys. Rev. 37, 507 (1931). In Bartlett = 
notation, A =7;+ 477 and C=12+ 
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The effects of this potential of carbon atom J on 
an electron on the neighboring nucleus J/, Q, as 
well as on an electron shared between atoms J 
and JJ, R, were calculated. Since the potential 
H, decreases rapidly with distance, its effect on 
other than neighboring electrons was found to be 
completely negligible. 


VI. RESULTS 


The levels calculated in this manner are shown 
in Fig. A. The numbers at the right of the lines 
give the energy in volts of the levels above the 
ground state. Since only the excitation of one 
electron from /=1 to 1=2 has been considered, 
only the lowest levels are obtained. The next 
higher ones, however, arising from the excitation 
of one electron from /=1 to /=3, one electron 
from /=0 to /=2 or two electrons from /= +1 to 
l= +2 have different symmetry properties from 
any of those obtained here (for instance I) 
and therefore do not interact with the states 
considered. 

Figures B, and Bz show for comparison the 
level scheme obtained from the Heitler-London 
method of calculation.! Figure B, contains the 
levels as calculated from nonpolar terms, Be 
those obtained by the inclusion of some polar 
configurations. The contribution of the polar 
terms is somewhat uncertain since not all 
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possible polar terms, but only the lowest cnes, 
were considered, and since the energy depends on 
two parameters which can only be roughly 
determined empirically.” 

The energy of the lowest singlet state, of 
symmetry 'IT'jo agrees very well in the two 
methods. Experimentally, a band of oscillator 
strength f=10- is observed at }\=2600A or 
energy about 4.7 volts. The small transition 
probability indicates that this band is due to a 
forbidden electronic transition. Analysis of the 
band structure shows that the observed vibra- 
tional progressions are in agreement with a 
symmetry character 'I'yo(Bo.) for the excited 
electronic state. 

Of the other singlet states, the energies of 
'T,o(E,.) agree well in both calculations. A 
glaring discrepancy occurs however in the posi- 
tion of the level 'IT'9. In the Heitler-London 
calculation both 'T, and 'Ty2 are purely polar 
terms. They are obtained in that method really 
as a second approximation, and it would not be 
surprising if the calculated energy would turn 
out to be too high. On the other hand, in the 
method of this paper, polar configurations enter 
with too large coefficients into all states, causing 
the polar states to be relatively too low, as will be 
discussed later. But the agreement for the level 
must then be regarded as purely 
accidental. 

The experimental evidence in the region of the 
far ultraviolet is not completely conclusive. The 
measurements of Henri! in solution, evaluated 
by Chako” show a band of oscillator strength 
one-tenth at about 2070A, or 6.0 ev. One must be 
inclined to attribute this to the electronic level 
11,2, since this is an allowed transition. On the 
other hand, the measurements of Carr and 


10 The two parameters are the difference between the 
Coulomb energy of a polar and a nonpolar state and the 
resonance integral fwa(1)y¥o(2)Hya*(1)Wa*(2). Their values 
were estimated from the absorption frequency of pentene 2 
(since cyclohexene has not been measured), and from the 


+ 
resonance of cyclohexene. States like 


were omitted since they introduced new variables and 
complicated the determinant. 

VY. Henri, J. de. Phys. 3, 181 (1922). 

12 N. Chako, J. Chem. Phys. 2, 644 (1934). 
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Stiicklen® in the gas seem to indicate the 
presence of two bands at 2040A and 1850A, of 
6.0 ev and 6.7 ev respectively, the band at longer 
wave-length being the weaker one. This would be 
in reasonable agreement with the calculations of 
this paper, which predict a state forbidden in 
combination with the ground state, therefore a 
weak band, at 5.8 ev, and an allowed band at 
8.0 ev. Since in solution the 2600A band of 
benzene is considerably shifted to longer wave- 
length, it is conceivable that a strong band at 
1850A might be shifted sufficiently in the liquid 
to mask a weak band at 2040A, which would 
correlate the experiments of Henri with the 
appearance of the bands in vapor. 

The energy of the triplet states turns out to be 
uniformly lower in the molecular orbital than in 
the Heitler-London calculation. The best agree- 
ment is obtained for *I'y, the worst for *I'y9. That 
the energy, calculated on the basis of molecular 
orbitals, is lower for the triplet states than the 
singlet states can be seen qualitatively. Since the 
eigenfunctions used are antisymmetrized they 
allow no more than two electrons of opposite spin 
to be simultaneously at one nucleus, but they 
allow this to happen very frequently. Terms 
corresponding to two electrons on one carbon 
atom give rise to a large repulsion and therefore 
an increase in the energy. Now in the triplet 
states, where four electrons have spin a, two 
spin 8, this piling up of two electrons of opposite 
spin will not occur as often as in the singlet 
states, causing the calculated energy of the 
triplet states to be lower than that of the 
singlets. For the true electronic eigenfunctions 
two electrons will not accumulate as freely on one 
carbon atom, so that one must expect the true 
difference between the energies of the excited 
singlet and triplet states to be smaller than that 
calculated on the basis of molecular orbitals. In 


1%E, P. Carr and’ H. Stiicklen, J. Chem. Phys. 6, 55 
(1938). 
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other words, the absolute value of the energies of 
all states computed on the basis of molecular 
orbitals are too high, but those of the singlet 
states more too high than those of the triplet 
states. Since we calculate the energies of all 
excited states above the ground state, which is a 
singlet state, it may be expected that the energies 
of the singlet states are more reliable than those 
of the triplet states. 


VII. SUMMARY 


The eigenfunctions used for the description of 
the states of the benzene molecule are built from 
products of molecular orbitals. They are made 
completely antisymmetric. The orbitals in turn 
are assumed to be linear combinations of the 27 
atomic eigenfunctions of the six carbon atoms. 
No assumptions other than these were made. 
The energy is obtained as a sum of integrals over 
atomic eigenfunctions. Except for the Coulomb 
term, interactions between electrons farther apart 
than on neighboring atoms were neglected. This 
approximation is not necessary, since the other 
occurring integrals can be estimated very easily. 
There exist, however, a great number of them. 
and the calculation of the coefficients with which 
they enter into the energy of the various levels 
becomes extremely tedious. We took into account 
the effect of quite a few of them and found that 
they tend to cancel. No other approximations 
were made and no empirical data were used. 
except the C-C distance in benzene. 

The energies depend then on the form of the 
atomic eigenfunctions. The integrals which are 
most sensitive to a change in the screening 
constant are the overlapping integrals s, and they 
also enter rather critically into the energy. The 
screening constant used here was 3.18. 

The authors wish to express their thanks to 
Professors K. F. Herzfeld and J. E. Mayer for 
their interest and for many valuable discussions. 
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LETTERS TO THE EDITOR 


not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per page) 
will not be made and no reprints will be furnished free. 


This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
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this section must reach the office of the Managing Editor 


The Absorption Coefficients of the Acetyl Halides 


The absorption coefficients of acetyl chloride, bromide 
and iodide for some of the lines of the mercury arc spec- 
trum have been measured with a thermopile and mono- 
chromator. The values obtained are listed in Table I. 
The long wave-length limits for the absorption, as esti- 
mated by an extrapolation of these data, are 2900, 3100, 
and 4000A for the chloride, bromide, and iodide, respec- 
tively. The absolute values of the absorption coefficients 
and these limits may be in error a few percent because of 
the reactivity of the compounds and, in the case of the 
chloride and iodide, the separation of condensed phases as 
the result of the reaction occurring. 

Preliminary analyses of the products formed in the de- 
composition of these compounds caused by the absorption 
of ultraviolet light indicate that the reactions involved are 
quite complex. No hydrogen, ethane, or ethylene was 
detected by suitable tests for each.! The non-condensable 
(in liquid air) gases consisted of carbon monoxide and 
methane. The non-condensable portion obtained from 
acetyl chloride on illumination with A=2537A was 60 
percent carbon monoxide and that obtained from the 
bromide or iodide was approximately 90 percent carbon 
monoxide. The condensable products contained the 
corresponding methyl halide but otherwise showed marked 
variations. 

In the case of acetyl chloride a liquid product separated 
out during the course of the illumination. This substance 
appeared first as a white mist which became yellow as it 
settled on the walls and finally it changed to a dark brown. 
The presence of some hydrogen chloride and unchanged 
acetyl chloride interfered with the analysis of the liquid 
but the carbon and hydrogen percentages were found to be 
in fair agreement with the hypothesis that this substance 
is a polymer of either diacetyl or ketene. Tests of water 


TABLE I. Absorplion coefficients. 


solutions of the condensable gases for carbonate and 
oxalate (from (COCI:) ions) were negative. 

Acetyl bromide did not yield a liquid such as was ob- 
tained from the chloride. In addition to methyl bromide 
the condensable products contained some bromine and 
compounds containing more than one atom of bromine per 
molecule, probably ethylene bromide. The products from 
acetyl iodide contained considerable quantities of free 
iodine and some liquid which was not analyzed on account 
of the small amount obtained. 

Detailed studies of the decomposition of these three 
compounds are in progress. 

D. H. EtTz_er 


G. K. ROLLEFSON 
Depart ment of Chemistry, 
University of California, 
Berkeley, California, 
August 25, 1938. 


1 Blacet and Leighton, Ind. Eng. Chem. Anal. Ed. 3, 266 (1931); 
Blacet, MacDonald and Leighton, ibid. 5, 272 (1933); Blacet and 
MacDonald, ibid. 6, 334 (1934); Blacet and Volman, ibid. 9, 44 (1937). 


Raman Frequencies of Hexamethylene Glycol 


In the course of preparation of a series of glycols, Mr. 
George W. Anderson of the Department of Chemistry of 
this University supplied the writer with a specimen of 
hexamethylene glycol (HO-(CHz2)s-OH), with a melting 
point of 40.9 deg. C, and a boiling point range of 135 5 
to 136.5 deg. C at a pressure of 8 mm Hg. 

In determining the Raman frequencies for this com- 
pound, mercury arcs in pyrex were used for excitation. 
In some runs a filter consisting of a saturated solution of 
sodium nitrite was used to eliminate the Raman lines 
resulting from the mercury 4047A group of lines. An iron 
spark spectrum adjacent to the Raman spectrum was used 
for the purpose of interpolating the wave numbers of the 
Raman lines. 

The values are as follows, with estimated relative in- 
tensities and mercury line sources according to the Kohl- 


Wave-l 
2537 216 “as — 814 (2-e), 852 (2-e), 874 (2-e), 915 (1-e), 1000 (2-e), 
2652 10.3 46.8 324 1052 (3-e), 1084 (3-e), 1298 (4-ek), 1435 (4-ekc), 1472 (4-ekc), 
2783 Ae 128 216 2855 (8-ek), 2905 diffuse (8-ek). 
2.9 5.0 R. C. WiLtiaMson 


Department of Physics, 
University of Florida, 
Gainesville, Florida, 

September 15, 1938. 
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Structure of Liquid Aliphatic Amines 


The Table I below gives the results of x-ray diffraction 
study in five liquid di-acyl amines and in n-amyl and 
tri-n-amyl amines. The usual technique of production of 
x-rays and of obtaining diffraction haloes in liquids was 
employed. The second column gives the Bragg spacing d, 
calculated from the diameter of the diffraction halo. The 
third column gives the constant k, in Raman! and Ram- 
nathan’s theory, i.e., the ratio of d to (m/p)!, which for 
spherical molecules should have a value between 0.8 and 
unity. 

It is seen that the value of Bragg spacing has the con- 
stant value 4.46A for all the normal di-acyl amines ex- 
amined. This value is close to 4.60A, which is the diameter 
of a carbon-atom-chain in aliphatic liquid molecules found 
by previous workers.? This constant value for the spacing 
cannot represent the distance between two carbon-chains 
in the same molecule of an acyl amine, since the diffraction 
haloes are far from sharp. The possible explanation of 
this value of the spacing is, therefore, that it represents the 
distance between two carbon-chains belonging to two 
adjacent molecules of the amine. 

We may assume that the flat pyramidal shape of the 
NH; molecule is maintained in the acyl amine molecule, in 
which one, two or three carbon-chains replace the hydro- 
gen atoms. Taking into account the tendency of long chains 
to arrange themselves parallel to one another and adopting 
Stewart's? idea of ‘“‘Cybotactic state’’ in liquid, the re- 
quired condition for diffraction is obtained if we suppose 
that in the liquid, the molecules are piled up one above the 
other, the vertices of the adjacent pyramids fitting into 
each other and the corresponding edges lying close and 
parallel. Here, the Bragg spacing will be the mean distance 
between two carbon-chains lying along corresponding 
edges of adjacent pyramids and this will be nearly equal to 
the thickness of the chain, since the packing can be as close 
as possible. The whole arrangement of the molecules will, 
therefore, somewhat resemble a row of rather flat tripods, 
placed one above the other, the three legs of one lying close 
and alongside those of its neighbor. It is found that in the 
haloes of the tri-acyl amines, there is considerable intensity 
within the diffraction halo but the outer edge is sharply 
defined. This is readily explained on the above picture, as 
the outer limit corresponds to the closest approach of one 
edge of a pyramid towards the corresponding one of its 
neighbor. Any disarrangement of the row of pyramids will 
increase the mean distance between corresponding carbon- 
chains and give the inner scattering of the halo. 

It will be also seen from the table that the spacing for an 


TABLE I. 


Liquid é= 2 sin 6/2 


Di-n-propyl amine 
Di-iso-propyl amine 
Di-n-butyl amine 
Di-n-amyl amine 
Di-iso-amyl amine 


n-amyl amine 
Di-n-amyl amine 
Tri-n-amyl amine 


iso-amine is larger than that for the normal one—an effect 
of the side linkage in the chain. Also, the value of the 
constant & in the last column is, in all cases but one, 
below 0.7, showing that the molecules are not spherical 
in shape. Only in the case of di-iso-propyl amine, it is more 
than .8. This must be due to the fact that this molecule 
is more nearly spherical in shape, since there is a side 
linkage attached to a short carbon-chain. 
B. V. THosar 


Physics Department, 
College of Science, 
Nagpur (C.P.), India, 
August 16, 1938. 


1 Raman and Ramnathan, Proc. Ind. Ass. Cult. Sci. p. 127 (1923). 
2 Stewart and Morrow, Phys. Rev. 30, 232 (1927). 
3 Stewart, Rev. Mod. Phys. 2, 116 (1930). 


Mayer’s Theory of Condensing Systems 


Mayer’s general results! can be derived almost im- 
mediately from the theory of dissociating assemblies, thus 
avoiding much of his preliminary argument. Fowler’s gen- 
eral method? for gaseous assemblies with any number of 
components and reactions is at once applicable to a gas of 
one component and any number of reactions—clustering 
processes or collisions. The dissociating molecule of type 
S is now a collision or cluster of S atoms. Fowler's equations 
(479) and (484), giving the law of mass action and the 
external pressure on the gas, are then, with classical 
partition functions, almost obviously identical with 
Mayer’s equations (43) and (37) when recast into the 
notation there used. The present more general argument is 
applicable to any type of partition function, proving very 
directly the generality of Mayer’s results. 

WILLIAM BanpD 

Department of Physics, 

Yenching University, 


Peking, China, 
August 1, 1938. 


1 Joseph E. Mayer, J. Chem. Phys. 5, 67 (1937). 
2R. H. Fowler, Statistical Mechanics, 2nd Edition, §5 -3. 


On the Symmetries of the Fields About Ions in Solution. 
Equilibrium between Forms of Different Symmetry 


The absorption spectra of solutions of trivalent europium 
are so sharp that it is possible to resolve the multiplets 
induced by the electrical fields about the europium ions. 
It was found! that the multiplets consisted of different 
numbers of components depending upon the nature of the 
solutions. Such differences were sought as a parallel to the 
different structures which arise in the spectra cf crystals 
when a rare earth is subject to crystal fields of various 
symmetries. 

In this note we shall describe new features in the spectra 
of aqueous solutions of europium nitrate and of europium 
chloride, especially in the blue region which contains 
lines rather easy to follow. This group consists of four 
lines in the chloride and three lines in the nitrate of differ- 
ent spacing, sharpness and relative intensities? These 
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4650A Band 


4647 4654 


Fic. 1. 1-Fe arc. 2-EuCls .0025 M. 3-EuCls .95 M. 4-Eu(NOs)s .0025 M. 
5-Eu(NOs3)3 .009 M. 6-Eu(NOs)3 .95 M. 7-Fe arc. 


separate patterns remain unchanged from about 1.5 M to 
about 0.01 M as observed with a 21-ft. grating. On diluting 
to 0.0007 M, we could discover no change in the pattern of 
the chloride but when the nitrate was diluted to 0.0007 M, 
its group which had consisted of three lines now had 
changed over to four, indistinguishable in character from 
those of the chloride. At concentrations between 0.01 M 
and 0.0007 M, the nitrate showed both spectra and with 
increased dilution the chloride pattern gained intensity 
at the expense of the nitrate pattern. (See Fig. 1.) 

Europium nitrate in aqueous solution has then dis- 
tributed itself between two discrete configurations in 
equilibrium corresponding to different symmetries of the 
fields about the europium ions. There is some resemblance 
here to the situation in which two allotropic modifications 
of a crystal are present together, with this difference, that 
the equilibrium in the solution of the nitrate takes place 
within one phase. 

In a forthcoming note it will be shown that europium 
nitrate and europium chloride are strong electrolytes 
similar to and of the same order of strength as lanthanum 
nitrate and lanthanum chloride which have served as 
standard salts in thermodynamic studies of trivalent ions. 
Because of the close resemblance of europium to other 
rare earths, we regard the existence of discrete forms of 
different symmetry as no rare phenomenon—among 
trivalent ions, at any rate. 

SIMON FREED AND H. F. JACOBSON 


H. Jones Laboratory, 
University of Chicago, 
Chicago, Illinois, 
September 2, 1938. 


- Ireed and S. I. Weissman, J. Chem. Phys. 6, 297 (1938). 
there, the publication of the previous note, we have discovered 
t G. Urbain (J. Chem. Phys. 4, 243 (1906)) had observed differences 


tween the spectrum of europium nitrate and of europium chloride in 
aqueous solution. 
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Symmetries of the Fields About Ions in Solution and 
Electrolytic Dissociation 


Recent experiments! have shown that the fields about 
the ions of trivalent europium in aqueous solution are 
sharply defined in their intensities and orientations and 
that these differ with the nature of the solution. However, 
with a solution of a given salt, the structures in the spectra 
which respond to changes in the fields are not affected 
when the solution is diluted. The wave-lengths of the 
lines, their sharpness, relative intensities, the spacings 
within the multiplet, all remain unchanged, from concen- 
trated to dilute solutions. Europium chloride has four 
lines in a particular group in the blue region. No change 
could be observed in the character of the lines when the 
concentration was varied from about 1.5 M to 0.0007 M. 
Corresponding to this group of lines, europium nitrate has 
three lines, altogether different in quality from the four 
of the chloride. These lines experience no change when the 
solution is diluted, even though a new additional form 
makes its appearance? at about 0.01 M. The spectrum of 
the new form is indistinguishable from that of the chloride. 
As dilution proceeds the ‘‘chloride-structure’’ of the ni- 
trate becomes more intense relative to the original ‘‘ni- 
trate-structure.’’ Nevertheless, even though the total 
intensity of each structure varies, the relative intensities 
within the structure, the sharpness, and the intervals 
remain invariable throughout dilution. 

The theory of Bethe® prescribes into how many com- 
ponents a particular energy level will be decomposed by 
an electrical field of given symmetry. The magnitude 
of the intervals between the components is a measure of 
the strength of the fields. We have found that the spec- 
trum of a crystal of hydrated europium nitrate contains 
three lines in its blue group, having roughly the same 
intervals as originate from a solution of the nitrate. We are 
then to conclude that the fields about the europium ions in 
the crystal are roughly of the same intensity as the fields 
about the ions in solution. Furthermore this intensity does 
not change with dilution. 

A natural assumption to make at this point“is that the 
salts of europium are revealing bonds which are perma- 
nently oriented. If that description is to be correct, it must 
not exclude the electrolytic dissociation of the salts since, 
as we shall show, salts of europium are strong electrolytes 
in about the same measure as the salts of lanthanum. 
The latter have repeatedly served to illustrate and test the 
inter-ionic attraction theory. 

There is such a mass of different data on the close re- 
semblance among the ions of the rare earths‘ that there 
can be little doubt concerning the resemblance with 
respect to ionization also. As a final“check on these infer- 
ences, which, to be sure, involve little interpolation, Mr. 
I. M. Klotz of this laboratory measured the electrical con- 
ductivities of solutions of europium nitrate and of euro- 
pium chloride and found them to be nearly of the same 
magnitude as those of the corresponding lanthanum salts. 
Mr. Klotz will publish his data shortly. As to thermody- 
namic behavior—the freezing point lowering of lanthanum 
nitrate at 0.01 M is 3.35 times® that of the same con- 
centration of mannose and at 0.001 M it is 3.75 times. 
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Judging from the spectra we have been led to assign 
rather permanent, discrete geometrical configurations to 
the fields about the ions of europium notwithstanding the 
fact that such configurations have long been imagined 
characteristic of non-electrolytes. This behavior is not 
limited to salts of europium. It is not due to some 
singular electronic distribution which europium possesses,® 
since €uropium ion acts in its salts and solutions like 
the spherically symmetrical gadolinium ion (8S). Salts of 
these two ions are almost inseparable by fractional crystal- 
lization; they are next to each other in the periodic table. 

The model for the structure of these electrolytes which 
the spectra seem to place before us has little resemblance 
to the model evolved by the interionic attraction theory. 
There is, for example, no sign of a diffuse ionic atmosphere 
changing continuously in extent and field intensity with 
dilution. The spectra give no suggestion that at greater 
dilution discrete structures will disappear. The fact is that 
when one discrete structure disappears another as discrete 
and fixed takes its place. 

There may, of course, be somewhat similar phenomena 
with bivalent ions but there seems little likelihood that 
they are prominent in univalent ions where the Debye- 
Hiickel theory has been so conspicuously successful. 
However, with simple trivalent ions in the range of con- 
centrations we have investigated, there are, it seems, 
sufficient reasons for expecting the Debye-Hiickel theory 
to be in disagreement with the experimental results. 


FREED 


G. H. Jones Laboratory, 
University of Chicago, 
Chicago, Illinois, 

September 9, 1938. 


1S, Freed and S. I. Weissman, J. Chem. Phys. (1938). 

2S. Freed and H. F. Jacobson, J. Chem. Phys., in this issue. 

3H. Bethe, Ann. der Physik 3, 133 (1929). 

4 The osmotic coefficients of the salts of those rare earths which have 
been measured (C. M. Mason, J. Am. Chem. Soc. 60, 1638 (1938)) 
differ extremely little. The isomorphism of the crystals, their miscibility 
with each other, the slight differences in the ionic sizes and in the molal 
volumes and especially the great difficulty of separating them by 
fractional crystallization. 

5 International Critical Tables. 

6 The exceptionally sharp lines, that is, the exceptionally weak 
coupling which Eu*** has with the oscillations, rotations, etc., of the 
environment, may, perhaps, be related to the zero value of the angular 
momentum of the electronically activated state, i.e., J=0. The basic 
state of Eu*** is actually J =0 and the thermally activated states are 
J =1 etc. Transitions between these and an activated state J =O are 
very probable. 


The Separation of the Carbon Isotopes by Diffusion 


The results to be described were obtained with a battery 
of 51 Hertz-type pumps as modified by Professor W. 


’ Bleakney of Princeton.! Preliminary experiments on mix- 


tures of argon-helium, carbon dioxide-nitrogen and carbon 
dioxide-argon permitted the determination, for these mix- 
tures, of the best operating conditions in respect to applied 
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heat and working pressure. With respect to the latter, our 
results confirm those of Barwich.? As for the question of 
the presence or absence of a capillary between the pump 
units we have found that the yield is better without 
a capillary. For the two mixtures carbon dioxide-argon 
and carbon dioxide-nitrogen we have obtained separation 
factors of 2.4 and 5.4 respectively. From these values it is 
possible to calculate, for methane, a theoretical separation 
factor of 1.18. Experimentally we found that, for methane, 
better separation of the carbon isotopes could be obtained 
with higher heats than in the cases of the above mixtures, 
In spite of this, however, the experimental separation 
factor is less than the theoretical value, in agreement with 
the data of Sherr.' At the present time, with an adequate 
system of purification and circulating a large volume of 
methane, we have obtained in thirty hours of pumping a 
yield of 300 cc of methane at 1.8 mm Hg pressure and 
containing from 30 to 32 percent of carbon 13. This is 
indicated by the relative intensities of the Swan bands 
C8—C and reproduced in Fig. 1. 

We wish particularly to thank Messrs. Breyre and 
Coppens of the Institut National des Mines at Paturages, 
who placed at our disposal specially purified methane. 

P. CAPRON 

J. M. DELFossE 
M. DE HEMPTINNE 
H. S. TayLor 


Institut de Physique, 
Université de Louvain, 
Louvain, Belgium, 

September 15, 1938. 


1R. Sherr, J. Chem. Phys. 6, 251 (1938). 
2 Barwich, Zeits. f. Physik 100, 166 (1936). 


The Diamagnetism of Gaseous Nitrosyl Chloride* 


In a recent paper Jahn! suggested that nitrosyl chloride 
(NOCI) exists in a triplet normal state, perhaps with a 
low-lying singlet state, in order to explain the apparent 
excess entropy of R In 3 or R In 4 calculated from equi- 
librium measurements over that calculated from struc- 
tural and spectroscopic information. This triplet ground 
state would lead to paramagnetism of the substance. 
Wilson? made magnetic measurements on liquid nitrosy! 
chloride and found it to be diamagnetic. There remains, 
however, the possibility that polymerization of the mole- 
cules in the liquid state might reduce the paramagnetism 
of the substance, as is the case with oxygen at sufficiently 
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low temperatures. We have therefore carried out measure- 
ments on the gaseous compound. 

The susceptibility was measured at 25° by the Gouy 
method, making comparison with oxygen as the substance 
of known paramagnetism (x= +0.142X10~ c.g.s.m. at 
20°).8 A glass tube of 19 mm diameter was separated into 
two compartments by a glass septum and equipped at 
either end for sealing off. This was suspended from the arm 
of an analytical balance, and the apparent change in weight 
(Aw in mg) of the tube was determined when the magnet 
was excited with two different currents. The absolute 
values of the magnetic fields obtained were 8080 and 9350 
oersteds at 9.00 and 14.00 ampere readings. With both 
ends evacuated there were found for the corresponding 
readings Aw=+0.10 and 0.16 (average of at least two 
determinations); with the upper end evacuated and dried 
tank oxygen at 701 mm pressure at 25° in the lower end, 
Aw=—8.29 and —11.08; with the upper end evacuated 
and pure nitrosyl chloride at 790 mm pressure at 25° in 
the lower end, Aw = +0.11 and 0.14. The probable error in 
the determination of Aw is of the order of +0.05 mg. The 
nitrosyl chloride was prepared from nitric oxide and 
chlorine which had been carefully purified for use in 
equilibrium measurements. A large excess of nitric oxide 
was allowed to react with the chlorine, and the nitrosyl 
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chloride was then frozen out in the tube and pumped free 
from volatile impurities. 

The predicted values of Aw for nitrosyl chloride in a *Z 
state at the pressure and temperature of measurement are 
—9.3 and —12.5; the predicted values for equipartition 
between three *2 and one 'Z states (quantum weight four) 
are —7.0 and —9.3; the predicted values for a singlet state 
(diamagnetic), taking into account the estimated molecular 
diamagnetism of —30X10-*, are +0.19 and 0.28 re- 
spectively. Correction for an estimated dissociation of 0.5 
percent‘ into chlorine and paramagnetic nitric oxide would 
lower the diamagnetic values to +0.15 and +0.24 re- 
spectively. The observed Aw values are in satisfactory 
agreement with this last prediction. We conclude therefore 
that nitrosyl chloride is diamagnetic and in a singlet state. 

CarRROoL M. BEESON 
CHarRLEs D, CoRYELL 
Gates and Crellin Laboratories of Chemistry, 
California Institute of Technology, 


Pasadena, California, 
August 23, 1938. 


* Contribution from the Gates and — \ oe of Chemis- 
try, California Institute of Technology, a 
1F, P. Jahn, J. Chem. Phys. 6, (193 
2 E. B. Wilson, Jr., J. Am. Chem. Soc. ~ Yond (1934). 
3 International Critical Tables, V1, 354 (19 
of unpublished of C. M. Beeson and 
. . ost. 
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